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PREFACE 


The  20TH  Annual  DoD  Atmospheric  Radiance  and  Transmittance  Meeting  was  held 
at  the  Air  Force  Phillips  Laboratory  Geophysics  Directorate  Science  Center  at  Hanscom 
AFB  from  10-12  June  1997. 

47  papers  and  presentations  were  given  including  rememberances  of  the  late 
Francis  (Frank)  X.Kneizys  who  was  one  of  the  primary  movers  in  the  Air  Force's  entry 
into  this  field. 

Updates  of  the  various  consumer  codes  such  as  MODTRAN,  MOSART  etc  and  their 
applications  in  modelling  for  sensors  and  systems  were  made.  Other  topics  addressed 
included  Spectroscopic  Techniques,  Imaging  Spectroscopy,  Aerosols,  Clouds,  Earth 
Limb,  Solar  Studies,  Climate,  Lidar  and  Microwave  techniques. 

This  volume  contains  the  abstracts,  grouped  together  at  the  front  for  easy 
perusal  and,  where  available,  hardcopy  furnished  by  the  authors.  No  significance 
should  be  accorded  to  variations  between  abstract  and  hard  copy  titles  and/or  author 
names . 

An  Author  Index  and  a  partial  list  of  Attendees  are  also  included. 

The  Editors  thank  the  authors,  and  attendees  for  their  participation. 

U  Gail  P. Anderson 

Background  Clutter  Mitigation  Branch 
Battlespace  Environment  Division 
Air  Force  Research  Laboratory 
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Research  and  Engineering. 

Col  Shaffer  was  bom  in  Dunkirk,  NY,  on  24  Jan  1955.  He  earned  a  B.  S.  in  Mathematics  from  the 
University  of  Vermont  and  Air  Force  Commission  through  St  Michael’s  College  Reserve  Officer 
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University  of  Utah,  and  subsequently  a  M.  S.  (with  Distinction)  from  the  Naval  Postgraduate  School, 
Monterey,  CA,  and  a  M.  S.  in  National  Resource  Management  from  National  Defense  University. 
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of  Technology  graduate  education  program  at  the  Naval  Postgraduate  School. 
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AFB,  OH  as  an  advanced  weather  officer  to  the  Air  Force  Foreign  Technology  Division,  where  he ' 
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foreign  weapons  systems.  In  1986,  Col  Shaffer  served  an  extended  temporary  duty  as  Officer  in 
Charge,  Weather  Support  Element,  Palmerola  Air  Base,  Honduras.  WMe  assigned  to  Det  1,  Col 
Shaffer  was  selected  as  die  1986  Air  Force  Weather  Company  Grade  Officer  of  the  Year. 

In  1987,  he  assumed  command  of  Det  2,  7th  Weather  Squadron,  Hanau  Army  Installation, 

Germany  with  concurrent  duty  as  Weather  Officer  to  the  Commanding  General,  3d  Armored  Division. 
During  this  tour,  Det  2  was  selected  as  the  1989  William  Award  winner,  given  to  the  most  outstanding 
Air  Force  Weather  unit.  In  1990,  Col  Shaffer  was  reassigned  to  5th  Weather  Wing,  Langley  Air  Force 
Base,  VA,  as  Chief,  Operational  Support  Division.  Upon  restructure  of  Air  Force  Weather,  he  became 
Chief,  Plans  and  Program  Branch,  and  Deputy  Chief,  Weather  Division,  Hq  Air  Combat  Command. 

In  1993,  Col  Shaffer  attended  the  Industrial  College  of  the  Armed  Forces,  Ft  McNair,  DC.  Upon 
graduation  in  1994,  he  was  assigned  to  the  Hq  Air  Force  Directorate  of  Weather,  Deputy  Chief  of 
Staff,  Plans  and  Operations.  While  at  Air  Staff,  he  was  Deputy  Chief,  Weather  Plans  Division;  then 
Chief,  Interagency  Division;  Chief,  Weather  Policy  Division;  and  Chief,  Weather  Resources  Division. 
He  assumed  his  present  position  and  rank  in  June  1996. 

Col  Shaffer’s  decorations  include  the  Meritorious  Service  Medal  with  du’ee  oak  leaf  clusters  and  the 
Air  Force  Commendation  Medal  with  one  oak  leaf  cluster.  He  is  married  to  the  former  Jacqueline  C 
Sklenar,  and  they  have  one  son,  Eric. 
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The  MOSART  Code:  1997 


by 


Dr.  William  M.  Comette 
National  Imagery  and  Mapping  Agency 
Bethesda,  MD  20816 


Sally  J.  Westmoreland 
Photon  Research  Associates,  Inc. 
San  Diego,  CA  92121 


Dr.  Prabhat  K.  Acharya,  Dr.  Alexander  Berk,  Dr.  David  C.  Robertson 

Spectral  Sciences,  Inc. 

Burlington,  MA  01803 

Gail  P.  Anderson,  Dr.  William  A.M.  Blumberg,  Dr.  Laila  S.  Jeong 
and  Dr.  Francis  X.  Kneizys 
Air  Force  Research  Laboratory 
Hanscom  AFB,  MA  01731 

Abstract 

The  Moderate  Spectral  Atmospheric  Radiance  and  Transmittance  (MOSART)  code 
is  a  U.S.  Government  standard  code  with  the  ability  to  support  scene  and  signature 
simulations.  MOSART  uses  the  same  band  model  and  1  cm''  resolution  band  parameter 
data  base  found  in  MODTRAN,  which  have  been  derived  from  the  1996  HITRAN  line 
atlas.  In  the  ultraviolet  and  visible  regions,  MOSART  contains  additional  molecular 
absorption  bands  at  a  coarser  resolution.  The  line-of-sight  geometries  include  all 
MODTRAN3  geometries,  plus  some  new  geometries,  including  Observer-Source, 
Observer-Background,  Observer-Source-Background,  Earth  limb.  Horizontal,  At-Source 
(i.e..  Observer  and  Source  are  co-located).  An  X-Windows-  and  Motif-based  MOSART 
Input  GUI  is  available,  as  is  a  simple  scrolling  input  file  builder. 

An  extensive  set  of  global  data  bases  is  incorporated  into  the  code,  including 
climatologies  (e.g.,  diurnal  surface  air  temperature,  three  etage  cloud  cover,  cirrus 
presence),  terrain  elevation  (10  arcmin  resolution),  water/snow  composition,  and 
ecosystem  type  (10  arcmin  resolution).  Coupled  with  the  ecosystem  and  terrain  elevations 
are  composite  terrain  scene  types  with  appropriate  boundary  layer  aerosol  composition. 
Each  scene  includes  appropriate  parameters  for  determining  the  terrain  temperature  for  each 
material  in  the  scene,  togeAer  with  the  reflected  and  emitted  radiances  for  determining  the 
mean  and  standard  deviation  for  the  scene.  The  Global  Data  Base  Model  allows  a  full  3-D 
representation  of  the  earth’s  atmosphere  (i.e.,  altitude,  latitude,  and  longitude)  with 
radiative  transfer  varying  as  the  line-of-sight  moves  within  the  spatially  changing 
atmosphere. 

The  MOSART  code  is  approximately  185,000  lines  of  ANSI  FORTRAN  77  code, 
including  about  1 1  utility  codes  and  manuals  (4  volumes:  Installation,  User's,  Technical, 
and  Software  Reference  Manuals).  Version  1.41  of  MOSART  was  released  in  September 
1995.  Since  then,  a  number  of  modifications,  upgrades,  and  "bug"  fixes  have  been  made. 
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COMPUTATIONAL  METHODS  APPLIED 
TO  THE  DEVELOPMENT  OF  HITRAN  AND  HTTEMP 


Laurence  S.  Rothman 


PL/GPOS 

Geophysics  Directorate 
Hanscom  AFB  MA  01731 


Almost  three  decades  have  passed  since  the  GOAT  (Group  On  Atmospheric 
Transmission)  team  convened  to  create  the  atmospheric  molecular  absorption  database, 
and  it  is  informative  to  review  the  direction  of  methodology  that  has  led  to  the  current 
databases,  HITRAN  and  HITEMP.  Initially,  the  techniques  of  perturbation  theory,  and 
specifically  contact  transformations,  were  applied  to  the  molecular  Hamiltonians  to 
achieve  fourth-order  accuracy  of  energy  levels.  Closed-form  solutions  were  sought,  and 
these  techniques  were  applied  especially  to  the  bent  triatomic  molecules,  H20,  03,  and 
S02.  Progress  was  accomplished  with  the  discovery  of  simplifying  generating  functions 
and  the  reduction  of  sets  of  operators.[l]  The  dipole-moment  of  water  vapor  based  on 
these  methods  is  still  the  accepted  value.  [2] 

However,  with  the  computer  power  now  available,  the  calculational  methods  of  ab 
initio  [3]  and  Direct  Numerical  Diagonalization  (DND)[4]  are  becoming  more  widely 
used  for  simple  molecular  systems.  These  methods  offer  the  capability  to  predict  spectral 
properties  not  accessible  in  controlled  laboratory  observations,  which  is  very  desirable  for 
long-path  atmospheric  transmission  and  hot-gas  problems.  In  addition,  the  computational 
methods  of  least-squares  have  played  a  vital  role  in  the  development  of  HITRAN  and 
HITEMP. 

Another  major  factor  in  HITRAN  and  HITEMP  development  is  the  accuracy  of 
laboratory  data.  Fourier  transform  spectrometer  observations  form  the  basis  for  a 
majority  of  the  line-transition  portion  of  the  compilations;  it  is  also  the  basis  of  the 
increasing  number  of  cross-section  files  on  the  compilation. 

1.  F.X.  Kneizys,  J.N.  Freedman,  and  S. A.  Clough,  J.Chem.Phys.  44, 2552 
(1966);  L.S.  Rothman  and  S.A.  Clough,  J.Chem.Phys.  54,  3246  (1971). 

2.  S.A.  Clough,  Y.  Beers,  G.P.  Klein,  and  L.S.  Rothman,  J.Chem.Phys.  59, 2254  (1973). 

3.  H.  Partridge  and  D.W.Schwenke,  accepted  J.Chem.  Phys.  (1997). 

4.  R.B.  Wattson  and  L.S.  Rothman,  J.Quant.Spectrosc.Rad.Transfer  48, 763  (1992). 


HISTORY  OF  EOSAEL 


Richard  C.  Shirkey 

Army  Research  Laboratory 
Information  Science  and  Technology  Directorate 
Battlefield  Environment  Division 
White  Sands  Missile  Range,  New  Mexico  88002-5501 


The  Electro-Optical  Systems  Atmospheric  Effects  Library  (EOSAEL)  is  a  state-of-the- 
art  computer  library  of  computer  codes  describing  various  aspects  of  electromagnetic 
propagation  and  battlefield  atmospheres  and  their  effects  on  electro-optical  sensors  and 
systems.  This  paper  presents  the  history  of  EOSAEL  and  describes  events  leading  up  to 
the  first  version,  distributed  in  1979,  to  the  current  version,  EOSAEL  92. 
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THE  APPLICATION  OF  A  VUV-FT  SPECTROMETER 
AND  SYNCHROTRON  RADIATION  SOURCE  TO 
MEASUREMENTS  OF  OXYGEN  AND  NO  MOLECULES 


K.  Yoshino,  J.R.  Esmond,  and  W.H.  Parkinson 

Harvard  Smithsonian  Center  for  Astrophysics 
Cambridge,  MA  02138,  U.S.A. 

A.P.  Thorne,  J.E.  Murray,  G.  Cox,  and  R.C.M.  Learner 
Imperial  College  of  Science,  Technology,  and  Medicine,  London,  U.K. 

K.  Ito,  T.  Imajo,  and  T.  Matsui 
Photon  Factory,  KEK,  Tsukuba,  Ibaraki  305,  Japan 

A.S-C.  Cheung  and  K.W-S.  Leung 
University  of  Hong  Kong,  Hong  Kong 


We  moved  the  IC~VUV  FT  spectrometer  from  Imperial  College,  London  to  the 
Photon  Factory,  Japan  to  exploit  the  bandwidth-limited  synchrotron  radiation  as  a 
background  source  for  FT  absorption  spectroscopy.  We  have  measured  all  absorption 
bands  of  NO  in  the  wavelength  region  195  to  160  nm  with  the  resolution  of  0.06  cm  ' 
(about  a  half  of  the  Doppler  widths).  With  resolution  of  0.06  or  0.12  \wn,  we  extended 
the  measurements  of  the  Schumann-Runge  bands  of  02  with  v'  >  12  in  the  wavelength 
range  181  to  175  nm. 

We  will  present  some  preliminary  results  of  \o2  and  NO  measurements. 

This  work  is  partly  supported  by  NSF  Division  of  Atmospheric  Sciences  grant  ATM- 
94-22854  to  Harvard  College  Observatory,  and  by  NASA  Upper  Atmospheric  Research 
Program  under  Grant  No.  NAG5-484  to  the  Smithsonian  Astrophysical  Observatory. 

We  also  acknowledge  the  support  of  the  UK  Science  and  Engineering  Research  Council 
and  the  Paul  Instrument  Fund  of  the  Royal  Society  for  the  development  of  the  UV-FT 
spectrometer. 
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02  PRESSURE  INDUCED  CONTINUUM  IN  THE  IR 


Eli  J.  Mlawer 

Atmospheric  and  Environmental  Research,  Ind. 


A  recent  high-resolution  measurement  of  surface  solar  radiance  show  clear  indication 
of  continuous  absorption  due  to  the  three  strongest  02  X3Sg  ->  a3Dg  transitions.  The 
differences  between  these  measurements  and  a  line-by-line  calculation  were  used  to 
determine  the  properties  of  these  continua  and,  for  two  of  the  bands,  led  to 
parameterizations  of  the  spectral  behavior  of  the  absorption  coefficients.  The  results 
indicate  that  these  continuous  bands  remove  approximately  0.8  W/m2  from  the  solar 
beam. 
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ACCELERATED  LINE-BY-LINE  CALCULATION  OF  ABSORPTION 
COEFFICIENTS  WITH  HIGH  NUMERICAL  ACCURACY 


Martin  Kuntz 

Forschungszentrum  Karlsruhe  GmbH 
IMK  Karlsruhe,  Germany 


An  efficient  method  for  the  calculation  of  spectral  absorption  coefficients  on  a  line-by¬ 
line  basis  is  presented.  For  each  spectral  line  retained  in  the  model,  this  method  evaluates 
a  contribution  to  the  absorption  coefficient  to  a  predefined  absolute  level  of  accuracy.  The 
algorithm  requires  the  definition  of  a  sequence  of  successively  finer  grids  and 
interpolation  from  coarser  to  finer  grids.  An  advantage  of  the  method  is  that  lines  do  not 
necessarily  contribute  over  the  entire  frequency  range  of  interest.  In  addition  the 
treatment  of  the  Voigt  line  shape  is  put  onto  a  more  rigorous  basis.  The  absorption 
coefficients  are  finally  interpolated  onto  a  non-equidistant  frequency  grid  which  offers  the 
prospect  of  improved  performance  modeling  the  radiative  transfer. 
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IMPROVED  ACCURACY  INFRARED  ABSORPTION  DETERMINATIONS 
USING  THE  MODIFIED  KUBELKA  MUNK  MODEL 


W.  G.  Egan,  and  Q.  Alexis 

Natural  Sciences  Department 
York  College/CUNY 
Jamaica,  NY  11451 


In  order  to  further  improve  the  accuracy  of  infrared  aerosol  optical  properties 
determinations,  a  modification  to  the  original  Kubelka-Munk  scattering  model  was  made 
to  include  the  effects  of  sample  surface  reflection  on  the  determination  of  absorption 
(i.e.the  Modified  Kubelka  Munk  Model:  M-K-M).  The  solutions  to  the  radiative  transfer 
equations  involved  for  the  M-K-M  model  are  presented  graphically.  The  results  indicate 
that  the  M-K-M  model  generally  produces  somewhat  less  absorption  than  the  K-M 
model,  but  comparable  to  the  dispersion  model  in  high  absorption  regions.  A  plot  of 
absorption  and  scattering  in  the  complex  plane  reveals  problems  with  the  dispersion 
theory  when  significant  scattering  occurs. 
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FASE/FASCOD4:  CURRENT  STATUS  AND  FUTURE  DIRECTION 


H.E.  Snell 


Atmospheric  and  Environmental  Research,  Inc. 
840  Memorial  Dr. 

Cambridge,  MA  02139 

G.P.  Anderson,  J.H.  Chetwynd,  and  S.  Miller 


Geophysics  Directorate 
Phillips  Laboratory 
Hanscom  AFB,  MA  01731 


The  Optical  Physics  Division  of  the  Air  Force  Phillips  Laboratory  with  support  from 
the  Department  of  Energy  (DoE)  Atmospheric  Radiation  Measurement  (ARM)  Program 
is  developing  a  state-of-the-art  line-by-lineatmospheric  radiative  transfer  model  as  the 
successor  to  FASCODE  (Fast  Atmospheric  Transmittance  Code).  The  goal  of  this 
project  is  to  create  a  computationally  efficient  model  which  contains  the  most  up-to-date 
physics.  The  new  model,  known  as  FASCODE  for  the  Environment,  or  "FASE", 
combines  the  best  features  of  FASCODE  and  LBLRTM  (Line-by-Line  Radiative 
Transfer  Model),  the  DoE's  standard  radiative  transfer  model.  Upgrades  to  FASE  include 
the  addition  of  a  solar  spectrum  module  to  compute  the  attenuated,  line-of-sight  solar 
radiance,  improvement  of  the  cloud  and  aerosol  descriptors,  based  on  changes  made  to 
MODTRAN  (Moderate  Resolution  Transmittance  Model),  and  the  ability  to  incorporate 
the  new  heavy  molecule  absorption  cross-section  data  found  on  the  HrrRAN96  database, 
which  is  of  slightly  different  format  from  the  previous  HTTRAN  cross-section  data.  This 
paper  addresses  changes  which  have  been  made  to  FASCODE  and  LBLRTM  to  create 
FASE,  gives  an  overview  of  the  new  capabilities,  and  outlines  current  and  future 
development  plans. 


FASE-CALCULATED  CONTRIBUTION  AND  SENSITIVITY  FUNCTIONS  FOR 
MSX  SPIRIT  m  MWIR  BANDS  FOR  BTH  MEASUREMENTS 


W.  Gallery,  H.  Snell,  C.  Sarkiasian,  D.  Hogan 
AER,  Inc. 

A.  Zachor 

Atmospheric  Radiation  Consultants,  Inc. 

J.  Winick 

Phillips  Laboratory 
Geophysics  Directorate 


The  radiance  measured  at  a  satellite  viewing  the  earth  originates  from  the  surface  and 
the  various  levels  of  the  atmosphere.  The  "contribution  function"  describes  the  relative 
contribution  each  level  makes  to  the  outgoing  radiance  while  the  "sensitivity  function" 
describes  the  sensitivity  of  the  outgoing  radiance  to  changes  in  an  atmospheric  parameter 
(e.g.,  temperature  or  water  vapor)  at  each  level.  Both  functions  are  important  in 
understanding  the  variations  in  the  measured  radiance  along  the  satellite  path.  In  this 
paper,  both  functions  are  defined  rigorously.  FASE  is  used  to  calculate  both  functions 
efficiently  for  the  MSX  SPIRIT  HI  MWIR  bands  (both  the  radiometer  bands  B 1  and  B2 
and  the  interferometer  channel  2.)  These  functions  are  presented  for  various  angles  from 
nadir  out  to  the  limb  and  for  both  LTE  and  Non-LTE  (contribution  function  only)  cases. 
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FASE,  FASCODE  AND  RADTRAN  ATTENUATION 
DIFFERENCES  IN  THE  MICROWAVE 


Michael  K.  Griffin  and  Alan  E.  Lipton 

Optical  Effects  Division 
Phillips  Laboratory 

Hanscom  Air  Force  Base,  Massachusetts,  01731-3010 

Alan  G.  Ling 

Radex,  Inc. 

3  Preston  Ct 

Bedford,  Massachusetts,  01730 


Modeling  radiative  transfer  in  the  microwave  has  become  increasingly  important  with 
the  wide  usage  of  the  DMSP  SSM/I  and  SSM/T-2  microwave  sensors.  Current  on-orbit 
sensor  frequencies  range  from  19  GHz  to  183  GHz;  experimental  airborne 
instrumentation  will  include  frequencies  up  to  650  GHZ.  The  increasing  interest  in  the 
full  spectrum  of  microwave  frequencies  and  the  selection  of  radiative  transfer  models 
available  to  the  user  prompted  this  comparative  study.  The  basic  output  (attenuation) 
from  three  models,  RADTRAN,  FASCODE  and  FASE  were  compared  over  the  range  of 
frequencies  described  above.  Differences  were  noted  and  the  effects  these  variations  had 
on  our  modeling  efforts  were  documented.  Specifically,  how  these  differences  might 
affect  our  study  of  the  sensitivity  of  computed  microwave  brightness  temperature  to  cloud 
liquid  water  content.  This  presentation  will  highlight  the  model  differences  and  discuss 
how  these  model-to-model  variations  might  affect  specific  parameter  estimations. 
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CHARACTERIZATION  OF  OPTICAL  EXTINCTION  PROFILES 
USING  RAMAN  LIDAR  TECHNIQUES 


C.  Russell  Philbrick  and  Michael  O'Brien 


The  Pennsylvania  State  University 
Electrical  Engineering  Department 
University  Park,  PA  16802 


Profiles  of  Raman  scattering,  due  to  vibrational  and  rotational  transitions,  from  a  the 
laser  transmitter  of  a  lidar  system  provide  a  signal  which  can  be  directly  interpreted  as 
optical  extinction.  We  have  been  using  the  Raman  scatter  from  a  Nd:YAG  laser  beam 
from  the  2nd,  3rd  and  4th  harmonics  at  532,  355  or  266  nm  to  obtain  profiles  of  the 
optical  extinction.  The  measurements  have  been  made  using  the  molecular  nitrogen 
profiles  at  607,  387  and  285  nm,  which  are  the  first  Stokes  vibrational  Raman  shifts  from 
the  Nd:YAG  laser.  The  rotational  Raman  scattered  radiation  at  530  nm  has  also  been 
used  to  obtain  optical  extinction  profiles.  The  vertical  profiles  of  optical  extinction  have 
been  measured  during  several  time  sequences  which  allow  a  study  of  the  evolution  of 
cloud  development.  The  growth,  dissipation  and  advection  of  the  cloud  particles  has 
been  investigated  using  the  technique.  The  fact  that  the  profiles  are  simultaneously 
measured  at  several  wavelengths  allows  some  inference  of  the  distribution  of  particle  size 
distribution  within  a  cloud  layer.  Significant  variations  in  the  relative  sizes  of  the 
particles  have  been  observed  between  the  upper  and  lower  regions  of  a  cloud  layer. 
Example  results  are  used  to  show  several  of  the  features  which  have  been  observed 
during  measurement  periods.  The  lidar  instruments  which  we  are  using  to  obtain  the 
optical  extinction  profiles  are  simultaneously  measuring  the  meteorological  profiles  of 
the  water  vapor  and  temperamre.  The  data  permit  investigations  of  the  physical  process 
involved  in  the  evolution  of  cloud  optical  properties. 
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MULTI-WAVELENGTH  /  TUNABLE  LASER 
RADAR  SIMULATION  CODE  LIDAR-PC 


William  E.  Wilcox 


W.  J.  Schafer  &  Associates 
321  Billerica  Road 
Chelmsford,  MA  08124-4191 


Dennis  K.  Kiliinger 

University  Of  South  Florida 
4202  E.  Fowler  Ave. 
Tampa,  FL  33620 


We  have  developed  a  Laser  Radar  (LIDAR)  simulation  code  which  models  multi¬ 
wavelength  or  tunable  system  performance  over  slant  paths  in  the  atmosphere.  Both  hard 
target  and  aerosol  LIDAR  systems  may  be  modeled  with  the  software. 

The  program  combines  high  resolution  spectral  data  from  HITRAN-PC  with  Aerosol 
models  from  LOWTRAN  7  to  predict  LIDAR  S/N  as  function  of  wavelength  and 
altitude/range.  Various  LIDAR  system  parameters  such  as  power,  wavelength,  telescope, 
and  detector  noise  are  entered  into  a  system  worksheet.  Calculations  are  then  made  based 
on  the  worksheet  and  plotted  in  an  interactive  Windows  environment. 


13 


TEMPERATURE  AND  TRACE  GAS  RETRIEVALS 
FROM  THE  CIRRIS  1-A  STRATOSPHERIC  LIMB  DATA 
UTILIZING  THE  OPTIMAL  ESTIMATION  METHOD 


Steven  M.  Miller  and  Jeremy  Winick 

Optical  Effects  Division 
USAF  Phillips  Laboratory 

Jean-Luc  Moncet  and  Hilary  E.  Snell 

Atmospheric  and  Environmental  Research,  Inc. 


A  general  purpose  retrieval  code  has  been  developed  utilizing  a  derivative  of 
FAS  CODE  as  the  forward  model  and  Rodger’s  optimal  estimation  method  as  the 
inversion  algorithm.  This  code  has  been  applied  to  stratospheric  limb  data  taken  with  the 
CIRRIS-IA  interferometer  flown  aboard  the  space  shuttle  in  April  1991.  The  output  of 
the  code  provides  simultaneously  retrieved  temperature  and  volume  mixing  ratio  profiles 
of  several  important  trace  gases  and,  as  important,  the  error  and  the  relative  contributions 
from  the  measurement  and  the  first  guess  profile  as  background  information.  The 
observation  of  non  local  thermodynamic  equilibrium  in  the  laser  band  of  C02  is  modeled 
using  the  Atmospheric  Radiance  Code,  ARC,  and  compensated  for  in  the  data  using  the 
NLTE  algorithms  in  FASCODE. 
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SCIAMACHY  LIMB  MODES,  PRESENTATION 
AND  PRELIMINARY  SIMULATIONS 


Christian  Muller',  John  Burrows*  and  Albert  Goede* 
'Belgian  Institute  for  Space  Aeronomy,  Brussels. 

*  IFE,  University  of  Bremen. 

I 

*Stichting  Ruimte  Onderzoek  Nederland,  Utrecht. 


SCIAMACHY  is  a  joint  project  of  Germany  ,the  Netherlands  and  Belgium  for  global 
atmospheric  measurements.  It  has  been  selected  by  ESA  to  fly  on  the  Envisat-1  polar 
platform,  to  be  launched  in  1998/1999.  SCIAMACHY  comprises  a  moderately  high 
resolution  (0.2-0.5  nm)  spectrometer  to  observe  transmitted,  reflected  and  scattered  light 
from  the  atmosphere  in  the  UV,  visible  and  near  infrared  wavelength  regions  over  the 
range  240-1700  nm,  and  in  2  selected  regions  between  2.0  and  2.4  micron.  The  goal  is  to 
allow  small  optical  absorptions  (as  small  as  2E-4  in  some  regions  of  the  spectrum)  to  be 
detected. 

While  the  SCIAMACHY  instmment  is  mainly  known  for  its  nadir  mode  where  it  will 
improve  on  the  existing  GOME  (Global  Ozone  Measuring  Experiment),  it  has  also  been 
designed  with  a  sophisticated  limb  mode  including  solar  occultation,  lunar  occultation 
and  measurement  of  limb  radiances  between  100  km  of  altitude  and  the  surface.  This 
mode  has  no  previous  equivalent  in  space-borne  earth  observation  experiments. 
Requirements  for  the  direct  models  allowing  to  simulate  this  limb  mode  will  be 
presented. 

Simulations  of  the  infrared  channels  occultation  modes  will  also  be  shown.  The 
spectral  intervals  were  chosen  essentially  for  nadir  sounding  of  tropospheric  columns  of 
N20,  CH4  and  CO  with  the  limitation  of  using  state  of  the  art  non  cryogenic  detectors. 
However,  simulations  show  that  occultation  data  will  be  relevant  for  the  determination  of 
the  lower  stratosphere  vertical  distribution  of  these  gases  and  water  vapor.  This  property 
makes  it  a  new  tool  for  the  understanding  of  the  troposphere-stratosphere  exchange 
mechanisms. 
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ADEOS/ILAS  MEASUREMENT  OF  THE  POLAR  STRATOSPHERE 


Tatsuya  Yokota 

National  Institute  for  Environmental  Studies 
16-2  Onogawa,  Tsukuba,  Ibaraki  305  Japan 


The  Advanced  Earth  Observing  Satellite  (ADEOS)  has  been  successfully  launched  on 
August  17,  1996.  Improved  Limb  Atmospheric  Spectrometer  (ILAS)  is  a  grating 
spectrometer  aboard  ADEOS.  BLAS  is  designed  to  measure  stratospheric  minor 
constituents  and  parameters  with  a  solar  occultation  technique  around  the  polar  region 
for  both  hemispheres. 

Vertical  profiles  of  temperature,  pressure,  and  aerosol  extinction  coefficient  are 
determined  from  the  visible  channel  covering  the  753  ~  784  nm  for  02  molecular 
absorption.  ILAS  has  also  observed  spectral  feature  of  ozone  Wulf  band  absorption  in 
this  spectral  range.  To  estimate  the  Wulf  band  contribution  to  the  transmittance  spectra, 
cross-section  data  from  MODTRAN3. 5-beta  are  used  in  the  calculation. 

Profiles  of  03,  CH4,  N20,  HN03,  H20,  N02,  CFC-11,  CFC-12,  and  aerosols  are 
simultaneously  obtained  from  the  IR  channel  (6.2  ~  1 1.8  mm).  In  the  ILAS  data  retrieval 
code,  HTTRAN  database  is  used  for  the  absorption  cross  section  calculation.  02 
continuum  calculation  around  6  mm  has  referred  to  the  method  in  the  FASCODE.  The 
method  of  CKD2. 1  from  Dr.  S.  A.  Clough  is  used  for  the  H20  continuum  calculation.  As 
for  the  cross-section  calculation  of  CFC-11,  CFC-12,  and  N205,  pseudo-line  parameters 
from  Dr.  G.  C.  Toon  are  used. 

ILAS  has  operationally  observed  polar  stratosphere  from  late-November  1996. 
Preliminary  results  of  the  ILAS  measurement  will  be  presented  here. 
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CIRRIS-IA  LIMB  SPECTRAL  MEASUREMENTS  ON  MESOSPHERIC  OZONE 


D.  K.  Zhou,'  M.  G.  Mlynczak,'  G.  E.  Bingham,' 
J.  O.  Wise,"  AND  R.  M.  Nadile" 

'Space  Dynamics  Laboratory 
Utah  State  University 
Logan,  Utah 

^Atmospheric  Sciences  Division 
NASA  Langley  Research  Center 
Hampton 

^Geophysics  Directorate 
Phillips  Laboratory 
Hanscom  AFB,  Massachusetts 


Infrared  spectral  radiance  measurements  were  made  by  the  Cryogenic  Infrared 
Radiance  Instrumentation  for  Shuttle  (CIRRIS  lA)  during  one  flight  of  the  space  shuttle 
Discovery  (28  April  to  6  May,  1991).  Observations  were  made  for  both  the  nighttime 
and  daytime.  The  radiance  from  the  ozone  9.6  mm  band  region  shows  a  significant 
difference  between  night  and  day  in  the  mesosphere.  Band-to-band  and  scan-to-scan 
spectral  analyses  (in  the  wavenumber  region  of  1010  -  1140  cm-1)  were  made  using  the 
Strategic  High-Altitude  Radiance  Code  (SHARC)  model  to  separate  the  contributions 
from  ozone  fundamental  bands,  hot  bands,  and  the  C02  band  in  this  spectral  region.  A 
spectrally  integrated  ozone  fundamental  u3  band  radiance  profile  was  used  to  retrieve  its 
excited  state  density.  The  excited  state  density  in  conjunction  with  the  modeled 
vibrational  temperature  was  used  to  estimate  the  total  ozone  density.  The  ozone 
abundance  derived  from  the  SHARC  model  and  the  Simplified  Single  Mode  (SSM) 
ozone  model  were  compared.  The  error  sources  and  accuracy  of  the  estimated  ozone 
abundance  were  discussed.  The  night/day  ozone  ratio  of  the  CIRRIS  lA  measurements 
demonstrates  evidence  of  a  non-LTE  mesospheric  ozone  diurnal  variation.  The 
photolysis  rate  affiliated  with  the  daytime  ozone  photo-dissociation  was  estimated.  The 
mesospheric  ozone  heating  rates  were  estimated  for  CIRRIS  lA  nighttime  and  daytime 
observations. 
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LIMB  AND  BELOW-THE-HORIZON  SPECTRA  MEASURED 
BY  THE  MSX  SPIRIT  m  INTERFEROMETER; 
COMPARISON  TO  SAMM  TEMPLATE  SPECTRA 


Alexander  S.  Zachor 

Atmospheric  Radiation  Consultants,  Inc. 

K 

R.  R.  O’Neil,  H.  A.  B.  Gardiner  and  J.  J.  Gibson 

Phillips  Laboratory 
Geophysics  Directorate 
Hanscom  AFB,  MA 


The  SPIRIT  IQ  interferometer  spectrometer  aboard  the  MSX  satellite  has  six  channels 
that  collectively  are  capable  of  measuring  terrestrial  emissions  over  the  wavelength  range 
2.6-28  mm  at  2  cm'‘  spectral  resolution.  Data  has  been  acquired  from  nadir  to  200  km 
tangent  height.  Automated  analysis  of  these  observations  includes  the  computation  of 
“template”  spectra  which  are  predictions  based  on  the  SAMM  code.  We  describe  the 
interferometer’s  capabilities  and  constmction  of  templates,  and  present  samples  of 
measured  spectra  and  comparison  template  spectra. 
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PURE  ROTATION  EMISSIONS  FROM  MESOSPHERIC  OH  OBSERVED 
USING  THE  SPIRIT  HI  INTERFEROMETER  ABOARD  THE  MSX  SATELLITE 


P.S.  Armstrong^  S.J.  Lipson*,  R.R.  0’Neil‘,  H.A.  Gardiner', 

A.  Zachor^  M.  Kendra^  J  J.  Gibson',  J.A.  Dodd^, 

D.L.  Vititoe',  W.A.M.  Blumberg',  K.  Farnham^  and  T.  Woolston® 

‘Phillips  Laboratory  Geophysics  Directorate,  Hanscom  Air  Force  Base,  MA  01731 

I 

^Stewart  Radiance  Laboratory,  139  Great  Road,  Bedford,  MA  01730 
*ARC,  Inc.,  59,  High  Street,  Acton,  MA  01720 
“Radex  Inc.,  3  Preston  Court,  Bedford,  MA  01730 
®Utah  State  University,  Logan,  UT  84322 


Pure  rotation  emissions  from  NLTE  (Non-Local  Thermodynamic  Equilibrium) 
OH(X2P)  have  been  observed  in  earthlimb-viewing  geometry  using  the  SPIRIT  III 
interferometer  aboard  the  MSX  satellite.  Emission  lines  from  rotational  levels  with  N»30 
have  been  identified  in  the  LWIR  (approximately  10-13  mm)  wavelength  region, 
representing  over  2  eV  of  rotational  energy.  These  observations  are  consistent  with 
earlier  measurements  made  during  the  CIRRIS  1 A  Space  Shuttle  experiment,  and  provide 
a  valuable  comparison  made  under  differing  aeronomic  conditions.  Simultaneous 
measurements  of  OH  fundamental  band  emissions  were  also  made  using  SPIRIT  IH,  thus 
permitting  the  characterization  of  vibrational  distributions  for  both  LTE  and  NLTE  OH 
rotational  populations.  The  calibrated  spectral  data,  obtained  using  detectors  with  vertical 
footprints  of  about  4-6  km,  permit  the  determination  of  absolute  line-of-sight  column 
densities  through  the  mesopause  layer.  These  observations  provide  important 
information  on  the  dynamics  of  OH,  a  significant  chemical  species  and  IR  emitter 
responsible  for  cooling  of  the  mesopause  region.  Results  from  a  limb  measurement  on  8 
October  1996  over  Bear  Lake,  Utah  will  be  reported.  Ground-based  OH  measurements 
were  also  recorded  at  Bear  Lake  and  Logan,  Utah. 
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OVERVIEW  OF  THE  15  MM  C02  LIMB  EMISSION  FROM  MSX 


J.O.  Wise,  R.  O’Neil,  H.  Gardiner,  J.R.  Winick,  and  R.H.  Picard 

Geophysics  Directorate 
Phillips  Laboratory 
Hanscom  AFB,  MA 

P.P.  Wintersteiner 
Arcon  Corporation 
Waltham,  MA 

M.  Kendra 
Radex,  Inc. 

3  Preston  Ct 

Bedford,  Massachusetts,  01730 


Sample  15  mm  limb  emission  data  from  the  MSX  SPIRIT  HI  sensor  is  presented  for 
tangent  heights  from  65  to  150  km.  The  data  were  measured  by  band  D  (13.5-15.9  mm 
FWHM)  of  the  radiometer.  Band  D  consists  of  4  columns  of  192  pixels  each.  The 
columns  are  oriented  vertically  in  the  limb  viewing  measurements  and  extend  over  an 
altitude  range  of  60  km  in  a  single  measurement.  Spatial  resolution  of  300  m  is  achieved 
in  the  limb  views.  The  spacecraft  has  superb  pointing  capabilities  and  nearly  constant 
tangent  height  altitudes  are  maintained  during  the  measurements,  nominally  30  minutes 
in  duration.  The  data  were  collected  for  various  geophysical  conditions  during  1996  and 
1997.  The  emissions  are  examined  for  their  dependence  on  latitude  and  season,  as  well 
as  from  dissipation  of  energy  from  tidal  and  gravity  waves.  Of  particular  interest  is  &e 
behavior  of  the  characteristic  knee  at  approximately  108  km.  At  this  altitude  the  emission 
can  fluctuate  significantly  due  to  variations  in  atomic  oxygen,  C02,  and  temperature,  and 
may  provide  important  clues  to  the  species  and  temperature  variability  between  100-150 
km. 
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STRATOSPHERIC  WAVE  OBSERVATIONS  FROM  MSX;  A  NEW  SOURCE 
OF  STRUCTURE  IN  43-MM  DOWNLOOKS  FROM  SPACE 


R.  H.  Picard,  E.  M.  Dewan,  R.  O’Neil,  H.AJB.  Gardiner,  J.  Gibson, 
J.  R.  Winick,  J.  Brovra,  G.  P.  Anderson, 

PhUlips  Lab  /  Geophysics3e<Mbrd,  Mass. 

E.  R.  Hegblom,  N.  A.  Grossbard,  E.  Richards, 

Boston  College  ISR,  Chestnut  Hffl,  Mass  W.  O.  Gallery, 
Atmospheric  Environmental  Research  Cambri(Jge,  Mass. 

A.  S.  Zachor,  Atmospheric  Radiation  Consultants,  Acton,  Mass. 

P.  P.  Wintersteiner,  ARCON  Corp.,  Waltham,  Mass. 

U.  B.  Makhlouf,  Stewart  Radiance  Laboratory,  Bedford,  Mass. 
J.  Gruninger,  Spectral  Sciences,  Inc.,  Burlington,  Mass. 


Much  MSX  radiometer  data  from  the  4.3-mm  spectral  region,  dominated  by  C02  n3 
band  emission  features,  shows  evidence  of  significant  stmcture  arising  in  the  cloud-free 
atmosphere.  This  occurs  in  both  below-the-horizon  (BTH)  and  above-the-horizon  (ATH) 
data-ga±ering  modes.  In  particular,  wavelike  structures  are  seen  in  these  MSX  mid-wave 
infrared  bands  in  BTH  mode  at  large  near-horizon  nadir  angles.  The  structures  are 
observed  in  both  the  narrow  B1  band  imagery  (4.22-  4.36  mm),  which  is  confined  to  the 
heart  of  the  strongest  n3  band,  and  the  wide  B2  band  imagery  (4.24-4.45  mm),  which  has 
extended  response  on  the  long-wave  side  encompassing  hot  and  isotopic  bands.  Both 
filters  have  very  steep  edges  resulting  in  high  out-of-band  rejection.  This  wavelike 
structure  is  contrasted  with  the  structure  resulting  from  clouds.  Data  taken  on  a  number 
of  day  and  night  collections  over  East  Asia  and  the  Indian  Ocean  are  exannned  for  the 
presence  of  wave  structure.  In  particular,  the  days  of  17  Aug  and  13  Nov  and  the  night 
of  2  Oct  96  are  analyzed  in  some  detail.  The  wavelike  structures  appear  to  have 
fractional  radiance  modulation  of  a  few  percent  and  scales  of  15-200  km,  consistent  with 
atmospheric  gravity  wave  (GW)  scales.  Moreover,  power  spectral  densities  of  the 
radiometer  time  series  indicate  they  obey  power  laws  consistent  with  expectations  from 
the  theory  of  saturated  GW  spectra  and  radiative  response  theory.  The  2  Oct  and  13  Nov 
show  particularly  dramatic  circular  wave  patterns,  indicating  point-source  wave 
generation  events.  High  cloud  tops  observed  in  GMS  5  satellite  LWIK.  images  at  the 
appropriate  location  and  time  iden^  specific  convective  storms  as  likely  sources.  These 
MSX  results  constitute  the  first  definitive  infirared  observations  of  GW’s  from  space  and 
implicate  such  waves  as  important  sources  of  BTH  stmcture. 
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MODTRAN  SW  AND  VNIR  RADIANCES  AND  FLUXES 


A.  Berk,  L.S.  Bernstein,  P.K.  Acharya, 

Spectral  Sciences  Inc. 

99  S.  Bedford  Street,  Suite  No.  7 
Burlington,  MA  01803 

G.P.  Anderson 


Geophysics  Directorate 
Phillips  Laboratory 
Hanscom  AFB,  Massachusetts 


The  Correlated-k  (CK)  option  of  MMODTRAN4  partitions  band  mode  transmittances 
into  Beer's  Law  components  for  proper  implementation  of  multiple  scattering  algorithms. 
Comparisons  between  CK  and  standard  band  model  short-wave  (SW)  and  visible/near- 
infrared  (VNIR)  spectral  fluxes  at  the  surface,  the  tropopause  and  the  top  of  the 
atmosphere  will  be  presented  to  benchmark  the  level  of  error  introduced  by  neglect  of 
Beer's  Law.  Results  will  be  presented  at  two  calculation  step  sizes,  for  clear  and  cloudy 
skies,  and  using  simple  2-stream  and  discrete  ordinate  (DISORT)  N-stream  multiple 
scattering  methods.  An  upgrade  to  the  MODTRAN  band  model,  important  for  02  band 
at  0.76|im,  will  also  be  presented. 
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RECENT  UPGRADES  TO  MODTRAN  -  NEW  SOLAR  IRRADIANCE 
DATABASES,  INCLUSION  OF  NOVAM,  VERSATILE  FREQUENCY  INPUTS 
AND  INCLUSION  OF  INSTRUMENT  SCANNING  FUNCTIONS 


P.  Acharya,  A.  Berk,  L.  Bernstein  and  D.  Robertson 

Spectral  Sciences  Inc. 

99  S.  Bedford  Street,  Suite  No.  7 
Burlington,  MA  01803 

G.P.  Anderson  and  J.  Chetwynd 

Phillips  Laboratory/GPO 
Hanscom  AFB,  MA  01731 


Recent  upgrades  are  intended  to  make  MODTRAN  a  more  accurate,  convenient  and 
versatile  atmospheric  modeling  tool.  The  current  default  solar  irradiance  database  due  to 
Kurucz  was  corrected  in  the  310-340  nm  region  using  measured  data.  Several  other 
datasets  are  also  made  available  to  the  user.  The  user  can  also  use  his/her  own  data 
residing  in  a  file.  Examination  confirms  that  the  corrected  Kurucz  data  is  still  the  best 
overall  among  all  currently  available  datasets.  MODTRAN  can  now  be  run  using 
frequency  inputs  which  are  in  wavenumber,  nm  or  micron.  Additionally,  there  is  now  the 
capability  of  obtaining  radiance  and  transmittance  which  are  convolved  with  an 
appropriate  instrument  scanning  function.  Several  choices  of  scanning  functions  are 
available  including  an  user-defined  option.  This  upgrade  is  helpful  to  users  who  would 
like  to  model  satellite  data  which  are  gathered  as  a  function  of  wavelength  and  have 
undergone  instrument  degradation.  All  upgrades  are  accomplished  while  maintaining 
strict  compatibility  of  TAPES  with  earlier  versions  of  the  code.  Work  is  currently  in 
progress  to  incorporate  NOVAM  (Navy  Oceanic  Vertical  Aerosol  Models)  into 
MODTRAN. 
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SOLAR  ULTRAVIOLET  IRRADIANCE  AND  ITS  VARIATION 


Giuiiana  de  Toma,  Gary  J.  Rottman,  Thomas  N.  Woods 

Laboratory  for  Atmospheric  and  Space  Physics 
Boulder,  Colorado 


Solar  ultraviolet  irradiance,  especially  at  wavelengths  below  300  nm,  is  strongly 
absorbed  by  the  Earth's  middle  atmosphere,  where  is  the  dominant  source  of  energy. 
Knowledge  of  the  amount  of  solar  energy  deposited  in  the  atmosphere  and  its  variation 
with  time  is  important  for  atmospheric  studies. 

The  Solar  Stellar  Irradiance  Comparison  Experiment,  SOLSTICE,  on  the  UARS 
satellite  makes  daily  observations  of  solar  spectral  irradiance  in  the  wavelength  region 
from  120  to  420  nm.  The  UARS  observations  start  in  late  1991,  near  the  peak  of  solar 
cycle  22,  and  continue  to  the  present  solar  minimum  condition.  During  this  period  two 
types  of  variation  are  most  notable  in  SOLSTICE  data:  the  27-day  variations  related  to 
the  rotation  period  of  theSun  and  the  longer  term  variability  associated  with  the  11 -year 
solar  cycle. 

In  this  paper,  the  SOLSTICE  irradiance  measurements  are  reviewed  with  particular 
emphasis  on  variability.  For  the  period  1991-1996,  we  report  a  variation  of  about  a  factor 
of  two  at  Lyman  alpha  and  a  variation  of  about  6%  near  200  nm,  gradually  decreasing  to 
less  than  1%  near  300  nm. 
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COMPARISON  OF  RSS  MEASUREMENTS  WITH  MODTRAN 


Lee  Harrison,  QELONG  MIN,  and  Joe  Michalsky 

Atmospheric  Sciences  Research  Center 
State  University  of  New  York,  Albany 


The  Rotating  Shadowband  Spectroradiometer  (RSS)  is  a  new  instrument  that  provides 
spectral  direct,  diffuse,  and  tot^  horizontal  irradiances  with  a  512  CCD-arxay  over  the 
wavelength  range  from  350  nm  to  1075  nm.  We  compare  the  measured  spectra  from  the 
RSS  to  that  t^en  by  collocated  MFRSR  and  MODTRAN-3.5  calculations.  The 
measurements  from  both  instruments  (RSS  and  MFRSR)  agree  within  2%.  The  model 
calculated  direct-horizontal  irradiances  show  very  good  agreement  with  the 
measurements.  However,  the  model  substantially  overpredicts  the  diffuse-horizontal  if 
the  standard  "rural  aerosol"  model  is  applied,  instead  we  use  MODTRAN's  "urban 
aerosol"  properties  for  the  asymmetric  factor  and  single  scattering  albedo,  the  discrepancy 
can  be  eliminated.  The  spectral  signatures  of  discrepancy  suggest  that  the  scattering 
treatment  of  the  diffuse  irradiance  and  uncertainties  about  aerosol  microphysics  might  be 
the  potential  reasons  for  "clear-sky  anomaly". 
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SOLAR  SPECTRAL  MEASUREMENTS:  A  UNIQUE  DATASET 


ALEXANDER  F.  H.  GOETZ‘,*  Kathleen  B.  Heidebrecht and  Bruce  Kindel ' 

'  Center  for  the  Study  of  Earth  From  Space/CIRES 

^  Department  of  Geological  Sciences 
University  of  Colorado 
Boulder,  CO  80309 


A  database  has  been  created  containing  high  temporal  (one  per  minute)  solar 
irradiance  and  radiance  measurements  taken  in  Boulder,  CO  during  September  and 
October,  1996  and  Petersham,  MA  in  June,  1996.  The  measurements  were  made  under  a 
variety  of  sky  conditions.  The  data  were  measured  with  an  ASD-FR  instrument  that  has  a 
spectral  range  from  .4-2.5  microns,  resolution  of  3-10  nm,  and  a  sampling  interval  of  2 
nm.  The  data  quality  is  good  and  is  in  high  agreement  with  MODTRAN  simulations.  A 
large  effort  was  made  to  get  an  absolute  calibration  for  the  instrument  and  to  validate  the 
data  against  independent  measurements.  Digital  photographs  of  the  sky  will  be  taken 
coincident  with  the  spectral  measurements  planned  for  the  spring  and  summer  of  1997. 
The  data  will  be  stored  in  a  commercial  database  and  available  to  other  researchers. 
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MWIR  ATMOSPHERIC  CORRECTION  ALGORITHM 


Hsiao-hua  K.  Burke,  J.  William  Snow,  Kristine  E.  Rhoades 
John  Kerekes,  Mike  P.  Jordan  and  Carolyn  Upham 


MIT  Lincoln  Laboratory 
244  Wood  Street 
Lexington,  MA  02173 


Currently  there  are  mid-wavelength  infrared  (MWIR)  space  assets  being  used  for 
surface  characterization.  An  effort  was  undertaken  to  reconstruct  surface  emitted  MWIR 
radiance  from  satellite  measurements  using  readily  available  atmospheric  information.  A 
new  surface  radiance  retrieval  algorithm  based  on  MWIR  satellite  data  (GOES  Imager 
Channel  2  alone)  was  developed.  It  is  demonstrated  that  retrieval  accuracy  can  be 
improved  by  using  existing  NOAA  operational  products,  specifically  the  'Eta'  prediction 
model  output,  combinedwith  the  state-of-science  atmospheric  radiative  transfer  model 
(MODTRAN).  Various  test  cases  were  investigated.  Significantly,  the  retrieval  of  sea 
surface  temperature  over  ocean  at  nighttime  is  shown  to  have  a  bias  of  0. 1  K  and  rms  of 
0.5  K.  In  this  presentation,  additional  cases  of  land  surfaces  as  well  as  daytime 
conditions  are  also  illustrated.  Analysis  of  environmental  parameter  sensitivities  and 
comparisons  of  results  between  LWIR  and  MWIR  will  be  discussed. 
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APPLICATION  OF  MODTRAN  TO  SEA  SURFACE 
TEMPERATURE  RETRIEVAL  SIMULATION 


D.  Hogan,  W.  Gallery 
AER,  Inc. 


Sea  surface  temperatures  (SST)  are  traditionally  measured  from  satellites  using 
radiometer  bands  in  the  8  to  12  and  the  3  to  5  micrometer  window  regions.  In-situ 
measurements  of  SST's  are  correlated  with  simultaneous  satellite  radiance  measurements 
to  develop  regression  equations  relating  the  measured  radiances  to  the  SST.  These 
regression  equations  are  then  applied  to  routine  satellite  radiance  measurements  to  obtain 
operational  SST's. 

We  have  used  MODTRAN  to  simulate  the  radiance  for  a  proposed  satellite 
radiometer.  A  dataset  of  globally  representative  atmospheric  profiles  of  temperature  and 
water  vapor  (the  TIGR  dataset)  was  split  into  training  and  validation  datasets.  The 
training  set  was  used  with  MODTRAN  to  calculate  the  simulated  radiances  and  the 
coefficients  for  several  multi-linear  regression  models.  The  validation  set  was  used  to 
derive  the  error  statistics  for  the  regression  models.  Several  models  were  tested, 
including  various  linear  and  non-linear  terms.  Finally,  the  models  were  tested  against  a 
simulated  scene  calculated  using  an  independent  model.  Results  of  these  simulations  will 
be  presented. 
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USE  OF  MODTRAN  IN  PERFORMANCE  AND  ANALYSIS  CAPABILITY  FOR 
EARTH  OBSERVATION  SYSTEMS  (PACEOS)  SIMULATION  SOFTWARE 


Steve  Marusa  and  Rafael  Rivera 
Lockheed  Martin  Missiles  and  Space 


I 

PACEOS  has  been  under  development  over  the  past  ten  years  by  heritage  Lockheed 
Martin  personnel  primarily  at  Valley  Forge  PA  under  IR&D  funding  PACEOS  is  an 
engineering  tool  used  primarily  in  the  design  and  development  of  satellite  remote  sensing 
systems.  PACEOS  consists  of  a  series  of  modules  that  are  run  in  sequence.  The  modules 
simulate  the  data  flow  from  its  source  at  the  earth’s  surface  and/or  atmosphere  through  the 
sensor  and  focal  plane,  the  platform  and  onboard  processing,  the  downlink,  ground 
processing,  and  extraction  of  information  by  the  end  user.  Once  the  procedures  are  set  up 
for  a  given  system,  the  effects  of  changes  in  any  of  the  scene,  sensor,  platform,  or  data 
handling  parameters  can  easily  be  assessed  by  comparison  of  the  values  extracted  from 
algorithms,  such  as  EDR  algorithms,  with  the  known  input  scene  conditions. 

This  paper  discusses  the  use  of  MODTRAN  in  the  scene  generation  portion 
of  PACEOS.  It  includes  sections  on: 

-  PACEOS  generation  of  MODTRAN  input  files 

-  cloud  layer  radiance  and  transmission  models  within  a  scene 

-  variations  in  scene  conditions  using  non-homogeneous  atmospheric 
profiles  across  the  scenes 

-  generation  of  scenes  from  material  data,  varying  albedo/emissivities 
or  varying  temperatures 

-  options  to  select  different  environmental  conditions  while 
running  PACEOS 

-  output  consistency  comparisons  on  the  Sun  and  DEC  Alpha  workstations 

-  generation  of  the  scene  at  the  sensor  aperture 

Sample  results  for  a  variety  of  conditions  will  be  shown. 

We  are  currently  upgrading  PACEOS  to  use  MODTRAN3.5  version  1.2  from 
MODTRAN3. 
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BOMB  DAMAGE  ASSESSMENT  THROUGH  SPECTRAL/RADIOMETRIC 
ANALYSIS  OF  PLUME  C02  EMISSIONS 


William  E.  Wilcox  &  Larrene  K.  Harada 


W.  J.  Schafer  &  Associates 
321  Billerica  Road 
Chelmsford,  MA  08124-4191 


We  have  developed  a  technique  to  assess  underground  facility  bomb  damage  by 
measuring  C02  emissions  from  the  plume  of  the  bomb's  explosion.  Spectral 
measurements  of  various  explosions  have  been  made  using  a  CIGARS  spectrometer  with 
a  high  rep  rate  (60Hz)  and  moderate  resolution  (3cm- 1)  spectrometer  in  the  3-5  micron 
band. 

The  resulting  spectra  show  a  strong  C02  emission  feature  which  correlates  to  the 
combustion  of  bomb  fuel  products  and  provides  information  on  the  type  of  explosion 
being  observed. 

The  analysis  requires  MODTRAN  simulations  of  atmospheric  transmission  and 
background  radiance,  as  well  as  HTTRAN/HTTEMP  calculations  of  C02  emission  at  high 
temperature. 
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ANALYSIS  OF  VISIBLE  TRANSMISSOMETER 
DATA  IN  THE  COASTAL  SURF  ZONE 


Michael  Zugger  and  C.R.  Philbrick 

The  Pennsylvania  State  University 
Electrical  Engineering  Department 
University  Park  PA  16802 


The  optical  propagation  conditions  in  coastal  environments  are  highly  variable  due  to 
the  complex  interaction  of  the  maritime  conditions  with  either  urban  or  rural 
atmospheres.  The  breaking  surf  is  a  generator  of  aerosol  particles,  which  brings  further 
uncertainty  to  the  description  of  the  optical  propagation  path  over  the  surf  zone.  This 
propagation  path  may  affect  the  performance  of  several  types  of  sensors  and  operational 
systems. 

In  this  presentation,  the  data  obtained  using  an  optical  transmissometer  during  a  recent 
campaign  at  Scripps  Oceanographic  Institution  will  be  described.  The  instmment  was 
operated  along  a  diagonal  path  across  the  surf  zone  between  the  shore  and  locations  on 
the  Scripps  pier.  Measurements  taken  using  the  two  instruments  allow  a  comparison  and 
correlation  of  the  transmittances.  These  transmittances  will  be  correlated  with 
temperature  profiles  from  a  vertical  array  of  sensors  in  the  surf  zone,  as  well  as 
information  on  wind  and  surf  conditions.  Transmission  data  will  also  be  compared  with 
predictions  of  atmospheric  models.  This  first  report  on  the  results  will  also  discuss  the 
design  of  the  visible  transmissometer.  The  data  was  taken  as  part  of  the  Navy’s  Electro- 
Optical  Propagation  Assessment  in  the  Coastal  Environment  (EOPACE)  program  during 
the  period  24  March  to  1 1  April  1997. 
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TEMPORAL  AND  SPATIAL  CHARACTERISTICS  OF  SURF-ZONE  PLUMES 


Christophe  F.  Bas  and  C.  Russell  Philbrick 


The  Pennsylvania  State  University 
Electrical  Engineering  Department 
University  Park  PA  16802 


One  of  the  goals  of  the  investigation,  Electro-Optical  Propagation  Assessment  in  the 
Coastal  Environment  (EOPACE),  is  to  enhance  our  understanding  of  the  impact  that  surf 
generated  aerosols  have  on  the  propagation  of  optical  radiation  through  coastal  regions. 
A  bi-static  lidar  system  uses  a  fan  beam  to  illuminate  the  atmospheric  region  above  the 
surf  which  permits  measurements  of  the  vertical  and  horizontal  extent  of  aerosol  plume 
structures.  The  scattered  radiation  is  imaged  by  two  sensitive  CCD  cameras.  Three 
measurement  phases  of  the  experiment  have  provided  data  on  the  surf  generated  coastal 
plumes  at  two  locations  and  under  a  range  of  different  meteorological  conditions.  The 
spatial  and  temporal  characteristics  of  these  "coastal  plumes  are  described. 
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USE  OF  ATMOSPHERIC  AND  SIMULATION  MODELS 
FOR  IMAGING  SENSOR  EVALUATIONS 


S.J.  Westmoreland,  F.C.  Mertz 

Photon  Research  Associates  Inc. 
5720  Oberlin  Dr. 

San  Diego,  CA  92121 

i 

W.M.  Cornette 

National  Imagery  and  Mapping  Agency 
4600  Sangamore  Rd.,  MS  D-33 
Bethesda,  MD  20816 


As  part  of  evaluating  remotely  sensed  data,  atmospheric  effects  must  be  removed  or 
otherwise  accounted  for.  In  this  paper  we  document  a  process  for  integrating 
atmospheric  model  results  with  sensor  spectral  response  to  (1)  remove  atmospheric 
effects  from  image  data  or  (2)  add  atmospheric  effects  to  surface  measured  spectra  for 
comparison  with  at-sensor  measurements.  The  environment  developed  for  this  process 
permits  interactive  visualization  and  assessment  of  atmospheric  effects  on  both  measured 
and  modeled  data,  Examples  from  sensor  evaluation  and  calibration  analyses  are 
presented  including  comparisons  of  measured  and  modeled  radiance  and  tranmission  data 
generated  by  MOSART  and  MODTRAN,  effectiveness  of  atmospheric  removal  for 
sensor  calibration,  and  sensitivity  analyses  of  spectral  response  to  variation  in 
atmospheric  parameters. 
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INFRARED  SPECTRAL  DATA  ANALYSIS  AND  REMOTE  SENSING 
USING  PATTERN  RECOGNITION  ALGORITHMS 


S.J.  Lipson,  R.B.  Lockwood,  D.L.  Vititoe,  C.L.  Allred,  and  W.A.M.  Blumberg 
Phillips  Laboratory,  Geophysics  Directorate,  Hanscom  AFB,  MA 

P.S.  Armstrong 

Stewart  Radiance  Laboratory,  Bedford,  MA 

M.P.  Jacobson,  S.L.  Coy,  and  R.W.  Field 
Massachusetts  Institute  of  Technology,  Cambridge,  MA 


In  remote  sensing  of  atmospheric  spectra,  the  complexity  of  analyzing  data  may  be 
compounded  by  extreme  departures  from  thermodynamic  equilibrium  and  by  spectral 
overlaps,  temporal  variation,  and  lack  of  knowledge  of  line-of-sight  effects.  We  have 
developed  algorithms  for  the  analysis  of  infrared  spectra  taken  under  these  conditions, 
and  applied  them  to  the  study  of  highly  vibrationally  and  rotationally  excited  carbon 
monoxide  observed  in  atmospheric  simulation  experiments.  The  data  were  taken  using 
the  cryogenic-background  LABCEDE  facility  at  the  Phillips  Laboratory.  Mixtures  of 
CO,  Ar,  and  N2  were  irradiated  with  a  pulsed  electron  beam,  and  the  resulting  time- 
resolved  infrared  CO  spectra  were  obtained  with  a  Michelson  interferometer  capable  of  2 
cm-1  resolution.  The  spectra  were  highly  self-absorbed  in  the  v=l-0  band,  but  not  in  the 
highly  vibrationally  and  rotationally  excited  bands.  The  complex,  overlapped,  self- 
absorbed  spectra  were  analyzed  using  pattern  recognition  algorithms,  and  excellent 
agreement  between  data  and  fit  was  obtained.  Vibrational  basis  sets  (patterns)  were 
obtained  both  from  synthetic  spectral  models  and  from  the  XCC  (Extended  Cross- 
Correlation)  method.  A  global  fit  using  this  combined  model  was  performed,  in  which 
the  time  dependence  of  the  known  basis  sets  and  the  XCC  patterns  was  determined. 
Applications  of  this  method  include  remote  sensing  of  environments  in  which  the  optical 
properties  of  the  line  of  sight  are  poorly  known,  where  spectroscopic  basis  sets  are 
unavailable,  or  where  spectra  change  with  time  or  viewing  angle,  such  as  in  the  aurora  or 
other  highly  structured  atmospheric  scenes. 

This  work  was  supported  by  the  Air  Force  Office  of  Scientific  Research. 
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COMPARISON  OF  THE  INFRARED  RADIATION  MEASUREMENTS 
WITH  CALCULATIONS  USING  LOWTRAN  CODE 


Beiying  Wu*  * 

Institute  of  Atmospheric  Physics 
Beijing  100029 
China 


A  field  experiment  was  carried  out  in  Xinxiang,  China  in  1993  to  validate  the 
LOWTRAN  code.  An  pyrgeometer  measured  The  infrared  downward  and  upward 
radiative  fluxes  at  the  surface  level  for  the  wavelength  greater  than  4  micron.  An  infrared 
thermometer  measured  the  brightness  temperature  of  the  sky.  Atmospheric  temperature, 
relative  humidity,  wind  speed  and  direction,  visibility  and  cloud  type  and  amount  at  the 
surface  level  as  well  as  the  profiles  measurements  by  a  tether  balloon  system  were 
recorded.  Radiosonde  was  also  released  in  the  morning,  afternoon  and  evening  at  the 
Zhengzhou  station  about  100km  south  of  Xinxiang.  Observations  were  conducted  every 
hour  for  one  week.  The  LOWTRAN  code  was  used  to  calculate  the  downward  md 
upward  radiative  fluxes  and  the  brightness  temperamre  of  the  sky  at  various  observation 
direction  using  the  observed  atmospheric  properties.  This  paper  compares  the  calculation 
with  observations  and  analyzes  the  cause  of  the  difference  between  them. 

*  Visitor  of  Atmospheric  Chemistry  Division,  National  Center  for  Atmospheric  Research, 
Boulder  30807. 
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MICROWAVE  TRANSMISSION  IN  LIQUID  WATER  CLOUDS: 
UNCERTAINTIES  AT  LOW  TEMPERATURES 


Alan  E.  Lipton  and  Michael  K.  GrifTm 

Optical  Effects  Division 
Phillips  Laboratory 

Hanscom  Air  Force  Base,  Massachusetts,  01731-3010 

Alan  G.  Ling 

Radex,  Inc. 

3  Preston  Ct 

Bedford,  Massachusetts,  01730 


For  many  years,  the  standard  for  modeling  the  attenuation  of  microwave  radiation  by 
cloud  liquid  water  has  been  to  rely  on  the  Ray  model  for  the  dielectric  properties  of 
water,  while  employing  the  Rayleigh  absorption  approximation.  The  Ray  model  is  used 
in  the  FASE  and  RADTRAN  radiative  transfer  codes.  Recently,  Liebe,  Manabe,  and 
Hufford  have  developed  models  for  the  dielectric  properties  of  water  using  two  relaxation 
terms,  in  contrast  to  the  single  term  used  by  Ray.  The  models  of  Ray  and  of  Liebe, 
Manabe,  and  Hufford  yield  very  similar  values  of  attenuation  for  temperatures  above  10 
degrees  Celsius,  within  the  frequency  range  of  0  to  250  GHz.  At  lower  temperatures, 
however,  there  are  substantial  spectral  differences  between  attenuation  values  produced 
by  the  different  models.  Considering  that  liquid  water  can  occur  at  temperatures  as  low 
as  -40  degrees  Celsius,  the  model  differences  indicate  significant  uncertainties  in  cloud 
remote  sensing  results  derived  with  the  aid  of  the  models. 
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ECRIN'96  :  AN  AIRBORNE  CAMPAIGN  FOR  THE  RADIANCE 
MEASUREMENT  OF  STRATOCUMULUS  CLOUDS 


P.  Siinoneau‘,  C.  Malherbe’,  A.  Boischot’,  G.  Durand',  B.  Rosier’,  A.  Delannoy' 
P.  H.  FIamant^  J.  Pelon^  L.  Sauvage*,  H.  Dirand* 

'  Office  National  d'Etudes  et  et  Recherche  Aerospatiale 
Palaiseau,  France 

*  Laboratoire  de  Meteorologie  Dynamique 
Ecole  Polytechnique,  Palaiseau,  France 

*  Service  d’Aeronomie  du  CNRS 
Universite  Pierre  et  Marie  Curie,  Paris,  France 


ECRIN'96  (etude  du  Champs  Radiatif  Infrarouge  des  Nuages)  has  been  achieved  for 
the  correlative  measurement  of  the  infrared  radiation  coming  from  low  altitude  clouds,  in 
situ  measurements,  and  lidar  measurements  on  clouds.  The  experiment  took  place  in 
Britanny  near  Brest  (West  coast  of  FRANCE)  in  October  96  with  the  participation  of  3 
aircraft  :  a  Caravelle  operated  by  the  Centre  d'Essai  en  Vol,  a  turboprop  Fokker-27  so 
called  ARAT  (Avion  de  Recherche  Atmospherique  et  Teledetection)  and  a  turboprop 
MERLIN  4  operated  by  Meteo-France. 

The  aircraft  were  instrumented  with  the  following  sensors  : 

-  an  airborne  infrared  cryogenic  spectrometer  and  an  infrared  imager  loaded  on  the 
Caravelle  for  radiance  measurements 

-  an  airborne  Lidar  for  the  sounding  of  high  altitude  clouds  and  microphysic 
measurements,  loaded  on  the  Fokker-27 

-  microphysic  and  thermodynamic  sensors  for  in-situ  measurements,  mounted  on 
the  Merlin  4 

-  fluxes  sensors  mounted  on  the  Fokker-27  and  the  ARAT 

The  measurement  strategy  and  preliminaries  results  (radiance  spectra,  lidar 
measurements  and  in-situ  microphysic  measurements)  are  presented. 
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ATMOSPHERIC  CORRECTION  THROUGH  TRANSMISSIVE  CLOUDS 


J.  William  Snow,  Hsiao-hua  K.  Burke,  Michael  P.  Jordan, 
Daniel  C.  Peduzzi  and  Kristine  E.  Rhoades 


MIT  Lincoln  Laboratory 
244  Wood  Street 
Lexington,  MA  02173 


From  space,  the  extraction  of  surface  parameters  by  using  IR  sensors  has  long  been 
impeded  by  clouds.  When  clouds  are  present,  retrieval  is  not  performed  and  thus  surface 
information  is  unavailable.  The  high  frequency  of  occurrence  of  clouds  calls  for 
technology  development  to  minimize  this  shortcoming.  In  particular,  transmissive  clouds 
offer  the  greatest  potential  since  statistically  their  global  coverage  extends  to  more  than 
20%,  and  typically  they  are  less  than  0.5  in  optical  depth,  which  makes  the  correction  of 
surface  radiances  more  tractable. 

A  methodology  for  the  determination  of  the  atmospheric  conditions  indicative  of 
cloud  presence  and  derivation  of  its  physical  thickness  was  developed  based  on  water 
vapor  and  temperature  profiles.  Transmissive  clouds  are  identified  and  parameterized, 
enabling  correct  surface  radiance  to  be  retrieved.  A  surface  radiance  retrieval  algorithm 
based  on  GOES  MWIR  satellite  data  which  utilizes  readily  available  atmospheric 
information  and  MODTRAN  is  used.  The  results  of  surface  radiance  retrieved  through 
these  transmissive  cloud  cases  are  encouraging  -  the  derived  surface  temperature  bias  and 
rms  were  shown  to  be  0. 1  K  and  1.0  K,  respectively. 

This  work  was  supported  under  the  auspices  of  the  Lincoln  Laboratory  Innovative 
Research  Program.  The  Innovative  Research  Program  is  supported  principally  by  the 
Department  of  the  Air  Force  under  Contract  F19628-95-C-0002. 
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ENVIRONMENTAL  TRANSMISSION  EFFECTS  WITHIN 
A  DISTRIBUTED  SIMULATION 


Joel  B.  Mozer,  Thomas  R.  Caudill,  Capt.  Mark  M.  Edwards 

USAF  Phillips  Laboratory 
Geophysics  Directorate 
Optical  Effects  Division 
Hanscom  AFB,  MA  01731  , 

Steven  A.  Ayer  and  Guy  P.  Seeley 

Radex,  Inc. 

3  Preston  Ct 

Bedford,  Massachusetts,  01730 


Distributed  simulation  is  an  increasingly  important  part  of  the  modeling  and 
simulation  activities  in  the  DoD  due  to  diminishing  development  and  training  budgets 
relative  to  the  sophistication  of  modem  weapon  systems.  A  crucial,  but  often  ignored, 
component  of  this  type  of  simulation  is  the  effect  of  the  natural  environment  on  an 
engagement  scenario.  For  example,  the  impact  of  atmospheric  constituents  such  as 
clouds,  wind,  fog,  and  haze  are  not  typically  included  in  a  weapon  system  simulation. 
However,  the  existence  of  these  constituents  in  a  real  system  can  mean  the  difference 
between  a  weapon  being  successful,  or  completely  useless  in  a  given  situation. 

We  have  developed  a  prototype  capability  to  generate  and  disseminate  atmospheric 
data  within  a  Distributed  Interactive  Simulation  (DIS)  environment.  This  capability 
utilizes  a  variety  of  models  and  simulations  to  generate  the  appropriate  3D  atmospheric 
constituent  data  using  a  combination  of  satellite  analysis  algorithms,  mesoscale  numerical 
weather  prediction  model  diagnostics,  and  Phillips  Laboratory's  Cloud  Scene  Sunulation 
Model  (CSSM).  The  transmission  effects  of  the  simulated  environment  on  a  weapon 
system  are  determined  via  a  modified  version  of  the  Electro-Optical  Tactical  Decision 
Aid  (EO/TDA)  program.  This  paper  details  the  construction  of  this  prototype  capability 
and  outlines  the  needs  for  this  type  of  activity  to  be  effective  in  future  simulations. 
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GLOBAL  WARMING  STABILIZATION  MODELS: 
VEHICLE  EMISSIONS,  REnGERANTS,  AND  DEFORESTATION 


Steven  Dupree,  Paul  Yeaman,  Robert  Alwine,  and  Jefferey  Campbell 

Worcester  Polytechnic  Institute 
Worcester,  MA 


Evidence  clearly  indicates  that  global  warming  is  an  observable  phenomenon  likely  to 
be  linked  to  human  activity.  If  this  phenomenon  continues  at  its  current  rate,  the  effects 
could  be  disastrous  for  humankind  and  all  life  on  Earth.  As  one  of  the  major  superpowers 
of  the  world,  and  one  of  the  most  technologically  advanced,  it  is  the  responsibility  of  the 
United  States  to  lead  the  efforts  to  reduce  global  warming.  This  project  considers  some 
possible  reduction  methods,  including  various  policy  changes  and  legislative  actions, 
evaluating  each  in  terms  of  effectiveness  as  well  as  cost.  A  final  comparison  leads  to  a 
recommendation  of  a  set  of  policies  which  will  curb  the  effects  of  global  warming  in  as 
cost  effective  a  manner  as  possible. 
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MODTRAN4:  Simulating  Atmospheric  Radiation 
Budgets  in  the  Middle  Atmosphere: 


G.P.  Anderson,  J.H.  Chetwynd,  A.  Berk, 
L.S.  Bernstein,  H.E.  Snell,  A.  Ratkowski, 

Phillips  Laboratory 
Geophysics  Directorate 
Hanscom  AFB,  MA 

E.P.  Shettle 

Naval  Research  Laboratory 
Washington,  DC 


The  new  MODTRAN4  band  model,  with  its  correlated-k  Beer's  Law  algorithm,  can 
efficiently  and  correctly  calculate  the  scattering  and  absorption  signatures  of  realistic 
molecular,  aerosol  and  cloudy  environments  in  the  lower  and  middle  atmosphere.  The 
current  approach  for  molecular  scattering  accommodates  line  overlap  and  partial 
correlations  between  both  molecular  species  and  the  solar  irradiance,  while  maintaining 
band  model  spectral  resolution  at  2  to  15  cm'".  This  new  level  of  evolution  and  validation 
will  permit  improved  syntheses  and  analyses  of  actinic  (direct  plus  scattered)  solar  and 
thermal  energy  sources  and  their  flux  divergence  within  the  stratosphere  and  troposphere. 

Rapid  calculation  of  IR  cooling  rates  for  a  variety  of  thermal  and  constituent  profiles 
demonstrates  the  flexibility  of  the  approach  and  the  sensitivity  to  latitudinal  and  dynamic 
variability.  Validation  is  provided  through  two  avenues.  The  first  involves  direct 
comparisons  with  line-by-line  calculations,  as  exemplified  by  FASE  (the  line-by-line 
algorithm  joindy  developed  from  FASCODE  by  DoD  and  DOE)  which  provides  the 
molecular  standard  for  layer  effective  optical  depths,  single  scattering  albedos,  and 
transmittances.  This  enables  the  more  primitive  MODTRAN4  algorithm  to  be  refined 
for  both  more  flexible  spectral  resolution  and  efficient/accurate  determination  of  these 
essential  layer  quantities  necessary  for  multiple  scattering  applications:  e.g.  DISORT. 
The  second  validation  step  centers  on  comparisons  against  a  variety  of  measurements, 
mostly  airborne  visible  and  IR  up-welling  radiances,  including  both  clear  and  clouded 
skies.  A  particular  attempt  will  be  made  to  examine  the  impact  of  polar  stratospheric 
clouds. 
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MWIR  Earthlimb  Scene 
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SHARC  and  SAMM  Validation 
Using  FASCODE  and  MODTRAN 
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Provide  statistical  structure  parameters  for  SWIR,  MWIR,  MLWIR,  and 
LWIR  background  scenes  containing  stochastic  structure 


SIG1.0 

SHARC  Image  Generator 
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Simulate  atmospheric  background  clutter  with  sufficient 
accuracy  to  meet  SBIRS  design  requirements 


Atmospheric  Background 
Model  Assessment 
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Temperature  Statistics 
Data  Base 


Atmospheric  Background 
Model  Validation  Using 
MSX  Data 


BTH  Broad  MWIR  Band  Data 
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Simulate  all-altitude  atmospheric  background  clutter  with  sufficient 
accuracy  to  meet  SBIRS  design  requirements 
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background  clutter  with  sufficient 
accuracy  to  meet  TMD  and  NMD 
design  requirements 


SBIRS  Model  Toolkit  Assessment 
Using  MSTI-3  and  ARES  Data 


SBIRS  Model  Toolkit  Uncertainty  Assessment 


SBIRS  Model  Toolkit  Assessment 
Using  MSTI-3  Data 


rxii:»«t: -  occe:  -.aKCtT 
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Measured  and  Modeled  SWIR  Scenes  Radiance  Histogram  Comparlsoir^ 
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Single  point  access  to  PL  IR  atmospheric  and  celestial  background  codes 
Embedded  expert  system  to  facilitate  code  applications  by  non-expert  users 


PLEXUS  3.0/3.0  Nl 
Phillips  Lab  Expert  Unified  Software 
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SHARC  4.0/SIG  1.0 
MWIR  Earthlimb  Scene 


New  Directions  for  Opticai 
Background  Codes 
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•  Important  role  in  defining  environmental  requirements  for 
systems  and  battiespace  simuiations 
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The  status  of  MOSART,  along  with  these  upgrades  will  be  presented.  Version  1.50  is 
planned  for  release  in  mid- August  1997,  with  a  major  upgrade  for  Version  2.00  planned 
for  March  1998.  The  plans  for  both  of  these  releases  will  be  presented,  together  with 
information  on  how  to  obtain  the  code  and  who  to  contact  for  technical  information  and 
"bug"  reporting. 


Introduction 


The  Moderate  Spectral  Atmospheric  Radiance  and  Transmittance  (MOSART)  code 
is  a  U.S.  Government  standard  code  with  the  ability  to  support  scene  and  signature 
simulations.  It  grew  out  of  a  merging  of  the  MODTR.AN  and  AP.ART  codes  wdth  some 
additional  capabilities  added.  Recommended  platforms  are  workstations  (SGI,  Sun,  HP, 
DEC  Alpha),  although  it  wdll  run  on  a  high  end  PC  (e.g.,  200  MHz.  64  Mb  RAM,  50-100 
Mb  of  free  disk  space).  It  has  been  thoroughly  tested  on  machines  that  use  the  IEEE  754 
floating  point  standard  (e.g.,  SGI,  Sun,  HP),  but  there  may  be  some  minor  numerical 
bugs  encountered  on  VAX/VMS  machines  due  to  the  hard  underflow  (i.e.,  one  knowm  case 
cr^hed  on  a  VAX  but  ran  on  the  other  IEEE  754-compliant  machines;  the  fix  was  simple 
and  was  purely  numerical). 

Except  for  DISORT,  MOSART  Version  1.50  contains  all  MODTRAN  3.5 
capabilities  and  has  been  compared  against  it  for  a  number  of  conditions  by  Spectral 
Sciences,  Inc.  and  Photon  Research  Associates.  Inc.  The  version  of  the  DISORT  code  in 
MODTRAN  3.5  is  still  under  beta  testing  and  further  work  is  on-going  at  Spectral 
Sciences,  Inc.  (SSI)  and  at  the  Air  Force  Research  Laboratory,  Hanscom  AFB,  M.A.  As 
soon  as  it  is  approved  by  Air  Force  Research  Laboratorv,  it  will  be  incorporated  into 
MOSART. 

As  with  MODTRAN,  MOSART  uses  a  1  cm"'  resolution  set  of  band  parameters 
derived  from  the  1997  HITRAN  line  atlas,  including  the  latest  change  to  partition  the  I  cm"’ 
region  into  four  sections  for  better  representation  of  the  continuum.  Based  on  comparisons 
with  degraded  spectral  created  by  the  line-by-line  code  FASCODE.  the  actual  resolution  of 
both  MODTRAN  and  MOSART  are  slightly  greater  than  1  cm'-  and  are  usually  advertised  as 
2  em'h  In  the  UV/  Visible  region,  a  few'  more  molecular  absorption  bands  are  included  at  a 
coarser  resolution,  and  the  multiple  scatter  model  is  different. 

The  viewpath  geometries  include  all  MODTR.AN3  geometries,  plus  some  new 
geometries,  including  Observ'cr-Source.  Obserx'er-Background.  Ob,server-Source- 
Background,  Earthlimb,  Horizontal,  At-Source  (i.e.,  Obsen'er  &  Source  are  co-located). 
The  user  can  specify  a  number  of  different  parameters  (e.g.,  observ'cr  altitude,  source 
altitude,  and  range;'obser\'er  altitudeAatitude/longitude  and  tangent  point  latitudc/longitude 
—  or  tangent  altitude,  for  a  couple  of  exam.ples),  any  potential  geometry  is  capable  of  being 
requested.  An  X- Windows/Motif  Input  GUI  is  available,  as  is  a  simple  scrolling  input  file 
builder.  However,  once  the  user  becomes  familiar  w'ith  the  input  file  stmeture  (it  is  self- 
documenting  with  full  text  descriptions  of  the  inputs  being  included  in  the  input  file),  an 
editor  like  vi  or  emacs  is  quite  adequate. 


Moderate  Spectral  Atmospheric 
Radiance  and  Transmittance 
(MOSART)  Computer  Code 
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MOSART  has  been  installed  and  tested  on  a  number  of  machine,  including: 

•  SGI  (IRIX) 

•  Sun  (OS  &  Solaris) 

•  HP  9000 

•  Data  General 

•  IBM  3094  (VM/CMS)  (Early  Version) 

•  High  End  PCs  (Lahey  &  Microsoft  compilers) 

•  DEC  VAXA-hMS  &  Alpha  3000 

•  Prime  (Early  Version) 

Although  installed  and  tested,  a  numerical  representation  problem  exists  on  DEC 
machines.  The  problem  is  related  to  DEC'S  non-compliance  with  IEEE  Standard  754,  and 
is  easily  fixable,  once  time  is  available  on  this  type  of  machine. 

To  install  MOSART  (values  are  for  a  SGI  machine),  50-100  Mbytes  of  Disk 
Space  is  required,  along  with  a  FORTRAN  77  Compiler  (ANSI  X3.9-1 978);  libraries 
satisfying  MIL-STD  1753  is  helpful,  but  not  necessaiy.  A  requirement  for  X- 
Window's/hfOTIF  is  optional  (for  GUI  only). 

To  date,  the  MOSART  code  has  been  utilized  by  a  number  of  simulations  and 
applications,  including: 

•  Synthetic  Scene  Generation  Model  (SSGM)  EMDO-NRE/PRA 

9 

•  Environment  Effect  in  Distributed  Interactive  Systems  (E“dis) 
DMSO./NRL,/PR.A/SSI 

•  Flexible  Infrared  Signature  Tools  (FIST)  RL.'N.41C/PR.4 

•  SPIRITS  (in  development)  PL/SSI 

•  SAMM2  (in  planning)  PUSSI 

•  Scene  Simulators: 

•  GENESSIS  (PRA) 

•  CLDSIM  (PR4) 

•  GCI  Toolkit  (PR.4) 

•  SensorVision  (PSI) 

•  CloudScape  (Visidyne) 

One  question  that  a  number  of  users  raise  is  why  two  codes  (i.e.,  M0DTR.4N 
and  MOSART)  and  which  one  should  be  used.  M0DTR.4N  is  a  small,  efficient  code 
with  representative  environmiental  conditions  (e.g.,  six  model  atmospheres)  that  performs 
point-to-point  or  point-to-space  calculations.  MODTRAN  uses  llne-of-sight  radiative 
transfer  (i.e.,  a  “soda  straw”  viewpoint).  MOSART.  on  the  other  hand,  has  access  to 

•  global  environmental  conditions  and  climatologies 
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•  global  terrain  altitude,  scene,  and  materia!  data  bases 

and  provides  output  to  support  scene  simulation  and  signature  models.  MOSART 
calculates  a  comprehensive  radiative  environment,  but  is  rather  large  and  is  not  very  fast. 
Each  code  has  its  own  purpose. 


Code  Status 

The  current  (i.e.,  as  of  June  1997)  version  of  the  MOSART  code  is  Version  1 .41 , 
which  was  delivered  in  September  1995.  This  version  of  MOSART  has  been  approved 
for  public  release,  distribution  unlimited  by  Air  Force  Public  Affairs  and  Ballistic  Missile 
Defense  Office;  release  by  Phillips  Laboratory  is  pending  code  modifications,  “bug”  fixes, 
and  upgrades.  The  MOSART  documentation  has  been  released  as  PL-TR-94-2244  (4 
Vols)  in  May  1996.  The  documentation  consists  of: 


1 .  Installation  Reference  Manual 

2.  Users’  Reference  Manual 

3.  Technical  Reference  Manual 

4.  Software  Reference  Manual 


As  of  June  1997,  Version  1.41  has  undergone  significant  testing  and  evaluation  by  several 
users,  resulting  in  some  general  clean-up,  “bug”  fixes  and  code  modifications, 
documentation  corrections,  clarifications  (resulting  fi-om  questions  from  various  users), 
and  an  explanation  of  FTNCHEK  w'amings.  The  major  upgrades  that  have  been 
implemented  since  the  initial  September  1995  release  are: 


1 .  calculation  of  ephemeris  time  improved 

2.  User-Defined  atmospheres  upgraded 

3.  snow7ice  in  temperate  climates  removed 

4.  consistency  in  fractions  vs.  % 

5.  oxygen  continuum  factor  (MODTRAN  3  compatibility) 

6.  molecular  concentration  normalization  (User-Defined ) 

7.  additional  checking  of  user  input 

8.  'To  Space"  output 

9.  numerical  "closeness"  on  machines  with  64-bit  register  arithmetic 


Code  Description 
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A  personal  message: 

While  we  all  miss  Frank  Kneizys,  we  all  miss  him  differently. 

His  calm  manner  and  great  intellect  propelled  the  codes  for 
over  a  decade.  Without  Frank  there  would  be  no  MODTRAN.  He  realized 
how  important  it  is  to  have  a  moderate  resolution  quick  code. 

Frank  also  co-claimed  responsibility  for  the  birth  of  FASCODE.  He 
even  had  back-of-the-napkin  sketches  of  the  original  deconvolution.  The 
collaboration  between  Tony  Clough  and  Frank  on  the  invention  of  FASCODE 
appeared  to  be  a  constant  struggle  towards  excellence  in  science,  and  sometimes 
even  a  constant  struggle.  Tony  and  Frank  worked  so  well  together. 

The  Geophysics  Lab  has  products  to  be  proud  of  because  of  their  work. 
Each  year  hundreds  of  new  people  discover  that  most  of  their  problems 
have  already  been  solved.  About  500  copies  of  MODTRAN  are  picked  up  from 
our  web  site  each  year. 

Frank’s  death  came  as  a  great  shock  to  us.  We  saw  him  not  long  before 
his  death  and  he  seemed  fine.  I  miss  him  most  as  a  fiiend.  Someone  to  walk 
with  during  lunch,  trading  jokes  and  stories  of  family.  He  was  one  of  my 
mentors  and  my  fiiend.  He  can  never  be  replaced.  I  will  always  remember  him. 

Jim  Chetwynd 
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Reduction  of  the  Fourth-Order  Asymmetric-Rotor  Hamiltonian 

F.  X.  Kneizys,  J.  N.  Freeduan,*  and  S.  A.  Clodcii 
An  I'oru  Cambridge  Research  Uboralories,  Bedford,  Massachusetts 
(Received  I  November  1905) 

for®moIeLtfo7orffrL“ ~  t  «ymmetric-rotor  Hamiltonian 

dents  has  been  rXed  to  19  TeTScTenu 

10  coefficients  of  sixth  power  in  th  ^  Power.  and 


TN  a  recent  ’  -  definition  of 

1  udno  papers,*  »  Chung  and  Parker,  second-order  distortion  constants  of  Kivelson  and 

^ilh  b^  formalism  of  Nielsen.  Amat.  and  Gold-  Wilson'*;  and  are  the  fourth-order  corrections  to 
r-  ^  y  derived  the  fourth-order  con-  For  the  case  of  orthorhombic  symmetry  there 

toniarFrn  b "®i^^‘"C-rotor  Hamil-  are  nine  independent  15  nonvaniS^ng  p i,*  and 
tomaii.  From  their  results,  it  is  now  possible  for  the  f05  MySe,].*  ^ 

iilf  theoretical  significance  to  coefficients  f^y  the  use  of  the  commutation  relations  for  angular 

up  to  the  sixth  power  in  the  angular  momentum.  The  momentum,  ^ 

results  of  Chung  and  Parker  were  obtained  by  apply-  PaP/s- PfP„=-ihP^  (O) 

mg  symmetry  arguments  to  the  coefficients  of  th^  r  i,  a  ■ 

asymmetric-rotor  Hamiltonian  In  this  nanp  b  (where  a,  y  are  in  cyclic  order),  and  with  the  appro- 
that  by  me  of  theT„X!;;omen,LT,iZu^^^  ““if'tfhle 

relations,  the  asymmelric-rolor  Hamiltonian  of  ortho  fono^'n  Hamiltonian.  The  trans- 

“KiS"  “  “3'"  '“™  -Ti  "“s  1c  t“fo«f  ;?er? 

x‘"ist:^„trprrcS  'Jf  5?:? 

..  HANn.X0H™__™H0RH0MBfC 


POINT  GROUPS 

The  Hamiltonian  of  an  asymmetric  rotor  of  ortho¬ 
rhombic  s^etry  in  the  IIP  representation*  is  of  the 
following  form; 

= A.  **-f2lPrr*-f  58P,*-f  (SP,* 

'f'*22(raPrJ+Pap^j)  PaP0PyPt 
afiyS 


+  MP«^P.*+P.*Py^)+^c(Py^Py*+  P*Py^) 

+  MPr^Py*+  Py*P.^)  +  i>i(Pf‘P/+  P/P.2) 

+  P.*Px*)  +  ‘I>10(Px'Pv*P.*+  P.^Py^Py^),  (3) 

+  Z)  Lo^y^vlPaPpPyPiP.Py,  (1) 

b  LtAov  ^=2H-J/i*(3r^2r/-2r,.+p,-p8-p, 

^'■^'^hose  of  Chung  ^  ^  -fL-fo, o 

tion  coefficients  of  Chung  and  Parker.^-^  The  term  h  ^  53+ 4ft  (Sn  2r9— 2r7*“P7+p8— pg 

.  -  Chi. 

27!  ^38  0957')“’’’  “•  “•  J-  Phys.  _  J 

:g-  and  H.  H.  Nielsen,  J.  Chem.  Phys.  27.  845  (1957)  ^p,.)+fi^'l>ii, 

T..i(..+2,.+n.+n,+|„.+j,.., 

II.  27  (I9«"^In  ^6-  T(^6+2rft+p6+p9+-|p^4+Xp^s) 

(.95?). -  J- 

the  (u"^)  A  67.  359  (196?).^"'"'  -  '"’• 

cussed  by  King,  Haincr,  and  Cross!  m.inncr  dis-  ”  P'cfcc,  N.  DiCiannI,  and  R.  H.  J.ackson,  J.  Chem.  Phys. 

oo»  730  (1963). 

2552 
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.  Phys.  20,  1575 
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Total  Measured  =  52.93  W 


1000  _  1200  1400 

Wavenumber  (cm 


Reports  of  Dispcrepancies  between  Measurement  and  Models  in  the  Shortwave 
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Absolute  Solar  Transmittance  Interferometer  (ASTI) 
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Difference:  ASTI  -  LBLRTM  April  18,  1996  2330  Z  Case 


iD 

iO 

o 

lO 

lO 

iO 

iO 

o 

lO 

vO 

O 

o 

o 

o 

o 

o 

o 

o 

X 

X 

X 

X 

X 

X 

X 

X 

0^ 

(^} 

c\? 

I  -  I 


((,.mo  JS  ^UlD)/f\Uj)  90U9J9jjia  aouBipB^ 


7000  7200  7400  7600  7800  8000  8200 

Wavenumber  (cm  ') 


ASTI  vs.  LBLRTM  April  18,  1996  2330  Z  Case 


CO 

o 

X 


in 

O  CO 

CO 

o 

o 

o 

X 

X 

X 

in 

q 

lO 

— 

js  ^uio)/^m)  aouBipe^ 

\Ci 

o 

X 

o 

lO 


o 


107 


8500  9000  9500 

Wavenumber  (cm"') 


iC  if5  O  lO  lO 


o 

o 

o 

o 

<•—4 

-.-H 

X 

X 

X 

X 

c\? 

I 


((,_uio  JS  ^mo)/MUi)  aaua-iajjiQ  aoueipieH 


108 


8500  9000  9500 

Wavenumber  (cm  ') 


ASTI  vs.  LBLRTM  April  18,  1996  2330  Z  Case 


CO 

CO 

lO 

o 

CO 

CO 

la 

o 

o 

o 

o 

o 

o 

X 

X 

X 

X 

X 

X 

lO 

o 

o 

iO 

q 

q 

id 

— i 

id 

((j_uio  JS  ^uio)/^ui)  aouBipB^ 


109 


7000  7200  7400  7600  7800  8000  8200 

Wavenumber  (cm'‘) 


Properties  of  the  continuum  and  discrete  absorption  bands  due  to  the  three  strongest  Oj  a' A  <—  Z”  transitions 
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20th  Annual  DoD  Atmospheric  Radiance  and  Transmittance  Meeting, 
Hanscom  AFB,  1 0-1 2  June,  1 997 


Accelerated  line-by-line  calculation 
of  absorption  coefficients  with  high 
numerical  accuracy 


Martin  Kuntz 

Institutfur  Meteorologie  und  Klimaforschung 
Forschungszentrum  Karlsruhe 


®  Introduction 

t  Efficient  calculation  of  the  Voigt  profile  function 
®  Determination  of  frequency  spacing 
®  Summation  of  spectral  lines 
t  Performance 
t-  Conclusion 
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introduction 


In  practice  calculating  absorption  coefficients  proves  to  be  the 
dominant  computational  expense  of  a  line-by-line  radiative 
transfer  model  (requires  up  to  90  %  of  execution  time). 


^lT,p,v)  =  AJT) 

"  - V - -  -v  -f  . 

■fee  strerigth  .  profile  tunctionD::;.' 


Three  questions  that  need  to  be  clarified 


®  How  to  calculate  the  profile  function  most  efficiently? 
(which  algorithm?) 

®  Where  to  calculate  the  profile  function  most  efficiently? 
(at  which  sampling  points?) 

®  How  to  perform  the  summation  most  efficiently? 
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Efficient  calculation  of  the  Voigt  profile  function 


Humlicek-Algorithm: 

Dividos  x,y-space  into  several  regions  \A/here  complex 
rational  polynomials  are  adequate  within  a  given  error 


Reduce  if-inquiries  applying  a  recursive  algorithm 
before;  ~2n  afterwards:  -2log.n 

Calculate  only  the  real  part  of  the  complex  Voigt  profile  function 
before:  ~  75  Flops  afterwards;  ~  19  Flops 

Automatically  switch  to  the  calculation  of  the  Lorentzian  function 
where  permissible 

Accelerated  calculation  ifj  remains  constant 
Performance: 

Factor:  --4.7  for  rel.  accuracy:  -2 1(F 
Factor;  ~  6  for  rel.  accuracy:  -1(F 


Determination  of  frequency  spacing 


idea:  Calculate  each  spectra!  line  as  accurate  as 

necessary  in  view  of  the  desired  absolute  accuracy 
and  interpolate  the  intermediate  values,  e.g.  by  a 
second  order  polynomial. 

problem:  Interpolation  error  50(y,jc,Ax)  depends  on  the  distance 
from  line  center  jc  and  the  frequency  spacing  Ax 


solution:  i'  Determine  the  highest  permissible  interpolation  error 
for  a  given  line  while  retaining  an  absolute  accuracy  5a 

®  Determine  Ax  in  such  a  way  that  the  permissible 
interpolation  error  60  is  not  violated! 
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Determination  of  frequency  spacing 
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i 


50  =  10''  50  =  10'- 
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Determination  of  frequency  spacing 


Voigt  cutoffs  for  different  parameters  y 


Approximation  by  Sparks 


120 


Forschungszentrum  Karlsruhe 
Technik  und  Umwelt 


:o 


o> 


CD 


£S  .  ^ 

H— ^  .CC.; 


o 


.S2  o 


.4^  .  CO 


121 


1100,0  ':1 100,2  1100,4  1100.6"  1100,8  1101,0 

'  A 


122 


Forschungszentrum  Karlsruhe 
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T-derivative  [cm  'K ']  Absorption  coefficient  [cm  '| 


Summation  of  spectral  lines 


Absorption  coefficient 


T-derivative  of  absorption  coefficient 
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Performance 


Microwindow  of  width  1cm'^ 

10  000  lines  +/-5  cm'^  around  the  microwindow 
of  which  ~  5000  lines  were  calculated 


Accuracy:  5a  =  10‘®cm‘'' 


— 

Av  =  510-^cm*‘’  • 

,  r  -Av  F  1  W  .cml 

height 

. .  new 

•  equidlist^f 

•new  • 

t2kmi 

liiljiigaMiiil 

143s- 

21 8.4s  • 

1rJ8S  . 

24km  i 

1485 

391 .7s 

1..59S 

36km  1 

158j0s  ^ 

1.54s 

745.3s 

2.23s 

48km 

•  148.7s  • 

1.43s 

:  632-7S  ’ 

2-15$ 

60km 

128.5s 

1.31s 

:  •432.7s 

1.71s 

72km 

-  . 

253.8s 

:  "1.31s 

Calculations  performed  with  different  accuracies 
(Av  =  1  1 0'^cm'^  height  =  72km) 


accuracy 

;  time  . 

memory 
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^  .3;46^  ^ 

i  98.0% 

.  lO^^cm” 

^  4.39s  • 

99.5% 
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Conclusion 


Advantages: 

t?  each  spectral  line  is  calculated  as  accurate  as  necessary  in  view 
of  the  desired  numerical  accuracy 

®-  numerical  accuracy  can  be  varied  continuously  during  running  time 
(e.g.  can  be  easily  increased  in  an  iterative  procedure) 

®  it  is  in  principle  possible  to  calculate  error  bars 

®  algorithm  may  be  readily  applied  to  alternative  line  shape  models 

t  increasing  the  spectral  resolution  does  not  necessarily  Imply  a 
higher  computational  workload  (only  some  additional  interpoations) 

®  improved  algorithm  for  the  calculation  of  the  Voigt  profile  function 

®  improved  frequency  spacing  for  the  representation  of  the 
Voigt  profile  function 

®  T-derivative  of  absorption  coefficient  is  calculated  simultaneously 
®  sampling  Intervals  can  be  varied  continuously 

nonequidistant  (quasi-equidistant)  frequency  spacing 

-  saving  of  valuable  computer  memory 

-  improves  computational  efficiency  of  the  radiative  transfer 

Disadvantages: 

®  rather  complicated  algorithm 
®  one  additional  parameter  to  specify 
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IMPROVED  ACCURACY  INFRARED  ABSORPTION  DETERMINATIONS 
IMPROVED  ACCU^i^  modified  KUBELKA  MUNK  MODEL 


W.  G.  Egan,  and  Q.  Alexis 

Natural  Sciences  Department 
York  College/CUNY 
Jamaica,  NY  11451 


In  order  to  further  improve  the  accuracy  of  infrared  aerosol  optical  properties 
determinations,  a  modification  to  the  original  Kubelka-Munk  scattering  model  was  made 
to  include  the  effects  of  sample  surface  reflection  on  the  determination  of  absorp 
i  e  Moafed  KuteSta  Mn^  Model:  M-K-M).  The  solufions  to  the  radiattye  tr^s  er 

theory  when  significant  scattering  occurs. 
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R,i  = 


^  I cos  g,  -  -  sin^ 

cos  t/j  +  iim^./nvy-  -  sin^T 


Erii/En 


i.ei* 


The.Kramers-Kronig  relation  links  In  r  to  cp, 


^  f”  In  '• 

H  Jn  -  V, 


\/ 


^  /  ^  /  \  /  \  ^  X  2.4 

^C./  /V'tC/^.C' 


phase  shift,  <5 


AMPLITUDE, 


Viewgraph  I:  a)  Vector  form  of  Fresnel's  reflection  equation 

b)  Kramers-Kronig  relation 

c)  Plot  of  Kramers-Kronig  relation  in  complex  plane 


•ANGLE  OF  INCIDENCE.  5. 


Viewgraph  II:  Idealized  polarized  reflectance  from  a  surface: 

Lower  curve  goes  to  zero  at  Brewster  angle; 

upper  curve  does  not  go  to  zero  because  of  scattering 
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Viewgraph  V: 


X  •  0 


X  •  d 


6-Sj 

<J'0 


FLUX  INPUT 


Upper:  Geometry  of  original  one-dimensional  Kubelka- 
Munk  (K-M)  theory 

Lower:  Geometry  of  6-flux  Modif ied-Kubelka-Munk 
(M-K-M)  theory. 
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Reflectance  (R|) 
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Viewgraph  VI:  Graphical  solutions  of  M-K-M  theory  for  incJex  of 
refraction  of  1.5;  Transmittances  (T)  and 
Reflectances  (R)  are  determined  on  powder  sample 
embedded  in  Potassium  Bromide  matrix: enter  graph 
with  T  _  and  R;  determine  6  and  r  ;  calculate 
scattering  (a)  and  absorption  (a)  from  lower 
equations . 


jTABLE  IV-2  Corning  CS1-64  Glass  Data  Analysis 
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Viewgraph  VII: 


Comparison  of  bulk  absorption  with  K-M  and  M-K-M 

samSVif  sintered  bubble  glLs 

?JsSitS.  CS1-64;M-K-M  model  gives  best 
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FASE/FASC0D4:  Current  Status  and  Future  Direction 


H.E.  Snell 

Atmospheric  and  Environmental  Research,  Inc. 

840  Memorial  Dr. 

Cambridge,  MA  02139 
phone:  (617)-547-6207 
fax:  (617)661-6479 
email:  snell@aer.com 

G.P.  Anderson,  J.H.  Chetwynd,  and  S.  Miller 
Geophysics  Directorate,  Phillips  Laboratory 
Hanscom  AFB,  MA  0173 1 
phone:  (617)377-2335 
email:  ganderson@plh.af.mil 

Abstract 

The  Optical  Physics  Division  of  the  Air  Force  Phillips  Laboratory  with  support 
from  the  Department  of  Energy  (DoE)  Atmospheric  Radiation  Measurement  (ARM) 
Program  is  developing  a  state-of-the-art  line-by-line  atmospheric  radiative  transfer  model 
as  the  successor  to  FASCODE  (Fast  Atmospheric  Transmittance  Code').  The  goal  of  this 
project  is  to  create  a  computationally  efficient  model  which  contains  the  most  up-to-date 
physics.  The  new  model,  known  as  FASCODE  for  the  Environment,  or  "EASE", 
combines  the  best  features  of  FASCODE  and  LBLRTM  (Line-by-Line  Radiative 
Transfer  Model^),  the  DoE's  standard  radiative  transfer  model.  Upgrades  to  EASE 
include  the  addition  of  a  solar  spectrum  module  to  compute  the  attenuated,  line-of-sight 
solar  radiance,  improvement  of  the  cloud  and  aerosol  descriptors,  based  on  changes  made 
to  MODTRAN  (Moderate  Resolution  Transmittance  Model),  and  the  ability  to 
incorporate  the  new  heavy  molecule  absorption  cross-section  data  found  on  the 
HITRAN96  database,  which  is  of  slightly  different  format  from  the  previous  HITRAN 
cross-section  data.  This  paper  focuses  on  the  new  heavy  molecule  cross-sections. 

Introduction 

Laboratory  measurements  of  the  absorption  properties  of  heavy  molecules  were 
recently  made  by  Varanasi  et  al.^  (subsequently  referred  to  as  "Varanasi's  cross  sections") 
and  are  included  on  the  1996  release  of  the  HITRAN  databases^  The  measurement  errors 
of  the  new  values  are  much  lower  than  the  old  data,  hence  the  inclusion  on  HITRAN-96. 
Because  of  the  complexity  of  these  molecules  and  the  fact  that  individual  spectral  lines 
cannot  be  resolved  at  temperatures  and  pressures  typical  of  the  earth's  atmosphere,  the 
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measurements  consist  of  relatively  low  spectral  resolution  information  in  the  form  of 
absorption  cross-sections  (i.e.  without  information  about  line  strengths,  positions, 
halfwidths,  etc.).  Many  of  these  molecules  have  been  included  on  previous  versions  of 
the  HITRAN  databases  in  that  form,  as  temperature  dependent  cross  sections  rather  than 
the  traditional  HITRAN  format  of  line  position,  line  strength,  halfwidth,  etc.  However, 
the  Varanasi  measurements  are  provided  as  a  function  of  pressure  and  temperature  as 
opposed  to  the  previous  HITRAN  format,  which  extrapolated  the  temperature-dependent 
values  to  zero  pressure.  To  compute  the  optical  depth  for  a  particular  (p,T)  combination 
with  the  'old'  data,  it  was  necessary  to  interpolate  between  temperatures  while  convolving 
with  the  appropriate  Lorentz  line  shape  for  the  pressure,  finally  combining  this  derived 
cross  section  with  the  density  weighted  path  amounts.  With  the  Varanasi  data,  one  is 
required  to  interpolate  and/or  extrapolate  between  the  tabulated  data  in  order  to  arrive  at 
the  correct  absorption  coefficient.  Thus  in  order  to  accommodate  this  new  data  directly  it 
would  be  necessary  to  re-configure  the  FASE  algorithms  to  bypass  the  convolution  and 
perform  the  necessary  interpolation  and/or  extrapolation.  Further,  one  must  develop  an 
appropriate  scheme  for  extrapolating  the  tabulated  values  to  pressures  and  temperatures 
outside  the  range  of  the  measurements. 

One  approach  through  which  the  Varanasi  cross-sections  may  be  utilized  by  the 
FASE  algorithms  is  to  perform  a  least-squares  fit  to  the  data  with  a  series  of  spectral 
lines,  solving  for  the  line  strength,  halfwidth,  and  lower  state  energy’.  The  pseudo-line 
approach  incorporates  all  the  (p,T)  combinations  of  data  into  the  least-squares  fit  and 
increases  the  accuracy  from  a  simple  interpolation  between  two  (p,T)  values.  The 
pseudo-line  data  is  easily  merged  with  the  spectral  data  taken  from  HITRAN  when 
creating  the  'TAPE3'  line  file,  and  the  radiative  transfer  calculations  are  performed  as 
usual,  with  the  model  using  the  spectral  line  routines,  rather  than  the  cross-section 
routines,  when  computing  the  absorption  of  these  species. 

One  drawback  to  the  pseudo-line  approach  is  that  new  sets  of  pseudo-line  fits  are 
required  as  more  molecules  are  added.  In  the  long  run  it  would  be  desirable  to  develop  a 
scheme  to  use  the  cross-sections  given  as  a  function  of  pressure  and  temperature.  Since 
the  problem  lies  in  interpolating/extrapolating  the  absorption  coefficient  data  to  specific 
values  of  the  pressure  and  temperature,  one  might  use  the  pseudo-line  approach  to  fill  out 
the  full  matrix  of  required  pressures  and  temperatures  (that  is,  use  FASE  to  compute  the 
absorption  coefficients  for  the  regions  not  represented  by  the  measurements).  For  the 
time  being  it  is  sufficient  to  investigate  the  implementation  of  the  pseudo-line  approach 
and  begin  to  explore  ways  in  which  the  actual  data  can  be  used,  perhaps  in  conjunction 
with  'pseudo-data'  based  on  calculations  with  the  pseudo-lines. 
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Pseudo-line  fits  have  been  made  available  for  FASE  for  a  number  of  heavy 
molecules  (see  Table  1).  However,  before  this  method  can  be  accepted  there  are  some 
questions  which  must  be  addressed:  (1)  does  the  calculation  with  the  pseudo-lines 
accurately  reproduce  the  measured  absorption  coefficients?  (2)  how  do  the  new  values 
for  heavy  molecule  absorption  impact  the  radiance  and  transmittance  computed  by 
FASE?  and  (3)  what  is  the  timing  impact  of  using  the  pseudo-lines  compared  to  the  (low 
spectral  resolution)  cross-sections?  An  evaluation  of  the  first  of  these  questions  was 
previously  presented  by  Snell  et  al.*  for  CFC-12  and  will  not  be  discussed  here.  The 
latter  two  questions  are  discussed  below. 


Table  1:  List  of  gases  with  pseudo-lines  currently  available  for  use  with  FASE.  Those 
species  with  cross-sections  given  on  the  HITRAN-96  database  are  indicated  in  the 
rightmost  column. 


Gas 

i  Interval  (cm"^) 

Grid  (cm'l) 

Measurer 

!  HITRAN-96 

CF4 

1265-1295 

0.005 

Varanasi 

CFC-12 

850-950  1 

0.01 

Varanasi 

X 

CFC-12 

1050-1200 

0.01 

Varanasi 

X 

CFC-11 

770-810 

0.01 

Varanasi 

X 

CFC-11 

1050-1120  1 

0.01 

Varanasi 

X 

CC14 

770-810  i 

0.01 

1  Orlando 

X 

CFC-22 

1  776-850 

0.00742 

1  Varanasi 

1  X 

CFC-113 

I  786-990 

0.01 

1  McDaniel 

SF6 

1  925-955 

0.01 

j  Varanasi 

1  X 

HCFC-142b 

1  870-1270 

0.01 

[  Newnham 

C10N02 

750-831 

0.00964 

Ballard 

X 

Radiance/Transmittance  Tests 

One  feature  of  the  new  heavy  molecule  cross  sections  is  that  the  measurements 
were  made  at  higher  spectral  resolution  and  with  a  better  signal-to-noise  ratio.  Thus  the 
new  values  exhibit  more  structure  and  a  larger  cross-section  at  band  center  than  the  cross- 
sections  tabulated  on  previous  versions  of  the  HITRAN  database.  While  the  heavy 
molecules  constitute  only  a  minor  fraction  of  gases  in  the  atmosphere,  the  impact  of  the 
new  values  on  different  calculations  must  be  determined.  The  test  scenario  was 
conducted  over  the  spectral  region  from  725  to  1245  cm'^.  The  atmosphere  consisted  of 
a  US  Standard  profile  with  the  first  seven  molecules  on  the  HITRAN  database  (those 


most  likely  to  contribute  to  a  typical  calculation)  along  with  HN03,  and  the  cross- 
sections  CFC-11  and  CFC-12.  The  CFC  profiles  were  consistent  with  those  obtained 
from  the  CIRRIS-1 A  experiment  (Bingham,  1997)  and  were  extrapolated  to  the  top  of  the 
atmosphere  using  the  shape  of  the  curve  given  in  Brasseur  and  Solomon  (1986). 

The  tests  were  conducted  for  input  atmospheres  with  (a)  no  CFCs,  (b)  CFCs  with 
the  old  cross-sections,  and  (c)  CFCs  with  the  new  cross-sections.  It  was  found  that  the 
new  cross-sections  have  a  negligible  impact  on  the  radiance  and  transmittance  for  zenith- 
and  nadir-looking  conditions.  There  are,  however,  noticeable  differences  for  certain  limb 
paths.  For  example.  Figure  1  shows  the  percent  difference  in  transmission  between 
calculations  using  the  old  and  new  values  for  the  cross-sections  for  a  zenith-looking  path 
from  the  surface  to  100  km.  One  can  see  that  the  differences  are  much  less  than  1%. 
Figure  2  shows  the  difference  in  transmission  with  a  tangent  path  from  100  km  to  10  km 
to  100  km.  The  region  around  1050  cm-1  is  opaque  due  to  ozone  absorption,  so  these 
differences  are  strictly  numerical.  However,  in  the  CFC  regions  (-800  -  950  cm”l  and 
-1070-1200  cm’l)  there  are  differences  of  the  order  of  2-5%  (though  the  differences  at 
the  center  of  the  CFC-12  band  reach  6%).  A  comparison  of  radiance  values  for  these 
types  of  paths  gives  similar  results.  The  above  tests  show  that  for  certain  paths  the  new 
values  of  the  heavy-molecule  absorption  cross-sections  give  answers  which  differ  from 
those  given  by  previous  versions  of  the  measurement  data.  Further  study  is  required  to 
assess  the  impact  of  these  new  values  in  trying  to  fit  measured  atmospheric  spectra. 
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Figure  1:  Differences  in  transmission  for  a  zenith  path. 


t 


Figure  2:  Differences  in  transmission  for  a  10  km  tangent  path. 
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Timing  Tests 

The  ability  to  read  the  new  cross-section  data  as  a  series  of  pseudo-lines  is  an 
inexpensive  alternative  to  developing  a  more  sophisticated  interpolation/extrapolation 
scheme.  However,  the  price  paid  is  one  of  algorithm  timing.  The  inclusion  of  pseudo¬ 
lines  can  add  anywhere  from  thousands  to  tens  of  thousands  of  spectral  lines  to  a  given 
calculation.  For  calculations  which  require  a  large  number  of  atmospheric  layers,  this  can 
add  a  significant  amount  of  time  to  the  calculation. 

A  series  of  tests  was  conducted  in  the  spectral  range  from  725  -  1245  cm‘*  using 
the  heavy  molecule  species  F-11,  F-12,  and  CC14.  These  species  were  chosen  because 
they  have  data  in  the  form  of  'old'  cross-sections  and  'new'  pseudo-lines.  Calculations 
were  done  for  the  uplooking  radiance  and  transmittance  from  0-40  km  and  0-100  km 
with  a  US  Standard  Atmosphere  (as  described  in  the  previous  section).  The  timing 
results,  given  in  Table  2,  represent  an  average  of  eight  individual  runs  for  each  case;  the 
standard  deviation  for  each  case  was  less  than  3  seconds. 
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The  calculations  from  100  km  to  the  surface  take  longer  than  those  from  40  km 
because  of  an  increased  number  of  atmospheric  layers,  and  because  the  upper  layers 
require  higher  spectral  resolution  than  the  lower  layers.  With  the  cross-sections, 
calculations  from  100  km  take  approximately  3.1  times  longer  than  those  from  40  km. 
The  difference  is  larger  with  the  pseudo-lines,  a  factor  of  about  3.5,  because  each  pseudo¬ 
line  must  be  accurately  computed  at  each  layer.  In  reality,  the  concentration  profiles  of 
many  species,  the  heavy  molecules  in  particular,  decrease  rapidly  with  altitude.  Using  the 
option  to  zero  the  contribution  from  weak  lines  when  they  no  longer  contribute 
significantly  to  the  overall  optical  depth  (NOZERO=0  on  input  RECORD  3.1) 
dramatically  decreases  the  computation  time  for  both  the  pseudo-line  and  cross-section 
cases  (however,  care  must  be  taken  for  calculations  which  require  very  high  accuracy, 
such  as  remote  sensing  retrieval  problems,  where  the  zeroing  of  absorber  amounts  can 
introduce  systematic  errors).  For  the  0-100  km  case  shown  here  there  is  a  20%  timing 
penalty  for  using  the  pseudo-lines  instead  of  the  cross-sections,  which  is  small  compared 
to  the  change  in  accuracy  with  this  new  data. 

Table  2;  Timing  test  results 


Case 

Old  Cross-Section 

New  Pseudo-Line 

0  -  40  km 

178.96 

401.11 

0  - 100  km 

562.93 

1413.07 

0  - 100  km 

263.02 

327.36 

(Zero  small  absorber 

amounts) 

Conclusion 

New  values  of  cross-sections  for  several  heavy-molecules  have  been  included  on 
the  1996  release  of  the  HITRAN  database  in  a  format  different  from  that  of  pervious 
versions.  A  pseudo-line  approach  was  adopted  whereby  "lines"  were  fit  to  the  cross- 
sectional  curves  to  allow  for  extrapolation/interpolation  to  temperatures  and  pressures  not 
given  on  HITRAN.  Depending  on  the  type  of  path  and  spectral  region,  the  new  values 
can  produce  significantly  different  values  (~5%)  for  the  calculated  transmission.  Tests 
have  shown  that  the  use  of  this  approach  can  add  a  significant  amount  of  time  to  the 
calculation.  Calculation  run-time  can  be  reduced  by  neglecting  the  calculation  at 
altitudes  where  the  species  of  interest  has  only  a  small  contribution  to  the  total 
radiance/transmission. 
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FASE-Calculated  Contribution  Functions  for 
MSX  SPiRiT-lii  MWIR  Bands  for  BTH  Measurements 
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Kaiser-Bessel  Apodization,  alpha=3.15 

Spectral  Resolution:  2.37  cm"'  Full  Width  Half  Height 


MSX  Radiometer  Filters  B1  &  B2  and  Interferometer  Scan  Function 
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Weighting,  Contribution  and  Sensitivity  Functions  Defined 
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LTE  Contribution  Function:  ZenithAngle  =  0 


o 


Watts/(cm*  sr '  cm  ’)/Km 


o 

o 

o 

o 

o 

t 

T— 

i 

-f 

o 

c 

o 

o 

o 

o 

o 

o 

q 

q 

o 

q 

<6 

ol 

d 

(•••subm) 


152 


2250  2300  2350  2400 

Wavenumber  (cm ') 


2300  2350  2400 

Wavenumber  (cm’) 
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Characterization  of  Optical  Extinction  Profiles  Using  Raman  Lidar  Techniques 

Michael  D.  O’Brien  and  C.  Russell  Philbrick 

Department  of  Electrical  Engineering 
Penn  State  University 


1.  Introduction 

The  presence  of  aerosol  particles  and  cloud  droplets  in  the  atmosphere  strongly 
affects  the  propagation  of  light  through  the  atmosphere.  These  particles  have  an  impor¬ 
tant  impact  on  the  Earth’s  climate  and  the  operation  of  electro-optical  sensors.  A  useful 
tool  for  measuring  optical  extinction  due  to  aerosols  is  a  Raman  lidar.  Raman  lidars  also 
have  the  capability  to  measure  atmospheric  temperature  and  humidity.  The  lidars  we  have 
at  Penn  State  also  have  the  capability  to  measure  ozone  in  the  lower  troposphere.  This 
profile  can  then  be  used  to  determine  the  extinction  due  to  ozone  absorption.  The  present 
lidar  system  provides  extinction  measurements  at  UV  and  visible  wavelengths.  The  tech¬ 
nique  has  been  verified  over  a  horizontal  path  with  independent  measurements  of  the 
particle  size  distribution  and  provides  measurements  with  a  high  degree  of  temporal  reso¬ 
lution. 

The  current  lidar  system  employs  Raman  scatter  from  a  ND;YAG  laser  beam  from 
the  2"*',  S''*,  and  4*’’  harmonics  at  532,  355,  and  266  nm  to  obtain  the  profiles  of  optical  ex¬ 
tinction.  The  measurements  have  been  made  using  the  molecular  nitrogen  profiles  at  607, 
387,  and  285  nm,  which  are  the  first  Stokes  vibrational  Raman  shifts  from  the  NDiYAG 
laser.  The  rotational  Raman  scatter  at  530  nm  has  also  been  used  to  obtain  optical  extinc¬ 
tion  profiles.  This  paper  vrill  begin  with  a  description  of  the  two  Raman  lidars  used  to 
acquire  the  data  to  be  presented.  Next,  the  details  of  our  analysis  will  be  discussed.  Then, 
the  capabilities  of  the  system  will  be  illustrated  with  several  examples  of  data  taken  during 
two  field  campaigns. 


2.  Background 

The  Raman  lidar  technique  is  a  powerful  way  of  measuring  a  variety  of  atmos¬ 
pheric  properties.  At  Penn  State  we  currently  have  two  Raman  lidars,  LAMP  and  LAPS. 
The  LAMP  system  began  operation  in  1991  and  has  been  used  to  make  a  variety  of  at¬ 
mospheric  measurements\  In  their  current  configuration,  both  LAMP  and  LAPS  are 
capable  of  measuring  water  vapor,  temperature,  ozone,  and  extinction  profiles  in  the  tro¬ 
posphere.  The  LAPS  lidar  system  is  a  second  generation  version  of  LAMP^.  The  LAPS 
system  has  a  much  greater  return  signal,  due  in  part  to  a  larger  receiving  telescope  and  a 
more  powerful  laser.  In  addition,  the  system  was  designed  as  an  operational  prototype 
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with  numerous  automated  systems.  This  makes  its  operation  easier  and  somewhat  more 
reliable  and  stable.  Table  1  summarizes  the  specifications  of  the  two  systems. 


Table  1.  Specifications  for  the  Penn  State  LAMP  and  LAPS  lidar  systems. 


LAMP 

LAPS 

Transmitter 

Continuum  NY-82  ND:YAG  laser 

20  Hz 

400  mJ  at  532  nm 

80  mJ 

Continuum  9040  ND:YAG  laser 

30  Hz 

600  mJ  at  532  nm 

120  mJ  at  266  nm 

Receiver 

0.41  m  diameter  telescope 

0.6  m  diameter  telescope 

Detector 

7  PMT  channels 

528,  530  nm  -  Temperature 

660,  607  nm  -  Visible  Water  Vapor 
295,  284  nm  -  UV  water  Vapor 

277,  284  nm  -  Ozone 

7  PMT  channels 

528,  530  nm  -  Temperature 

660,  607  nm  -  Visible  Water  Vapor 
295,  284  nm  -  UV  water  Vapor 

277,  284  nm  -  Ozone 

Data 

System 

100  MHz  count  rate 

75  m  range  resolution 

100  MHz  count  rate 

75  m  range  resolution 

For  both  of  the  systems,  water  vapor  profiles  are  determined  by  taking  the  ratio  of 
the  vibrational  Raman  return  from  water  vapor  to  the  vibrational  Raman  return  from  mo¬ 
lecular  nitrogen.  The  water  vapor  measurements  are  made  at  both  ultraviolet  and  visible 
wavelengths  using  the  295/284  and  660/607  ratios.  The  UV  measurement  is  corrected 
for  ozone  absorption  by  measuring  the  O2  vibrational  Raman  signal.  The  ratio  of  the  mo¬ 
lecular  oxygen  to  the  molecular  nitrogen  signal  provides  a  Raman/DIAL  measurement  of 
the  tropospheric  ozone  density^.  The  system  uses  the  temperature  dependence  of  rota¬ 
tional  Raman  scattering  to  obtain  profiles  of  atmospheric  temperature  using  the  ratio  of 
530/528  signals'*.  These  measurements  provide  useful  background  meteorological  data  for 
studying  changes  in  optical  extinction.  The  aerosol  extinction  coefficient  is  determined  by 
observing  departures  from  the  expected  gradient  of  the  Raman  lidar  profiles.  The  analysis 
used  for  this  method  will  be  presented  later. 

The  Raman  lidar  data  can  also  be  used  to  measure  the  optical  extinction  at  284, 
387,  530,  and  607  nm.  The  lidars  measure  the  scattered  light  corresponding  to  the  1st 
Stokes  vibrational  Raman  transition  for  molecular  nitrogen  to  obtain  a  profile  of  atmos¬ 
pheric  nitrogen.  Any  deviation  in  this  profile  jfrom  the  gradient  of  the  molecular 
atmosphere  is  due  to  aerosol  extinction.  Figure  1  shows  a  schematic  representation  of  this. 
The  figure  shows  the  range  corrected  signals  for  a  Raman  shifted  nitrogen  channel  and  a 
corresponding  elastic  backscatter  channel.  The  figure  demonstrates  when  there  is  a  large 
increase  in  the  backscatter  there  is  a  deviation  in  the  slope  of  the  molecular  nitrogen  Ra¬ 
man  signal.  We  are  trying  to  measure  the  change  in  this  slope,  which  is  proportional  to 
the  exponential  attenuation  term  in  the  Raman  lidar  equation.  Figure  2  shows  an  example 
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of  the  raw  photon  count  profiles  from  LAPS  on  9  October  1996  while  it  was  aboard  the 
U.S.N.S.  Sumner.  The  profiles  rapidly  increase  to  a  maximum  as  the  laser  fully  enters  the 
field  of  view  of  the  telescope  at  approximately  1000  m  ^  Above  this  altitude  the  profile 
shows  an  altitude  region  where  the  slope  of  the  profile  is  much  greater  than  the  expected 
gradient  for  the  molecular  atmosphere.  This  increased  slope  is  caused  by  an  increase  in 
the  optical  extinction  due  to  scattering  particles. 


Figure  1 .  This  figure  depicts  the  change  in  the  slope  of  the  Raman  signal  as  the  lidar  hits 
an  intense  scattering  layer. 
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LAPS  Raw  Data 


Accumulated  Counts 


Figure  2.  This  figure  shows  a  plot  of  the  raw  data  fi’om  LAPS  taken  on  9  October  1996  at 
03:16  GMT  aboard  the  U.S.N.S.  Sumner.  The  profiles  show  a  region  of  heavy  aerosol 
extinction. 


3.  Analytical  Equations 

It  is  possible  to  derive  a  relation  for  the  extinction  coefficient  fi’om  the  Raman  lidar 
equation  The  Raman  lidar  equation  is  given  as, 

P(2j  =  ^N(z}^ex,\^-f{ar'(z)  +  a^{2)  +  arU)  +  a--'(z))dzj  ,  (1) 

where  P(z)  is  the  received  signal  power  at  the  Raman  shifted  wavelength,  k  is  an  instru¬ 
mental  constant,  ^(z)  is  the  telescope  form  factor,  N(z)  is  the  molecular  number  density, 
doldCl  is  the  Raman  backscattering  cross  section,  and  are  the  molecular  extinc¬ 
tion  coefficients  for  the  transmitted  and  Raman  scattered  wavelengths,  and  are 
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the  aerosol  extinction  coefficients  for  the  transmitted  wavelength  and  Raman  scattered 
wavelengths,  and  z  is  the  range.  The  aerosol  extinction  coefficient  can  be  determined 
from  (1)  and  is  given  by, 


aV'  = 


dz 


In 


^n(^) 


-aT(z)-aT(z)-aT(z) 


(2) 


Only  data  above  the  altitude  where  the  laser  is  completely  in  the  field  of  view  of  the  tele¬ 
scope  is  used  in  this  analysis.  The  molecular  number  density  and  extinction  coefficients 
are  determined  from  a  standard  atmosphere  and  the  known  molecular  scattering  cross 
section.  The  aerosol  extinction  coefficient  for  the  transmitted  wavelength  of  532  nm  is 
determined  from  the  rotational  Raman  signal.  The  received  signal  for  the  rotational  Ra¬ 
man  channels  is  also  given  by  (1),  and  in  this  case  the  extinction  coefficients  of  the 
transmitted  and  backscattered  Raman  wavelengths  are  assumed  to  be  the  same  because 
they  are  only  separated  by  2  nm.  The  532  nm  aerosol  extinction  coefficient  is  then  deter¬ 
mined  with, 


a, 


'532 


dz 


m 


(3) 


This  value  of  then  used  as  the  value  of  in  (2)  so  that  can  be  calculated. 

This  procedure  assumes  that  the  optical  depth  is  not  too  great  and  that  multiple  scattering 
can  be  neglected. 

At  UV  wavelengths  a  rotational  Raman  signal  is  not  available  so  only  the  vibra¬ 
tional  Raman  signals  are  used.  The  extinction  due  to  ozone  absorption  is  determined  with 
the  ozone  profile  measured  by  the  lidar.  The  difference  between  the  extinction  coefficients 
between  the  transmitted  and  received  wavelengths  is  accounted  for  with  a  power  law  rela¬ 
tionship,  which  assumes  the  same  wavelength  dependence  for  scattering  as  was  found 
between  532  and  607  nm.  The  resulting  equation  is  given  as, 

i  [‘"wl  -  i< + O  -  «  + )  ■  No.  (z) 


where  cr^/g  and  <72°4  are  the  ozone  absorption  cross  sections  for  the  transmitted  and  Ra¬ 
man  scattered  wavelengths,  (z)  is  the  ozone  density  profile,  and  k  is  determined  by 
the  wavelength  dependence  of  particle  scattering. 
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4.  Results 


Over  the  past  years  our  research  group  has  gathered  data  with  the  two  lidar  sys¬ 
tems  under  a  variety  of  conditions  and  locations.  The  LAMP  lidar  was  used  in  September 
1995  to  measure  the  particle  size  distributions  and  extinction  properties  of  coastal  aerosols 
in  the  Coastal  Aerosol  Scattering  Experiment  I  (CASE  I)  campaign  conducted  at  NASA’s 
Wallops  Island  test  facility.  The  lidar  was  operated  in  vertical  and  horizontal  modes.  In 
the  vertical  operating  mode,  the  lidar  measured  water  vapor,  temperature,  and  ozone  pro¬ 
files  simultaneously  with  other  lidar  and  balloon  sensors^®.  In  the  horizontal  operating 
mode,  the  aerosol  extinction  coefficient  was  determined  from  the  Raman  lidar  profiles  and 
from  particle  size  measurements  made  with  a  bistatic  lidar.  The  bistatic  lidar  consists  of  a 
camera  that  focuses  an  image  of  the  horizontal  laser  beam  onto  a  linear  photodiode  array. 
This  allows  for  the  measurement  of  the  scattering  intensity  fi'om  the  aerosol  particles  at  a 
number  of  angles.  Data  from  the  ratio  of  the  cross  polarization  scattering  phase  function  is 
inverted  to  determine  the  particle  size  distribution^*.  The  aerosol  extinction  coefficient 
that  corresponds  to  this  distribution  is  then  calculated.  A  comparison  of  the  extinction 
values  determined  with  the  two  techniques  is  given  in  Table  2.  The  two  extinction  coeffi¬ 
cients  agree  within  the  experimental  error  for  the  methods. 


Table  2.  Comparison  of  aerosol  extinction  coefficients  measured  with  a  Raman  lidar  and  a 
bistatic  lidar. 


Wavelength 

Extinction  fi'om  Raman  Lidar 

Extinction  from  Bistatic  Lidar 

532  nm 

0.259  (km‘*) 

0.279  (km'*) 

607  nm 

0.209  (km'*) 

0.243  (km'*) 

During  the  Fall  of  1996  LAPS  was  deployed  aboard  the  USNS  Sumner  for  its  first 
field  tests.  The  system  performed  quite  well  and  provided  several  hundred  hours  of  data. 
One  example  of  an  extinction  profile  from  this  data  is  shown  in  Figure  3.  This  profile  of 
optical  extinction  for  532  ran  was  obtained  on  9  October  1996.  The  profile  shows  a  sev¬ 
eral  layers  of  strong  extinction.  At  the  same  time  we  determined  the  extinction  coefficient 
at  284  nm,  which  is  shown  in  Figure  4.  With  the  LAPS  lidar  system  we  have  the  capabil¬ 
ity  to  make  measurements  of  the  extinction  coefficient  with  high  time  resolution.  An 
example  of  this  capability  is  shown  in  Figure  5. 
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Figure  3.  This  Figure  shows  the  optical  extinction  for  532  nm  measured  vertically  with 
LAPS  on  9  October  1996  at  04:16  GMT.  The  data  was  acquired  with  a  30  minute  inte¬ 
gration.  The  molecular  extinction  profile  is  shown  for  reference. 


Extinction  Coefficient  (km'^) 


Figure  4.  This  Figure  shows  the  optical  extinction  at  284  nm  measured  vertically  with 
LAPS  on  9  October  1996  at  04:16  GMT.  The  data  was  acquired  with  a  30  minute  inte¬ 
gration.  The  profiles  left  to  right  are  ozone  extinction,  molecular  extinction,  and  aerosol 
extinction. 
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Figure  5.  This  Figure  shows  a  time  series  of  the  532  nm  extinction  coefficient  starting  at 
03:00  EDT  on  9  October  1996  from  LAPS  aboard  the  U.S.N.S.  Sumner.  The  intensity 
scale  is  in  km'^ . 


4.  CONCLUSIONS 

The  measurements  presented  in  this  paper  have  demonstrated  the  utility  of  the 
Raman  lidar  technique  in  measuring  optical  extinction.  These  measurements  have  been 
independently  verified  with  the  bistatic  lidar  measurements.  The  vertical  profiles  showed 
that  this  technique  can  be  successfully  employed  to  measure  extinction  coefficients  in  the 
lower  troposphere  for  aerosols.  In  addition  we  have  shown  that  the  time  resolved  meas¬ 
urements  of  extinction  provide  a  powerful  way  of  studying  aerosols  in  the  atmosphere. 
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CIRRIS-IA  limb  spectral  measurements  on  mesospheric  ozone  airglow 
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Abstract.  Shuttle-based  limb  infrared  spectral  measurements  on  mesospheric  ozone  airglow  were  made  by 
the  Cryogenic  Infrared  Radiance  Instrumentation  for  Shuttle  (CIRRIS-IA)  during  one  flight  of  the  space 
shuttle  Discovery  (28  April  to  6  May,  1991).  A  significant  difference  between  ozone  nightglow  and  day- 
glow  was  observed.  A  spectrally  integrated  ozone  U3  fundamental  band  radiance  profile  was  used  to  re¬ 
trieve  its  excited  state  density  which,  in  conjunction  with  the  modeled  ozone  U3  vibrational  temperature, 
was  used  to  infer  the  total  ozone  density.  The  night/day  ozone  profiles  demonstrate  evidence  of  a  non¬ 
local  thermodynamic  equilibrium  mesospheric  ozone  diurnal  variation. 


1.  Introduction 

Measurements  of  the  mesospheric  ozone  airglow  were 
first  made  by  the  Cryogenic  Circular-Variable  Filter  Spec¬ 
trometer,  a  rocket-borne  experiment  [Stair  et  al.,  1974]. 
Later,  more  advanced  sensors  were  used  to  observe  meso¬ 
spheric  ozone  [e.g.,  Stair  et  al.,  1985;  Rawlins  et  aL, 
1993].  Some  laboratory  studies  have  been  made  to  im¬ 
prove  our  understanding  of  ozone  kinetics  and  spectros¬ 
copy  [e.g.,  Menard  et  ai,  1992;  Upschulte  et  al.,  1994].  It 
is  well  known  that  ozone  is  produced  by  the  highly  exo¬ 
thermic  recombination  (O+O2+M  O3+M),  and  mainly 

destroyed  by  the  photolytic  reaction  (03+hv  O+O2).  A 

few  non-local  thermodynamic  equilibrium  (non-LTE) 
ozone  models  have  been  developed  based  on  the  knowl¬ 
edge  of  ozone  kinetics  and  spectroscopy  [e.g.,  Mlynczak 
and  Drayson  1990a,  b;  Manuilova  and  Shved,  1992; 
Sharma  et  al.,  1996].  Laboratory  studies,  model  develop¬ 
ments,  and  experimental  measurements  have  sought  to 
gain  the  information  necessary  to  understand  the  middle 
atmospheric  ozone  behavior  and  its  role  in  Earth’s  atmos¬ 
pheric  environment. 

This  letter  reports  on  the  shuttle-based  infrared  spectral 
measurements  of  the  Cryogenic  Infrared  Radiance  Instru¬ 
mentation  for  Shuttle  (CIRRIS-IA)  which  was  operated 
on-board  the  space  shuttle  Discovery  (28  April  to  6  May, 
1991).  This  broadband  interferometer  with  moderated 
spectral  resolution  (1,0  cm'^)  provided  the  data  needed  to 
analyze  spectrally  radiant  contributions  from  the  different 
gases  and  their  emission  bands;  and  the  limb  scan  profile 
allowed  the  retrieval  of  atmospheric  gas  densities.  Meso¬ 
spheric  ozone  night/day  measurements  are  reported  in  sec¬ 
tion  2.  The  retrieved  ozone  U3  excited  state  population  and 
inferred  total  ozone  abundance  are  presented  in  section  3, 
and  the  summary  and  remarks  are  given  in  section  4, 

2.  Spectral  Analysis  /  Data  Preparation 


The  CIRRIS-IA  data  presented  here  were  collected  on 
29  April  1991  by  a  Michelson  interferometer  with  detectors, 
preamplifiers,  and  collection  optics  (0.3-m  telescope) 
cooled  to  near-liquid-helium  temperatures.  The  CIRRIS-IA 
instrument,  measurements,  calibration,  tangent  height  de¬ 
termination,  and  other  details  of  this  experiment  are  docu¬ 
mented  elsewhere  [e.g.,  Ahmadjian  et  al.,  1990;  Bingham  et 
al,  1997].  Four  nighttime  limb  scan  profiles  and  two  day¬ 
time  limb  scan  profiles  were  used  to  provide  one  mean  pro¬ 
file  for  nighttime  and  one  mean  profile  for  daytime.  The 
data  collection  location,  solar  zenith  angle  (SZA),  and  local 
time  (LT)  are  listed  in  Table  1. 


Table  1.  CIRRIS-IA  Limb  Scan* 


Lat. 

Long. 

SZA 

LT 

60°N 

185®E 

490 

1342 

59°N 

168°E 

50® 

1419 

28°S 

353°E 

152® 

0142 

29°S 

13°E 

153® 

0132 

44°S 

25°E 

147® 

2249 

43°S 

316°E 

149° 

2304 

^Indicated  near  85  km. 


The  samples  of  spectra  (1010  -  1 140  cm’^)  obtained  from 
limb  observations  under  nighttime  and  daytime  conditions 
(near  85  km)  are  shown  in  Figures  la  and  lb,  respectively. 
The  estimated  sensor  noise  of  -2x10'*’  (W  cm‘^  srVcm"*) 
near  10  pm  through  the  mesosphere,  yields  a  signal  to  noise 
level  high  enough  to  identify  the  O3  Uj  band  emissions  at 
tangent  heights  below  ~90  km.  The  spectra  are  primarily 
generated  by  the  emissions:  (1)  O3  U3  and  Uj  fundamental 
bands  and  the  wealth  of  O3  hot  bands,  and  (2)  CO2  (0001 1- 
10002)  band  mostly  at  daytime.  The  O3  O3  emission  band 
contribution  was  dominated  in  this  region,  with  the  O3  Uj 
band  and  hot  bands,  as  well  as  the  CO2  band  making  minor 
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but  significant  contributions.  The  shaded  regions  shown  in 
Figure  1  are  the  Strategic  High-Altitude  Radiance  Code 
(SHARC)  [Sharma  et  al,  1996]  modeled  spectra  fit  to  the 
CIRRIS-IA  data. 


(a)  NIGHTTIME  SPECTRA 
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(b)  DAYTIME  SPECTRA 
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Figure  1.  Examples  of  CIRRIS-IA  mesospheric  ozone  spectra 
collected  at  (a)  nighttime  and  (b)  daytime.  The  spectral  contri¬ 
butions  from  different  bands  were  modeled  by  SHARC-3  model 
to  derive  the  band  radiance  profiles. 

SHARC  modeled  spectra  were  added  band  to  band  to 
fit  the  measurement.  In  the  nighttime  profiles,  the  CO2 
contribution  is  negligible  while  the  O3  U3  band  takes 
-91.5%  of  the  total  band  radiance  with  the  standard  devia¬ 
tions  of  -1.6%  through  tangent  heights  of  70  -  100  km. 
The  same  analysis  was  performed  for  the  daytime  profiles. 
The  averaged  contribution  fractions  of  the  CO2  (00011- 
10002)  band  and  O3  (001-000)  band  were  derived  using 
two  daytime  profiles  collected  under  almost  the  same  con¬ 
ditions  (SZA  of  -50®)  as  shown  in  Figure  2.  As  a  result, 
the  spectrally  integrated  band  radiance  data  were  averaged 
to  one  mean  nighttime  profile  and  one  mean  daytime  pro¬ 
file  as  shown  in  Figure  3.  Because  the  O3  U3  band  radiance 


has  a  higher  signal  to  noise  ratio  and  the  band  strength  is 
less  thermally  dependent  (-7%  from  296  K  to  200  K),  the 
spectrally  integrated  O3  U3  radiance  mean  profile  is  chosen 
to  invert  the  O3 1)3  excited  state  population  which,  in  con¬ 
junction  with  the  modeled  vibrational  temperature,  is  used 
to  infer  the  total  O3  abundance  profiles  in  the  next  section. 


Contribution  Fraction  (%) 


Figure  2.  Daytime  O3  (001-000)  and  CO2  (00011-10002)  band 
radiance  contributions  to  the  spectral  region  of  1010-1 140  cm‘l. 


Figure  3.  Nighttime  and  daytime  CIRRlS-lA  measured  band  ra¬ 
diance  profiles  with  the  mean  profiles  and  O3  03  band  radiance 
profiles. 


3.  Results 

The  retrieval  can  be  processed  using  the  line-by-line 
(LBL)  method  or  the  weak  line  (WL)  approximation  [e.g., 
Zachor  and  Sharma,  1985;  Solomon  et  al,  1986].  Before 
the  inversion,  both  the  LBL  and  WL  approximation  forward 
O3  U3  band  radiance  calculations  were  made  and  the  dis¬ 
crepancy  between  these  two  methods  was  small  (<  5%  at 
tangent  heights  above  70  km).  Therefore  the  WL  approxi¬ 
mation  was  used  in  this  work,  Rooi 
(l/47i)Z{hcVio[03(001)]A,o}dx,  where  Rooi  is  the  measured 
O3  O3  band  radiance,  h  is  the  Planck’s  constant,  c  is  the 
speed  of  light,  Vjq  is  the  frequency  of  the  transition  in  wave- 
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number,  [03(001)]  is  the  ozone  03  excited  state  population, 
A 10  is  the  Einstein  coefficient,  and  dx  is  the  elements  of 
limb  measurement  path  length  along  the  line  of  sight. 
However,  the  thermal  correction  and  bandpass  correction 
on  the  Einstein  coefficient  were  made  through  the  band 
strength  correction  (using  HITRAN  database  {Rothman  et 
al,  1992]).  Now,  the  [03(001)]  retrieval  is  an  Abel-type 
inversion.  The  retrieved  [03(001)]  profiles  from  nighttime 
and  daytime  are  shown  in  Figure  4. 


E 

J!*J 


D) 
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Figure  4.  Retrieved  ozone  nighttime  and  daytime  03  excited 
state  population  associated  with  03  radiance  shown  in  Figure  3. 

To  accurately  retrieve  O3  abundance  depends  on  the 
knowledge  of  O3  kinetics  and  spectroscopy  which  yield  the 
ratio  of  the  O3  excited  state  to  ground  vibrational  popula¬ 
tion  (i.e.,  the  vibrational  temperature).  The  chemical 
pumping  fraction  of  O3  recombination,  the  collisional 
quenching  rates,  etc.  become  the  critical  parameters  of  the 
O3  non-LTE  model  to  precisely  predict  vibrational  tem¬ 
peratures  {Mlynczak  et  al.^  1996].  Therefore,  the  following 
results,  viewed  as  an  estimation,  are  based  on  the  models 
using  current  knowledge  of  O3  kinetics  and  spectroscopy. 

Estimation  of  the  O3  ground  state  population  is  made 
utilizing  two  models:  (1)  the  SHARC-3  model  and  (2)  the 
Simplified  Single  Mode  (SSM)  non-LTE  O3  model 
[Mlynczak  and  Drayson,  1990a,  b].  The  model  atmos¬ 
phere  for  CIRRIS-IA  observation  was  obtained  from  the 
2-dimensional  (2-D)  photochemical  model  of  Garcia  and 
Solomon  [1983,  1985].  The  different  model  inputs  (e.g., 
the  quenching  rates)  largely  affect  the  vibrational  tem¬ 
perature  calculations.  For  example,  the  temperature  inde¬ 
pendent  quenching  rate  was  used  in  the  SSM  model  while 
a  Landau-Teller  type  temperature  dependent  quenching 
rate  was  used  in  the  SHARC  model.  It  was  found  that  the 
quenching  rates  initially  used  in  the  SSM  model  were 
larger  in  comparison  with  recent  laboratory  measurements 
[e.g.,  Menard  et  aL  1992;  Upschulte  et  al.,  1994].  The 
following  results,  related  to  the  SSM  model,  were  obtained 


using  the  modified  quenching  rates  (scaled  down  by  a  factor 
of  4). 

The  total  inferred  O3  abundance  is  dependent  on  the 
model  predicted  vibrational  temperature  which  is  shown  in 
Figure  5.  The  population  ratio  written  as  a  function  of  the 
vibrational  temperature  is  [03(00 1)]/[03]  =  exp(-C2Vio/Tooi), 
where  C2  is  the  second  radiation  constant  and  Tqoi  is  the  vi¬ 
brational  temperature.  Using  the  Tqoi  derived  from  the  ap¬ 
plications  of  SHARC  and  SSM  O3  models  to  2-D  modeled 
atmosphere,  the  ground  state  number  density  [O3]  profile 
was  readily  inferred  using  the  above  equation  which  applies 
retrieved  [03(001)]  data.  The  [O3]  profiles  resulting  from 
the  SHARC  and  SSM  modeled  Tqoi  and  the  [03(001)]  re¬ 
trieved  from  the  measurements,  including  the  discrepancies 
[([OslsHARc-  [OsIssm)  /  [OslsHARc],  ^Te  plotted  in  Figure  6. 


Temperature  (K) 


Temperature  (K) 


Figure  5.  The  2-D  photochemical  model  derived  local  kinetic 
temperature,  and  ozone  03  vibrational  temperatures  modeled  by 
(a)  SHARC-3  and  (b)  SSM  ozone  models. 
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Figure  6.  Inferred  ozone  density  profiles  using  CIRRIS-IA  O3 
03  excited  state  population  profiles  in  conjunction  with  modeled 
03  vibrational  temperatures  calculated  by  the  SHARC  and  SSM 
models.  The  discrepancy  from  the  SHARC  and  SSM  models, 
([OslSHARC  -  [OsISSM)  >  [OslSHARC^  was  also  plotted. 


167 


4 


4.  Summary  /  Remarks 

Mesospheric  O3  infrared  emission  data  from  CIRRIS- 
lA  shuttle-based  interferometer  limb  measurements  were 
analyzed.  This  high  quality  spectra  data  set  allowed  de¬ 
termination  of  the  O3  U3  and  Oj  fundamental  bands  and 
CO2  band  contributions.  The  limb  profile  measurement 
allowed  the  retrieval  of  the  O3  O3  excited  state  population 
used  to  infer  the  total  ozone  abundance  which  is  partially 
based  on  the  current  knowledge  of  ozone  formation 
(chemical  pumping  process),  quenching,  absorption,  and 
other  elements,  as  well  as  the  modeled  atmosphere.  Accu¬ 
rate  laboratory  measurements  of  ozone  kinetic  parameters 
will  be  beneficial  in  the  more  precise  modeling  of  non- 
LTE  ozone  and  hence  accurately  retrieve  mesospheric 
ozone  measurements. 

This  set  of  data  is  deemed  of  interest  because  the  en¬ 
ergy  loss  rate  through  the  O3  infrared  radiation  (860  -  1 160 
cm'*)  can  be  retrieved;  the  energy  absorption  gain  rate 
through  the  O3  U3  band,  the  chemical  potential  energy 
deposition  rate,  and  the  solar  energy  deposition  rate 
through  the  O3  Hartley  band  can  be  evaluated  based  on  the 
ozone  density  taken  from  CIRRIS-IA  measurements.  A 
more  complete  communication  which  addresses  some  of 
these  issues  is  in  preparation. 
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MIDCOURSE  SPACE  EXPERIMENT  (MSX) 
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SPIRIT  III  FIELDS  OF  VIEW 


THE  INTERFEROMETER’S  SIX  CHANNELS 
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For  tangent  heights  near  60  km 

Lower  value  corresponds  to  saturation  (max  S/N) ,  upper  one  to  to  max[spectral  S/N]  =  10 
(based  on  preflight  engineering  data) 

“Open”  =  “no  filter” 

“Prewhitened”  =  notch  attenuation  (1 0'^  over  the  1 5  pm  COj  band 


o 


< 

Q 

3 

m 

m 


CO  □: 
CO  h- 
LU  cc 
O  < 
O  ^ 

PF  ^ 

CL  111 
I — 
Lll 


Q 

Lll 

o 

h- 

D 

< 


O 

cn 

Lll 

LL 

cc 

LU 

h- 

z 

LU 

X 

K 

LL 

O 


Co 


e 

St 

'S 

St 

o 

D 

Si 


Co 

bjo 

•S 

S 

<3 


CN 

> 

o 


c/3 


03 

i-t 

bJ) 


-o 

<D 

-4— > 

o 

U4 

o 

o 


03 


03 


03 

<N 


t  t 


CL 


u 


T 


^  t  t 


S 

Q 

O 

I 

I 

CO 


CO 

< 

pci 

O 

O 

P< 

tL 

P< 

tu 


(N 

> 

cu 


Vi 

<o 

o 

C3 

c3 

Vi  ’.i 

4)  " 

o  CO 

a  ^ 

•2  2 

^  iJ 

2  o 

Sl-o  S- 


Pi 

o 


S 

CQ  a> 

.fi  S 

S  bO 

CA  O 

r  a 

6o 

(U 


'H- 

S 

CD 

^  C/5 

o  8 

i-i 

I 

«-w  ^ 

O  ^ 

®  ^ 
^  X) 

CQ 

V5 


"d 

O) 


V5 


C/3 

a 

o 

^  8 
•ri  > 
T3  CD 
(U 

H  o 


a>  ci-i  c/5 


< 

Pi 


PQ 

Pm 

CO 


^  ^  •T-i 

a  ,2  g 

S  rs  T-j  cJ 


SS  o-’o 

s  g  i 

S 


cd 


cd 


O 


173 


OVERVIEW  OF  AUTOMATED  ANALYSIS  AND  INTERACTIVE  GRAPHICAL  REVIEW 


174 


E-8 


CONSTRUCTION  OF  TEMPLATE  SPECTRA 
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TEMPLATE  SPECTRA  CONSTRUCTED  IN  THE  AUTOMATED  ANALYSIS 
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W AVENUM  BER 


DCATT  EXPERIMENT  DC3701 00005  ON  18  SEPTEMBER  1996 
SERIES  OF  CONSTANT  TANGENT  HEIGHTS;  LONG-STROKE  MODE  (2  cm  MAX  RESL  N)  - 


'4^02-000) 


CHANNEL-6  SPECTRA  FROM  DC37  ON  9/18/96 
H_  =  90  km,  Lat/Lon  =  4N/163W 
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CHANNEL-1  SPECTRA  FROM  DC37  ON  9/18/96 
=  80.6  km,  Lat/Lon  =  7S/1 68W 
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EARTHLIMB  EXPERIMENT  EL090200050  ON  18  OCTOBER  1996 
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GMS-5:  8:32  on  October  18, 1996 
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-  FIRST  SIDELOBE  IS  2.2  x  10''  (10^  RELATIVE  TO  SING  FUNCTION) 
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EXAMPLES  OF  THE  USE  OF  “INDICATOR”  HISTORIES 
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R  vs.  is  compared  to  a  threshold  value  that  indicates  03  depletion. 
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RED  EDGE  OF  FILTER  (^m) 


AN  INDICATOR  FOR  OZONE  DEPLETION  ("OZONE  HOLF;;) 


3  SO  100  ISO  ^  OBSERVED/PREDICTED  AMBIENT  (no  hole) 

OZONE  PARTIAL  PRESSURE  (nb)  SPECTRAL  RADIANCE  AT  12.85  ixm 


Pure  Rotation  Emissions  From 
Mesospheric  OH  Observed  Using  the 
SPIRIT  III  Interferometer 
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Sponsored  by  the  Ballistic  Missile  Defense  Organization  (BMDO) 


MSX  Earthlimb  Measurement: 
October  8, 1996 
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OH  Pure  Rotation  Emissions 
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Vibrational  Energy  (cm 
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Latitude  (deg.  N) 


Latitude  Dependence: 

Rotationally  Thermal  Population 
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Latitude  (deg.  N) 
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Strong  latitude  dependence  of  the  high-N  population 
Implications  for  odd  oxygen  chemistry  models 
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High  Latitude  (Auroral)  Scenes 
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enhanced  during  particle  and  Joule  heating  events. 


Low  to  Mid-Latitude  Scenes 
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EL06:  CONSTRNT  TRNGENT  HEIGHT  SCPN 
EL06G 1000 160 10101  RS-FLOWN  ON  09/14/96 
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depletions  in  15  um  C02.  This  is  probably  due  to  depletions  in  O 
during  geomagnetic  activity. 


Radiometer  Band  D  Radiance  Profile 
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Radiometef  Band  D  Radiance  Profile 
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Hext'  to  Conti  HUB, 
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Press  'Next' to  Continue, 


Radiometer  Band  B  Radiance  Profile 
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MSX:  Middle-Atmospheric  Structure 
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SX  structure  -  Problem  &  Approach 
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Higher-altitude  structure  is  consistent  with  origin  in 
internal  gravity  waves  (GW) 


MSX  Structure  -  Observations 


q 

h- 

< 

> 

(T 

LU 

CO 

QQ 

O 


■a 


13 

o 

0 

0 

0 

0 

o 

0 

•  mmmm 

"0 

•  ■MM 

0 

m  imtmm 

0 

JZ 

> 

> 

> 

> 

CD 

0 

0 

0 

0 

■  mammm 

I 

CD 

CO 

CM 

CM 

O 

•  m 

o 

M  W 

9^ 

O 

CO 

00 

03 

■  ■ 

o 

o 

T- 

o 

CM 

CD 

o 

CD 

CO 

CO 

LO 

'T— 

CM 

CM 

CM 

CO 

a: 

LU 

00 

O) 

00 

CJ3 

C33 

MhI 

Q. 

o 

1— 

o 

O 

X 

— J 

-J 

-J 

-J 

LU 

LU 

LU 

LU 

LU 

LU 

CD 

CD 

CD 

CD 

03 

(33 

<33 

CD 

C33 

03 

CL 

C33 

> 

0 

13 

0 

■+— • 

O 

LU 

< 

CO 

O 

Z 

< 

lO 

r-. 

o 

U 

CO 

D 

CM 

T“ 

CM 

• 

• 

• 

• 

• 

214 


^Experiments:  EL09  -  ATH  /  BTH  Clutter 

EL18  -  BTH  Clutter 


EL09  17  Auf»usi  1996 


tLl:19:32l7 


SOA  <  0 
SOAO-T2 
SOA12-1® 
SDA  »  10 


i7iM« 


ELI 8  25  May  1996  (part  3) 


)8:20:0. 3 


/••f--** — 1^03:119:0.3 


kl8:0  J3 


08:11:0.3 


08:16:0. 


105  107  109  111  ■  113  115  117  j2i  123  125 

- - - - - Longitude _ _ _ _ 

Hgurc  l  The  gromd  trace  of  the  MSX  spacecraft  and  the  target  point  at  1  minute  intervals'.  The  upper 
figure  IS  for  the  17  August  1996EL09  (nadir  angle  57.5  deg.)  data  coUection  evenn  In  ±e  lower  fi^ 
e  satellite  ^und  trace  and  the  target  pomt  trace  are  nearly  coincident  for  the  25  May  1996  ELi8 
(naoir  angle  1.5  deg.)  measurement. 
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October  2, 1996  -  NIGHT  (EL090200047010101) 

174  SECOND  SCENE  (21: 1:18  -  21:  4:13  UT)  .Nadir  =  60.5 
BAND  B1:  4.22  -  4.36  uirCAND  B2:  4.24  -  4.45  urn 
Vertical  Cross-Section  Lat:  26.5  -  17.1;  Lon:  90.1  -  89.4 

at  2:31Min;Sec  437km  X  1047km  -  IFOV2279X  253  Meters  ^  HO 
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MSX  Structure  -  4.3-\xm  Radiative 

Transfer 
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LTE  and  Non-LTE  Radiative  Transfer 
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-  radiative  loss  of  CO2  V3 

-Integrate  optically  thick/thin  4.3-pm  CO2  V3  lines 
line-of-sight 


The  LTE  and  Non-LTE  Radiative  Transfer 

Probiem 
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=  integrated  mean  intensity  =  average  of  /(v)  over  absorption  iine  and  over  soiid  angie 


SAMM  (SHARC  And  MODTRAN 

Merged) 
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aerosol  single  /  multiple  scattering,  single-scatter  solar 
radiance 

Path  averages  by  Curtis-Godson  approximation 


ARC  (Atmospheric  Radiance  Code) 
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Calculate  line-of-sight  radiance  spectra  including  opacity 
effects 


Production  &  Loss  Processes  -  002(^3) 
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AE=1 9.2  cm'^:  kw  =  5.0  x  10'i'^(300/Tr 


Production  &  Loss  Processes  -  N2(v) 
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-W  narrows 

-tropospheric  peak  disappears 


aaruLLTv 
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MSX  Structure  - 
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In-scan  PSDs  confirm  this  prediction 


O' 

G) 

O 

o 

o 


c 

X 


n 

O 

u 

o 


.n] 

C/^ 

a 


© 

«M 

««( 

0) 

o\ 

0 

CO 

X 

•-< 

U 

w« 

Wi 

o 

0 

0} 

e 

1M4 

w 

X 

O' 

u 

<w 

c 

Jj 

c 

< 

n 

3 

© 

<0^ 

u 

jj 

©  O 

a*. 

mmt 

3 

5 

X 

c; 

*3 

0 

< 

u 

c 

C 

2 

*c 

.* 

\ 

a; 

m 

Q} 

o 

o 

0 

© 

«»< 

o 

e 

h 

CO 

• 

(S 

•• 

5 

2 

c 

X 

b: 

X 

2 

c 

■Mr 

^( 

z 

c 

< 

o 

o 

© 

1 

•• 

2  ^ 

o 

in 

bt 

o 

C 

e 

b: 

r- 

•H 

q 

< 

^  c 

CO 

1 

«• 

W  1 

rn 

jC 

O 

c 

2  a 

1 

0> 

n 

c 

u  w 

2 

X 

fN 

CI 

H 

O  o 

u 

© 

h 

’O 

a  . 

2 

O 

o 

u 

E  •H 

2 

^o 

c 

Q) 

in 

0) 

2 

m  o  ^ 
o  XI  ♦  ^  JJ 
o  E  o  to  IS 

o  — I  t^  X 

•  X  CN  ts 

r*  w  *0  -. 

o 

•  X 


©  < 


O  m 

V  t  c 
o  E 

W  CN  3 


H  I  S 
U 
c: 

< 

© 
o 

CA 


< 

l: 


u  «  ..  .. 

1  H  -X 
2  O 
1  IN 


H 
O 

o  o 


o 

<!l 


E 

E 

2 

JJ 

3 

3 

o 

w 

-J 

•H 

o 

C 

IM 

o 

o 

m 

© 

H 

o 

CJ 

u 

O 

2 

© 

\o 

• 

r 

X 

o 

© 

m 

CO 

o 

O 

1 

z 

» 

b3 

• 

< 

n 

fS 

IN 

© 

u 

1 

VD 

1 

CO 

» 

b: 

b3 

in 

o 

2 

Z 

in 

2 

Cu 

«n 

< 

u 

O 

o 

X 

2 

CO 

Z\au*l  I 

-^<oo^  =  _  tn  O  cr^ 


c;:  a 


H  O 
o  W  M  tn 


o  cn 


z  <N  ©  u  •  O  < 

L>  •  ©  X  O  CC  Z  •-< 

<pH  c  u  u  c 

^2  £»«  O  •-*  >  — ' 

in*-iitZOr^  <-X 
X  <  O  CO  O' 

K  U  •  CJ  c:  •  *  IT. 

a  u  CO  CO  u  o  H 

>•  O  •>  £->  r  u.  tn 

H  U  2  X  X  U.  C 

r-  >«-•*-  X  e 

H  ©  <  b.  <  = 

O’  erne 

XWX22COO  - 
&2<<<  20Cu 
<2  I-  b.  Cu  fS  <n 

cn  I  •-<  2  e  X  s  «% 

«H  o  <  M  cn  CO  • 

•H  «-«  u  X  Cu  CO  tn 

CO  C  CO  t-r  tt  o 

»H  o  &3  U  ©  2  e 

ZC  t0C022  O  — 

e*  ^  CO  2  Ci2  »  r 

<  IN  o  or  X 

X  O  2  K  <  2  *-» 

—  u  IN  2  2  *-.  = 

Zb:  I  ik  o 

O  C  C  2  c  X 

o  b:  ic  2  X  < 

2  •  K  X  U  O  Z 

2  ©  M  2  b:  M  < 

©  2  >022 

^  t  X  H  b:  2 

(N  2  <  b*  b:  2  O 

2  X  b:  2  O  2  b. 

b2  CO  >  <  r 

2  2  2  X  b:  o  b. 

2HXC0<&2X 
2<C0r-*-^Ob.C 
2  2  X  Z  2  X  U 

o  ^  O  CO  X 

H  b2  o  Z  2  2^. 

OXb«X>.2CCZ 
X  -<  O  X  C  X  U 
2Ib  b«Ob.C0CfeX 

W  C  C-  C 

CO  CO  C  X  c  c 

•^*ii«22CSC02CX 

e-  u  >  2  b:  2 

C0C0M2C  222 

•-^f-rre-  rco 

ZZX2X*^C3C0 

c-t-<Z22HZ< 


x4 

n; 

Os  >0 
^  Os. 

Sj  ^ 

‘V 


] 

0 

O 

o 

{ 

\ 

K 

c 

4 

) 

op. 

S 

O 

K 

^  T  1 

■* 

\ 
'  • 

V. 

> 

0 

0 

i 

■% 

0 

n 

X 


J 


V 


4 

i\ 

Pi 


I 


< 

u 

CO 


o. 


2 

b: 

<N  B 

‘o  2 
2 
> 


CO 

2 

2 

o 

L. 


< 

> 


©CO 

5*:  DHnS  NV2S  KI  23YSSAY 


233 


o> 

a> 

O 

o 


3 

C 

\] 


Ol 

— M 

>* 

a-4 

a 

c 

03 

c 

a 

4mM 

u 

0 

t 

0 

.MJ 

a 

J3  •— 

m 

cn  ui 

o\ 

MM 

-M  fl; 

a\ 

0 

ce  Q 

q;  w 

w 

- 

CD 

0  © 

r- 

•H 

X  -4 

<m4 

*2 

cn 

•iMf  C 

.Jj 

0  -c 

n 

3 

—  Li 

Oft  o 

on—* 

3 

m2 

©  *3 

C 

X 

C  <C 

MM* 

**-*  X 

•« 

N. 

0} 

O 

© 

UJ 

O 

•3  E 

1C 

n» 

•• 

3  X 

a 

c 

o 

< 

o 

—4 

1 

•• 

C^  X) 

u 

o 

c 

•H 

0  c 

I  *3 

*• 

m2 

•4C 

o\  o 

c 

1 

31 

© 

ro 

-  « 

•H 

O 

<N 

©  H 

O 

W 

•3 

0 

o 

Li 

3  - 

to 

o 

©  oJ  in 

♦n 

o 

JLi 

o 

S3 

4mm 

Mt  t 

o 

£ 

o 

« 

a 

o 

«i4 

ro. 

m2 

< 

m2 

(N 

ui 

fN 

-Mt* 

•o 

M 

o 

3 

tn 

Jj 

O 

B 

•M 

o 

<  C 

4 

cS. 


I 


o 

W  -f 


w  • 

n  u 


I 

u 


o 

tn 


O  <  o  E  << 
to  r4  3  u 

^  C  -4  X 

'U  ■fl  ••  —  o 

I  H  mX  ^  U  H 

X  O  O  fN  c 

I  OJ  «  c  c  o  o 

0)  Wi  E  £  «H  c: 

V4  XJ  3  3  O 

«  ^  *M  ^  O 

CN  I  V-I  o  O  <n  c^ 

-H  H  O  U  U  •  O 

.Cm  O  to  • 

I  u 

H  OD  O  to 

to  o  o  in 


u 
(N  in 
in 


SI 


> 

X  to 
<  o 
u  • 

CO  u 
r* 
s  in 


in  *-1  B 
O 

a  u  • 

&  U  CO 
>  £  - 
H  U 
> 

H  o>  < 
O  • 

*2-  U  iJ 
Cm  C.  < 

<  X 

to  «  *-< 

1-1  03  t5 

»-4 

CO  X 

=  «° 
H  e- 

<  rs 
a  o 

X  u 

o  o 
o 


C  X 
U  CO 
e  X 
r  t- 

X  < 


fH  u 

o  u 
^  o 
a  u. 

u 

CO  CO  C 

M-  MM  o 

H 

CO  CO  *1 

S  X  !3 

H  H  < 


o 

I 

u 
o 

X 

-c 

-«•  O  .3 
O 

in  I  b4 
cn  X 
H  o 
u 

X  o 

-e  4a  . 
ffi  .3  O 

« 

Cm  O 

X  O 

<  o 

U  X  . 
CO  u  o 
^  H 

mm3 

Cm  < 
O  X 
X  X  CO 

<  < 

M  b« 

X  c 

•<  ca  CO 

u  X  u 
CO  M 

u  u 
CO  to  X 
CO  X  ca 
O  o 
X  ca 
U  <N  c 

I  Cm 
C  X 

ca  ca  X 

b3  C  U 

M  X  ca 
>•  o  X 
*3  H 

<  Cm  la 

z  O  c 

< 

m2  X 

CO  <  a 

t  r  Q 

M  O  CO 

X  »  X 
m2  O 
Cm  O  Cm 
O  X  ^ 

X  ca  LO 

ca  >  X 
c  o 

S  U  M 
XXX 


in 


ta 
z  z 

M  U  m 
-  U  O 
X  to 
<  b: 
CJ  X  \o 
CO  O  Cn 

CO  M  m 
to  X 
o  < 

X  X  — 4 
U  U  O 


o 


CO 

H 

X 

X  o 


H  O 

o  * 

X  o  u 

U  n  r-i 
m2  I  0\ 

to  * 
to  in 

u  c 

CO  X  I 
XU—* 


X  iC 
<  X  — 
X  M  X 
X  < 
C  m2  u 
C  <  CO 

o  X 

—  <  X 
X  M 
XX- 
X  O  r: 
^  X  X 

O  X  — 4 
X  X  c 
X  o 
U  X 

ca  H 

X  X  u 
OCX 

m2  X 

CO  X  .2 

o 

C  X 
XXX 
M  X  X 
X  X  X 

=  CO 
O  X  CO 
t*  X  < 


I 


I 


Oi 


X. 

\ 


X 

X 

r;  B 


c 

to 

X 

X 

c 

X 

CO 

X 


< 

> 


X 

> 

< 

3 


234 


c 

.2 

Vi 

o 

c; 

o 

O 


o 

0) 

c5 

> 

03 


0 

O 

C 

0 

t: 

tt  ■■Mi 
> 
0 

g 

O 

sz 

(/) 

0 
•i— < 

0 

“O 

X 
I — 
CD 

■D 

c 

0 


0 

13 

Z5 


0 


< 

X 

C/D  Q. 
^  >> 


0 

0 


0 

T3 

O 

O 

O 

-H-* 

0 

Z3 

T3 

C 

B 

'o 

0 
■  HMM 

0 

3 

3 

&i_ 

-i— ' 
0 

0 

0 

"O 

X 
{ — 
m 


0 


0 

0 

X} 

0  -C 
C  O 

0  z 

0  § 
“D 

c 
0 
J3 


C 

o 

0 

0 


0 


0 

0 


CD  O 

X  5 

CO  « 


0 


r-  — 

*nw  ^  N 


o 

0 

0 

0 

O 


0 
-i-J 

O 
c 
0 

JC 

E  m 

zL  3 
^  C 

Ti.2 

o  2 


0 

TD 

3 


0 


§;! 

■§  o 

t  - 

£-  0 

o  <15 

c  o 

0  ^ 

T  M— 

0  0 
*4—* 

^  'o)  0 


,C 

-»— '  o 
0 

0 

0  0 
>  CL 
0  CL 
^  0 


0 

0 
•  MMM 

0 

Q. 

O 

CL 

0 

-L;: 

■  MMi 

"0 

> 

0 


0 

0 


O 

3 

-♦— » 

0 

0 
a  M^M 

K 


235 


Radiative-transfer  calculations  allow  the  possibility  that  it 
is  due  to  atmospheric  gravity  waves  (GW) 

Radiance  power  spectra  consistent  with  GW  structure 


MODTRAN  SW  AND  VNIR  SPECTRAL 
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MODTRAN  FEATURES 
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10  TO  1000  TIMES  FASTER  THAN  LBL  MODELS 
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4.5  MINUTES*  (388  MINUTES  WITH  8-STREAM  DISORT) 

0.4  TO  2.5yt/m  SPECTRAL  FLUXES  AT  75  ALTITUDES 
ON  A  99  MHz  MACHINE. 


ISAACS  2-STR  :  MLS  :  60°  SUN  ;  15  CM'^  ->  10  NM 


(WN  /  g.wD  mH)  xnid  asnjdiQ  isn  +  iDauio 


1200  1600  2000  2400 

WAVELENGTH  (NM) 


m  /  g_w3  MTi)  xnij  UV10S  aavMdn 


GROUND  :  MLS  :  60°  SUN  :  15  CM"^  ->  10  NM 


NN  /  g_N3  MTl)  xnid  QHVMdn 


WAVELENGTH  (NM) 


ISAACS  2-STR  ;  MLS  ;  60°  SUN 


m  /  s-WG  MTl)  xnU  HV10S  QUVMdn 


WAVELENGTH  (NM) 


CONCLUSIONS 


LU 

Q 


< 

CC 

h- 

Q 

O 


CO  ui 
Z  CO 

o  < 


o 

Q 

LU 

CC 

0- 

O 


LU 

Q 

< 

CC 

CO 

0- 

D 

I- 


I 

I- 


H 


LU  — 
X  O 
I-  z 

^  CC 
O  DO 
I- 


O) 

LU 


LU  O 
LU  CO 


CC 

CO 

< 


O 

h— 

Q_ 

CC 

O 

CO 

00 

< 


—  LU 


LU 

CD 

I 

o 


LU 

QC 

o 


Q 

LU 

CO 

< 

LU 

-J 

LU 

CC 

LU 

00 


<  2  2 


CO 

LU 

X 


CO 


_  o 
5^0 


< 

CC 

h- 

Q 

O 


O 

O 

>■ 

CO 

DC 

< 

LU 

-J 

o 


CD 

z 
< 

DC 
LO  LU 

CO  Q 


< 

DC 

H 

Q 

O 


LU 

LU 

DC 

CO 

< 


LU  ^ 
H 

<  CO 

O  CC 

< 

o 

CO 


DC 

O 


I — 
Q 
O 


CO 

Q 


O 

UL 

o 


LU 

^  LU  o 

Q  -J  y 

LU  d.  2 

DC  X  ^ 

—  H  CO 
15  =J  LU 
O  ^  QC 

LU  ^  Q_ 

DC 


CO  O 

^  Z 


<  P 
P  O 


G 

LU 

DC 

CL 


CO 

LU 

X 


LL 


253 


ABSTRACT 


RECENT  UPGRADES  TO  MODTRAN  -  NEW  SOLAR  IRRADIANCE  AND 
INSTRUMENT  SCANNING  FUNCTIONS 

20th  Annual  Air  Force  Transmission  Meeting 
Phillips  Laboratory/GPO 
Hanscora  AFB,  MA  01731 

P.  K.  Acharya,  A.  Berk,  L.  Bernstein  and  D.  Robertson 
Spectral  Sciences,  Inc. 

99  S.  Bedford  Street,  Suite  No.  7 
Burlington,  MA  01803 

G.  Anderson  and  J.  Chetwynd 
Phillips  Laboratory/GPO 
Hanscom  AFB,  MA  01731 

H.  E.  Snell 

Atmospheric  and  Environmental  Research,  Inc. 

Cambridge ,  MA 

J.  Vail 

Mei  technologies,  Lexington,  MA 

Recent  upgrades  are  intended  to  make  MODTRAN  a  more 
accurate,  convenient  and  versatile  atmospheric  modeling  tool. 
The  current  default  solar  irradiance  database  due  to  Kurucz 
was  corrected  in  the  310-340  nm  region  using  measured  data. 
Several  other  data  sets  are  also  made  available  to  the  user. 
The  user  can  also  use  his/her  own  data  residing  in  a  file. 
Examination  confirms  that  the  corrected  Kurucz  data  is  still 
the  best  overall  among  all  currently  available  data  sets. 
MODTRAN  can  now  be  run  using  frequency  inputs  which  are  in 
wavenumber,  nm  or  micron.  Additionally,  there  is  now  the 
capability  of  obtaining  radiance  and  transmittance  which  are 
convolved  with  an  appropriate  instrument  scanning  function. 
Several  choices  of  scanning  functions  are  available  including 
an  user-defined  option.  This  upgrade  is  helpful  to  users  who 
would  like  to  model  satellite  data  which  are  gathered  as  a 
function  of  wavelength  and  have  undergone  instrument 
degradation.  All  upgrades  are  accomplished  while  maintaining 
strict  compatibility  of  TAPES  with  earlier  versions  of  the 
code.  Work  is  currently  in  progress  to  incorporate  NOVAM 
(Navy  Oceanic  Vertical  Aerosol  Models)  into  MODTRAN. 
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ATMOSPHERIC  TRANSMISSION  MODELS  (10-12  JUNE,  1997) 
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NEW  CHOICES  OF  BUILT-IN  In  FILES 
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STRUCTURE  OF  USER-SUPPLIED  In  FILE 
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PADDED  BY  DATA  FROM  DEFAULT  BUILTIN  FILE 
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ATLAS  SSBUV  IS  TAKEN  TO  BE  THE  BENCHMARK  FOR  THIS  TASK 
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(KR  =  corrected  Kurucz,  CH  =  corrected  Chance,  CB  =  SSBUV,  TH  =  SOLSPEC) 
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MAXIMUM  DIFFERENCE  IS  11%  AT  263  NM 
AGREEMENT  IS  BETTER  AT  SHORTER  WAVELENGTHS 
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KITT-PEAK  AND  CORRECTED  KURUCZ  (200-400  NM 
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KITT-PEAK  AND  CORRECTED  KURUCZ  (400-800  NM 
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KITT-PEAK  AND  CORRECTED  KURUCZ  (200-400  NM 
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INSTRUMENT  SCANNING  FUNCTION 
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IF  FWHM  IS  TOO  FINE  ERROR  MSG  IS  PRINTED  AT  THE  OUTSET 
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SOLAR  ULTRAVIOLET  IRRADIANCE  AND  ITS  VARIATION 


Giuliana  de  Toma,  Gary  J.  Rottman,  Thomas  N.  Woods 
Laboratory  for  Atmospheric  and  Space  Physics,  University  of  Colorado,  Boulder,  CO,  USA 


Abstract 

Solar  ultraviolet  irradiance,  especially  at  wavelengths  below  300  nm,  is  strongly  absorbed  by  the  Earth’s 
middle  atmosphere,  where  is  the  dominant  source  of  energy.  Knowledge  of  the  amount  of  solar  energy 
deposited  in  the  atmosphere  and  its  variation  with  time  is  important  for  atmospheric  studies.  The  Solar 
Stellar  Irradiance  Comparison  Experiment,  SOLSTICE,  on  the  UARS  satellite  makes  daily  observations 
of  solar  spectral  irradiance  in  the  wavelength  region  from  119  to  420  nm.  SOLTICE  observations,  started 
in  October  1991  and  have  continued  since  that  time.  They  cover  the  entire  declining  phase  of  solar  cycle 
22  from  solar  maximum  to  solar  minimum  providing  information  on  solar  ultraviolet  variability  over  both 
short-  and  long-time  scales.  In  this  paper,  the  SOLSTICE  irradiance  measurements  are  reviewed  with 
particular  emphasis  on  variability.  For  the  period  1991-1996,  we  report  a  variation  of  about  a  factor  of 
two  at  Lya,  a  variation  of  about  5%  near  200  nm,  gradually  decreasing  to  less  than  1%  near  300  nm. 


1  Introduction 

The  solar  irradiance  in  the  spectral  interval  120-420  nm  is  of  particular  interest  not  only  for  studying 
the  Sun  and  its  magnetic  activity  but  also  for  its  effects  on  the  terrestrial  environment  (Lean,  1991; 
Rottman  et  ah,  1996).  This  wavelength  region,  which  covers  the  far  ultraviolet  (FUV)  up  to  the  visible, 
originates  in  different  levels  of  the  solar  atmosphere  (Foukal,  1990).  The  solar  emission  at  wavelengths 
short  of  150  nm,  has  its  origin  in  the  outer  layers  of  the  solar  atmosphere:  the  chromosphere  and  transition 
region.  These  are  hotter,  less-dense  regions  where  the  magnetic  fields  of  the  Sun  play  an  important  role. 
Consequently,  radiation  emitted  from  these  layers  is  strongly  influenced  by  solar  magnetic  activity  and 
shows  significant  variability  with  the  solar  activity  cycle.  Longer  wavelengths  originate  in  lower  levels  of 
the  solar  atmosphere:  radiation  emitted  around  160  nm  comes  from  the  temperature  minimum  region 
and  even  longer  wavelengths  have  origin  in  the  photosphere.  Above  the  A1  I  ionization  edge  at  208  nm 
almost  all  the  solar  emission  is  from  the  lower  photosphere  with  the  exception  of  the  strong  absorption 
lines  of  Mg  II  at  280  nm  and  of  Ca  II  at  393  and  397  nm,  both  formed  in  the  chromosphere.  Unlike 
the  chromospheric  emission,  the  solar  photospheric  radiation  shows  much  less  variability  related  to  the 
magnetic  activity  of  the  Sun  (Foukal,  1990;  Lean,  1991). 

An  example  of  the  solar  FUV  spectrum  measured  by  SOLSTICE  is  shown  in  Figure  1.  We  note  that 
the  spectrum  below  160  nm  is  characterized  by  emission  lines  (some  of  them  not  completely  resolved  at 
SOLSTICE  resolution)  superposed  on  a  relatively  weak  continuum.  This  shorter  wavelength  region  is 
dominated  by  the  the  H  I  Lya  line  at  121.57  nm,  which  is  the  strongest  line  in  the  entire  UV  region. 
At  longer  wavelengths,  where  radiation  is  mainly  photospheric,  the  solar  spectrum  changes  appearance 
and  there  are  fewer  and  fewer  lines  in  emission.  Above  200  nm  the  spectrum  is  characterized  by  many 
absorption  lines  in  the  underlying  continuum.  The  intensity  of  the  solar  emission  is  also  a  function  of 
wavelength,  and  increases  by  several  orders  of  magnitude  from  the  H  I  Lya  to  the  visible.  The  ultraviolet 
region  is  only  a  small  fraction  of  the  total  radiative  output  of  the  Sun:  solar  irradiance  below  420  nm 
contributes  about  10%  to  the  total  solar  irradiance,  and  the  ultraviolet  region  below  300  nm,  which  is 
about  15  Watt/m^,  contributes  only  1%.  Yet,  because  the  largest  amplitude  variations  of  the  Sun’s 
output  occur  at  ultraviolet  and  radio  wavelengths  the  region  between  120  and  300  nm  contributes  to  a 
large  fraction  of  total  solar  variability. 

Solar  ultraviolet  radiation  is  very  important  because  of  its  influence  on  Earth’s  environment,  and 
solar  radiation  measurements  are  required  for  many  geophysical  studies  and  for  climate  models.  Different 
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wavelengths  affect  different  regions  of  the  atmosphere.  Solar  emission  at  wavelengths  between  300  and 
420  nm  reachs  the  ground  only  moderately  attenuated  by  the  Earth’s  atmosphere.  In  contrast,  the  spectral 
region  120-300  nm  is  almost  totally  absorbed  in  the  stratosphere  and  mesosphere  where  it  is  responsible 
for  radiative  heating  and  photodissociation  of  important  atmospheric  constituents  (O2,  O3,  H2O,  NO2, 
etc.)  (Brasseur  and  Solomon,  1986).  Changes  in  solar  ultraviolet  irradiance  can  therefore  result  in  changes 
in  the  temperature  structure  and  composition  of  the  middle  atmosphere  from  20  to  80  km  (Brasseur,  1993; 
London,  1994).  For  example,  the  ozone  distribution  is  strongly  influenced  by  the  solar  ultraviolet  radiation 
both  directly  (O3  photodissociation)  and  indirectly  (O2  and  H2O  photodissociation).  Changes  in  total 
ozone  caused  by  solar  variations  are  comparable  to  the  ozone  depletion  caused  by  the  chlorofluorocarbons 
(CFCs).  The  solar  influence  can  enhance  or  partially  cancel  this  anthropogenic  effect  (Chandra,  1991; 
Hood  and  McCormack,  1992).  Less  clear,  but  not  less  important,  is  the  effect  of  solar  variability  for  the 
wavelength  region  above  300  nm.  This  spectral  region  varies  less  than  1%  during  the  solar  cycle  and 
total  solar  irradiance  varies  only  0.1%  (Wilson,  1991).  This  small  solar  effect  seems  too  small  to  play  an 
important  role  on  the  climate  change,  yet  several  geophysical  records  seem  to  show  a  solar  signature  in 
the  past  climate  of  the  Earth  (Lean,  1994  and  references  therein). 

Measurements  of  the  solar  ultraviolet  irradiance  and  its  variability  are  thus  necessary  to  separate 
the  effects  of  solar  variability  on  the  Earth’s  middle  atmosphere  (and  possibly  the  troposphere)  from 
changes  induced  by  other  natural  phenomena  or  by  anthropogenic  activity.  This  separation  has  been 
difficult  in  the  past  because  accurate  solar  measurements  were  not  available.  Measurements  of  ultraviolet 
radiation  must  be  performed  from  space  which  makes  the  measurements  particularly  diifficult  due  to  the 
limited  lifetime  of  space  missions  and  also  to  the  degradation  of  instruments  with  time.  Until  recently, 
no  space  instrument  had  measured  the  ultraviolet  irradiance  and  its  variation  in  time  with  high  accuracy. 
The  Upper  Atmosphere  Research  Satellite  (UARS)  carries  two  ultraviolet  solar  spectrometers  and  is  the 
first  satellite  that  can  monitor  the  UV  solar  irradiance  with  high  precision  and  accuracy.  UARS  was 
launched  in  September  1991,  near  the  maximum  of  solar  cycle  22,  and  it  is  still  operating  normally.  It 
now  provides  6  years  of  reliable,  continuous  solar  irradiance  observations  from  solar  maximum  to  solar 
minimum  conditions.  In  this  paper,  the  solar  measurements  from  one  of  the  two  UARS  solar  instruments, 
the  Solar  Stellar  Irradiance  Comparison  Experiment  (SOLSTICE),  are  presented  (Section  2),  and  solar 
variability  as  a  function  of  wavelength  over  both  short-  and  long-time  scales  is  reviewed  for  the  wavelength 
range  120-300  nm  (Section  3).  Monitoring  of  the  solar  ultraviolet  radiation  will  continue  with  a  second 
generation  SOLSTICE  instrument  for  the  NASA  Earth  Observing  System  (EOS)  with  a  launch  in  the 
year  2002.  This  instrument  will  achieve  higher  accuracy  for  the  wavelength  region  300-420  nm  and  will 
extend  the  measurements  to  the  visible  and  near  infrared  regions  (Section  4). 


2  SOLSTICE  Solar  Observations 

The  SOLSTICE  instrument  on  the  UARS  satellite  is  a  3-channel  ultraviolet  spectrometer  which 
performs  full-disk  solar  irradiance  observations  in  the  spectral  region  119-420  nm  with  a  resolution  varying 
from  0.15  to  0.3  nm  (Rottman  et  al.,  1993;  Woods  et  al.,  1993,  1996).  SOLSTICE  observations  started 
in  October  1991,  about  a  month  after  UARS  launch,  and  have  continued  regularly  to  the  present  time 
with  only  a  few  gaps  (only  7%  of  the  data  are  missing).  Solar  spectra  are  measured  about  15  times  a  day 
during  the  daytime  portion  of  each  orbit,  and  reported  as  a  daily  average  with  a  1-nm  resolution.  A  list 
of  the  major  ultraviolet  emission  lines,  and  the  Mg  II  index  of  chromospheric  activity  are  also  provided 
daily.  The  data  in  irradiance  units  and  corrected  to  one  astronomical  unit  can  be  obtained  by  anonymous 
FTP  from  “ftp://daac.gsfc.nasa.gov/data/uars/solstice/”. 

The  primary  science  objective  of  SOLSTICE  is  to  measure  the  ultraviolet  irradiance  with  high  absolute 
accuracy  over  long  periods  of  time.  This  has  been  achieved  through  an  extensive  pre-flight  calibration 
(Woods  et  al.,  1993)  and  the  in-flight  capability  to  track  calibration  changes  (Rottman  et  al.,  1993).  To 
monitor  changes  in  instrumental  sensitivity,  a  group  of  bright  stars,  very  stable  at  UV  wavelengths,  are 
observed  on  a  daily  basis  during  the  nighttime  portion  of  each  orbit.  Using  these  stars  as  a  calibration  ref¬ 
erence,  it  is  then  possible  to  correct  for  instrumental  effects  and  determine  the  long-term  solar  variability 
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Figure  1:  Two  solar  FUV  spectra  from  SOLSTICE  are  shown  for  solar  maximum  and  solar  minimum 
conditions.  The  spectra  are  obtained  averaging  SOLSTICE  daily  measurements  over  three  solar  rotations, 
during  the  periods:  30  Dec.  1991  -  19  Mar.  1992  and  17  Feb.  1996  -  7  May  1996.  Several  of  the  major  solar 
lines  are  identified.  The  spectral  ranges  which  can  induce  dissociation  of  molecular  oxygen  and  ozone,  the 
two  major  absorbers  of  ultraviolet  radiation,  are  also  shown. 

with  1-2%  relative  accuracy. 

SOLSTICE  measurements  have  been  validated  by  the  observations  of  the  second  solar  instrument  on 
UARS,  the  Solar  Ultraviolet  Spectral  Irradiance  Monitor  (SUSIM),  which  covers  the  same  wavelength 
range  (Brueckner,  1993).  The  measurements  of  the  two  instruments  are  consistent  and  agree  within  their 
2-(j  calibration  uncertainty  of  about  10%  (Woods  et  al.,  1996). 

An  example  of  SOLSTICE  solar  observations  is  presented  in  Figure  1,  where  spectra  for  solar  maximum 
and  solar  minimum  conditions  are  shown.  These  spectra  are  the  average  over  81-day  (three  solar  rotations) 
during  the  winter  of  1992,  representing  solar  maximum  for  solar  cycle  22,  and  during  the  spring  of  1996, 
corresponding  to  the  recent  solar  minimum.  The  data  at  solar  minimum  are  still  preliminary  values 
because  the  correction  for  instrument  degradation  is  not  yet  final.  New  SOLSTICE  data  wdth  all  the 
corrections  included  will  be  publicly  available  in  early  1998. 
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Figure  2:  Solax  irradiance  as  a  function  of  time  in  three  wavelength  regions  important  for  the  terrestrial 
atmosphere.  Solar  variability  decreases  at  longer  wavelengths;  for  the  period  1992-1996  we  observe  a  change 
of  about  a  factor  of  two  at  the  H  I  Lya,  14%  in  the  band  180-185  nm,  and  8%  in  the  band  200-205  nm 
(notice  that  the  y-axis  for  the  center  and  bottom  plots  are  chosen  to  have  the  plots  axe  on  the  same  relative 
scale) 
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11— year  Solar  Cycle  Variability 
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Figure  3:  Solar  variability  associated  with  the  11-year  solar  activity  cycle  as  a  function  of  wavelength.  Data 
at  solar  nainimum  are  still  preliminary. 

3  Solar  UV  Variability 


In  Figure  2  the  solar  irradiance  observed  in  three  important  wavelength  regions  is  shown  as  a  function 
of  time.  The  Sun  exhibits  variability  over  several  time  scales,  from  minutes  to  years,  but  the  most  obvious 
variations  in  these  time  series  are  the  slow  change  in  the  average  intensity  related  to  the  11-year  activity 
cycle  of  the  Sun,  and  the  27-day  modulation  of  this  variation  induced  by  the  rotation  of  the  Sun  through 
the  appearance  and  disappearance  of  active  regions  on  the  solar  disk. 

In  the  top  panel  of  Figure  2  the  H  I  Lya  line  is  shown.  This  emission  line  is  very  important  because 
of  its  high  intensity  and  because  of  its  influence  on  the  mesosphere  where  it  dissociates  O2  as  well  as  CH^ 
and  H2O  with  consequent  ozone  destruction.  The  H  I  Lya  line  has  origin  in  the  upper  chromosphere 
and  is  quite  sensitive  to  solar  magnetic  activity.  We  notice  a  change  of  almost  a  factor  of  two  in  its 
intensity  from  solar  maximum  to  solar  minimum.  The  central  and  bottom  panels  of  Figure  2  show  the 
time  series  for  the  spectral  regions  180-185  and  200-205  nm  respectively.  These  emissions  correspond 
to  the  Schumann-Runge  bands  and  Herzberg  continuum  of  molecular  oxygen  and  also  influence  ozone 
densities.  In  contrast  to  H  I  Lya,  radiation  at  these  wavelengths  originates  in  lower  and  cooler  regions  of 
the  solar  atmosphere.  Variability  is  still  significant,  but  quite  smaller  than  for  the  H  I  Lya. 

3.1  Solar  Cycle  Variability 

Solar  cycle  variability  is  associated  wdth  the  evolution  of  the  solar  magnetic  fields  during  a  period 
of  approximately  11  years  from  solar  maximum  to  solar  minimum.  At  solar  maximum  several  active 
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Figure  4:  Solar  variability  associated  with  the  27-day  rotational  modulation  as  a  function  of  wavelength  for 
winter  1992,  during  solar  maximum  activity. 


regions,  which  are  regions  of  emergence  of  magnetic  fluxes,  are  present  on  the  surface  of  the  Sun,  while  at 
solar  minimum  the  solar  surface  appears  relatively  uniform  and  “quite”.  At  ultraviolet  wavelengths,  dark 
sunspots  are  not  very  important  and  radiation  is  enhanced  at  solar  maximum  by  the  presence  of  bright 
active  regions  (plages)  surrounding  the  sunspots. 

To  estimate  solar  variability  over  the  11-year  solar  cycle,  we  have  used  the  average  spectra  for  the 
winter  of  1992  and  spring  of  1996  shown  in  Figure  1.  Their  ratio  in  1-nm  bin  intervals  is  plotted  in  Figure 
3.  Solar  maximum  was  reached  in  1989-1991,  but  the  winter  of  1992  was  a  period  of  very  high  solar  activity 
and  well  represents  solar  maximum  conditions.  The  data  in  spring  of  1996  correspond  to  the  recent  solar 
minimum.  These  latter  SOLSTICE  measurements  are  still  preliminary  so  the  results  shown  here  must  be 
considered  only  a  first  estimate  for  variability  during  solar  cycle  22.  Nevertheless,  we  can  call  attention  to 
several  interesting  aspects  of  solar  variability  seen  in  Figure  3.  Solar  variability  is  considerably  larger  at 
the  shorter  wavelengths  (especially  in  correspondence  of  the  emission  lines)  and  decreases  at  the  longer 
wavelengths  as  we  move  from  hotter  to  cooler  emission  regions  of  the  solar  atmosphere.  The  H  I  Lya  varies 
70%  or  more,  the  wavelength  region  around  170  nm  varies  about  10%,  and  the  region  at  200  nm  about 
5-6%.  Then  variability  decreases  sharply  above  the  A1  I  ionization  edge  at  208  nm,  with  the  exception 
of  the  Mg  II  emission  at  280  nm,  which  is  of  chromospheric  origin.  For  wavelengths  longer  than  300  nm 
the  solar  variability  is  less  than  1%  and  is  beyond  the  accuracy  limit  of  the  SOLSTICE  instrument. 

3.2  Rotational  Modulation 

The  amplitude  of  the  27-day  modulation  caused  by  the  Sun’s  rotation  changes  throughout  the  solar 
cycle.  If  we  examine  Figure  2  we  notice  episodes  of  activity  which  indicate  the  emergence  of  new  bright 
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active  regions.  These  regions  dominate  for  several  rotations  (typically  4  to  8  rotations)  and  then  wane.  At 
solar  maximum,  when  the  Sun  is  very  active,  the  amplitude  of  the  27-day  modulation  is  generally  larger, 
and  then  gradually  decreases  with  decreasing  solar  activity.  This  change  from  maximum  to  minimum  is 
not  strictly  monotonic,  as  the  appearance  of  new  active  complexes  can  cause  a  new  pulse  of  activity;  for 
example  in  early  1994.  At  solar  minimum,  when  the  Sun  is  very  quite,  we  have  times  when  the  rotational 
modulation  is  almost  absent;  an  example  is  the  spring  of  1996,  when  it  is  very  difficult  to  identify  the 
rotational  signal  in  the  data. 

To  study  the  27-day  modulation  as  a  function  of  wavelength,  we  have  chosen  a  solar  rotation  near 
maximum  solar  activity,  when  rotational  modulation  is  quite  large.  We  ratio  two  SOLSTICE  solar  spectra 
taken  on  February  27  and  March  6,  1992.  These  days  are  identified  by  asterisks  in  Figure  2  and  correspond 
to  the  peak  and  the  valley  of  a  solar  rotation  period.  The  ratio  of  peak  to  valley  in  1-nm  bin  intervals 
are  shown  in  Figure  4.  This  plot  has  many  similarities  with  Figure  3.  Variations  are  larger  at  the  shorter 
wavelengths  and  decrease  at  the  longer  wavelengths  which  originate  in  the  photosphere.  Again  we  notice 
the  enhancement  in  variability  of  the  chromospheric  Mg  II  line  and  the  drop  associated  with  the  A1  I 
edge.  As  we  move  toward  visible  wavelengths  (not  shown  here)  the  effect  of  sunspots  starts  to  appear, 
and  we  sometimes  observe  a  slight  decrease  in  intensity  associated  with  their  passage  across  the  solar  disk 
instead  of  an  increase  associated  with  plages. 

4  Conclusions  and  Future  Missions 

Variations  in  the  ultraviolet  irradiance  can  significantly  affect  the  energy  balance,  photochemistry  and 
dynamics  of  the  middle  atmosphere  between  20  and  80  km.  Therefore,  knowledge  of  the  magnitude  and 
variability  of  solar  output  is  important  to  understand  the  solar  influence  on  the  Earth’s  atmosphere,  and 
to  separate  solar  effects  from  other  natural  and  anthropogenic  forcing. 

SOLSTICE  measurements  from  the  UARS  satellite  now  provide  a  6-year  set  of  ultraviolet  irradiance 
data  with  better  accuracy  and  precision  than  previous  measurements.  The  careful  pre-flight  calibration 
of  the  SOLSTICE  instrument  and  the  validation  of  SOLSTICE  measurements  with  SUSIM  measurements 
have  allowed  us  to  establish  the  value  of  solar  ultraviolet  irradiance  with  higher  photometric  accuracy, 
especially  at  the  shorter  wavelengths.  Moreover,  the  length  of  the  UARS  mission  which  now  covers 
the  entire  declining  phase  of  solar  cycle  22,  provides  new  information  about  solar  ultraviolet  variability. 
Because  of  the  SOLSTICE  in-flight  calibration  technique  it  is  now  possible  to  establish  solar  UV  variability 
over  long-time  scales  with  an  accuracy  of  1-2%,  and  with  an  even  higher  accuracy  over  short-term  scales 
of  the  order  of  a  solar  rotation. 

This  1-2%  accuracy  limit  is  appropriate  for  the  shorter  ultraviolet  wavelengths  where  solar  emission  is 
highly  variable,  but  can  only  set  an  upper  limit  for  the  photospheric  emission  at  wavelengths  longer  than 
300  nm,  where  variability  is  typically  less  than  1%.  In  the  visible  and  infrared  regions,  solar  variability 
is  known  to  be  small,  but  very  little  is  known  about  the  spectral  distribution  of  variations  since  only  the 
total  solar  irradiance  variation  has  been  recorded.  Because  this  portion  of  the  solar  radiation  reaches 
the  Earth’s  surface  and  affects  directly  the  troposphere,  its  long-term  variation  can  influence  the  Earth’s 
climate.  Accurate  measurements  are  therefore  important  to  establish  the  solar  contribution  to  climate 
change.  At  this  purpose,  a  new  improved  version  of  SOLSTICE,  SOLSTICE  II,  will  be  flown  in  the 
year  2002  as  part  of  the  EOS  Mission.  SOLSTICE  II  will  continue  to  monitor  the  ultraviolet  region 
between  120-300  nm  using  a  grating  spectrometer  similar  to  the  current  SOLSTICE.  In  addition,  an 
electrical  substitution  radiometer  and  prism  spectrometer  will  cover  the  visible  and  near  infrared  region 
up  to  2.0/im  with  a  relative  accuracy  of  0.01%.  New  solar  measurements  will  also  be  available  in  the  near 
future  at  the  shorter  ultraviolet  wavelengths.  The  Solar  EUV  Experiment  (SEE),  on  the  TIMED  mission 
planned  for  the  year  2000,  will  carry  two  solar  instruments  which  will  make  daily  solar  measurements 
from  0.1  to  200  nm.  They  will  provide  new  information  on  the  spectral  region  very  important  for  the 
upper  atmosphere  above  80  km  whose  long-term  variability  is  still  poorly  known. 
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Abstract 

The  Rotating  Shadowband  Spectroradiometer  (RSS)  is  a  new  instrument  that  provides  spectral 
direct,  diffuse,  and  total  horizontal  irradiances  with  a  512  CCD-array  over  the  wavelength  range 
from  350  nm  to  1075  nm.  We  compare  the  measured  spectra  from  the  RSS  to  that  taken  by  a 
collocated  MFRSR,  and  MODTRAN-3.5  calculations.  The  measurements  from  both  instruments 
(RSS  and  MFRSR)  agree  within  2%.  The  model  calculated  direct-horizontal  irradiances  show  very 
good  agreement  with  the  measurements.  However,  the  model  substantially  overpredicts  the  diffuse- 
horizontal  if  the  standard  “rural  aerosol”  model  is  applied.  If  instead  we  use  MODTRAN’s  “urban 
aerosol”  properties  for  the  asymmetric  factor  and  single  scattering  albedo,  the  discrepancy  can  be 
eliminated.  The  spectral  signatures  of  discrepancy  suggest  that  the  scattering  treatment  of  the 
diffuse  irradiance  and  uncertainties  about  aerosol  microphysics  may  be  the  reasons  for  “clear-sky 
anomaly” . 


Introduction 

The  two  prototype  RSS  instruments  were  operated  at  the  H2O  Intensive  Operation  Period 
held  at  the  Southern  Great  Plains  (SGP)  field  site  of  the  Atmospheric  Radiation  Measurement 
(ARM)  Program,  September  1996.  Both  instruments  operated  continuously  without  failure  for 
the  campaign.  Each  instrument  was  calibrated  at  least  once  per  day  to  assess  instrument  stability 
under  field  conditions.  The  instruments  were  stable  and  reproducible  to  approximately  0.3%  in 
irradiance  (the  limit  of  stability  of  our  calibrating  standard)  over  the  experiment  period,  with 
no  statistically  significant  trend.  For  the  RSS,  the  external  dimensions  of  the  fore-optic  diffuser 
and  associated  shadowband  mechanism  are  identical  in  scale  and  operation  to  that  of  the  Multi¬ 
filter  Rotating  Shadowband  Radiometer  (MFRSR)  [Harrison  et  al.,  1995].  However  the  internal 
geometry  of  the  fore-optic  is  substantially  different,  and  instead  feeds  a  prism  spectrograph  with 
a  512  element  linear  CCD-array  detector  [Harrison  et  al.,  1996]. 

We  will  compare  the  clear-sky  data  to  that  taken  by  a  collocated  MFRSR  and  model  calculations 
using  MODTRAN-3.5  [Berk  et  al.,  1989]  during  H2O  lOP  at  SGP  site.  After  the  lOP  campaign 
with  a  new  windowless  CCD  array  device,  the  instrument’s  out-of-band  rejection  has  been  im¬ 
proved.  We  further  compare  the  improved  measurements  observed  at  ASRC  with  the  MODTRAN 
calculations. 
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Comparison  with  MODTRAN  and  MFRSR 


The  modeled  spectra  were  produced  using  MODTRAN  3,5:  a  recent  release  that  includes  a 
discrete-ordinates  scattering  code  (DISORT)  [Stamnes,  1988],  so  that  scattering  atmospheres  can 
be  treated.  All  calculations  were  done  using  16  ordinates.  Spectra  were  computed  at  much  higher 
spectral  resolution,  and  then  convolved  with  an  estimate  of  the  RSS  slit-function  to  produce  the 
results  shown  here.  The  standard  fall  mid-latitude  H2O  column  was  used  without  modification;  in 
this  particular  case,  September  22,  1996,  shown  in  Figure  (1),  it  was  close  to  the  real  atmospheric 
state.  Calculations  were  done  for  two  aerosol  models:  MODTRAN ’s  standard  “urban”  and  “rural.” 
In  both  cases  the  visibility  was  set  to  54  km:  a  value  that  produced  best  agreement  in  the  direct- 
horizontal.  (Equivalent  to  using  the  direct-beam  measurement  to  determine  the  optical  depth.) 


at  SGP,  together  with  data  from  a  collocated  MFRSR  and  results  modeled  using  MOD¬ 
TRAN  3.5. 


In  the  top  panel  of  Figure  (1)  (the  direct-horizontal)  we  can  see  that  agreement  is  exceptional, 
for  all  wavelengths  less  than  900  nm.  At  wavelengths  greater  than  900  nm  the  RSS  shows  an 
irradiance  greater  than  that  seen  by  the  MFRSR  at  940  nm,  or  than  MODTRAN  predicts.  For 
wavelengths  beyond  1  fim  the  slope  is  aphysical.  The  RSS  data  from  SGP  are  affected  at  these 
longest  wavelengths  by  inadequate  out-of-band  rejection.  This  problem  is  due  to  reflections  asso¬ 
ciated  with  the  window  in  front  of  CCD  array  (Modifications  made  to  the  instrument  after  the 
field  campaign  have  eliminated  this  artifact,  and  arguably  we  can  recover  accurate  data  from  the 
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Figure  2  RSS  data  for  a  single  observation  near  local  solar  noon  on  Nov.  4,  1996  at 
ASRC,  and  results  modeled  using  MODTRAN  3.5. 


campaign  via  deconvolution).  However,  for  the  moment  we  should  simply  ignore  RSS  data  taken 
during  this  campaign  at  wavelengths  >  900  nm. 

The  second  two  panels  show  that  if  the  standard  “rural  aerosol”  model  is  applied  it  then  over¬ 
predicts  the  total-horizontal  shortwave  irradiance  by  roughly  8%.  This  result  is  reminiscent  of 
current  discussions  about  a  “shortwave  anomaly”  based  on  broadband  shortwave  data  and  models 
[Wild  et  al.,  1995;  Kato  et  al.,  1996].  They  see  comparable  agreement  of  direct-beam  irradiance 
between  measurements  (active  cavity  radiometer)  and  models,  and  model  overprediction  of  diffuse 
broadband  irradiance  compared  to  pyranometer  measurements.  The  situation  is  diflficult  to  un¬ 
tangle  solely  with  broadband  data  both  because  without  spectral  signatures  it  has  been  difficult 
to  ascribe  the  mechanism,  and  calibration  uncertainties  associated  with  the  pyranometers  (and 
hence  the  accuracy  of  derived  direct/diffuse  ratios)  are  difficult  to  rule  out. 

The  third  panel  demonstrates  that  the  discrepancy  seen  between  the  prediction  assuming  a  rural 
aerosol  and  the  observations  clearly  has  an  aerosol-like  signature.  Further,  the  discrepancy  can 
be  eliminated  by  choosing  the  “urban  aerosol”  model  that  decreases  the  assumed  single-scattering 
albedo  and  asymmetry  factor.  At  the  moment  the  microphysical  optical  properties  of  the  aerosol 
over  the  SGP  site  are  not  observed,  and  so  the  question  of  whether  the  assumption  of  an  “urban 
aerosol”  is  plausible  remains  moot.  The  issue  awaits  further  collection  of  data  (where  we  may  see 
a  range  of  aerosol  properties),  and  hopefully  aerosol  optical  measurements  aloft. 
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Figure  (2)  shows  the  comparison  of  MODTRAN  calculation  with  an  terrestrial  spectrum  taken 
with  RSS  at  ASRC  (in  Albany  NY)  after  the  windowless  CCD  array  device  was  substituted.  The 
long  wavelength  data  are  now  accurate. 

In  this  case  the  data  were  taken  at  ASRC  on  Nov  4, 1995,  and  again  the  default  standard  seasonal 
atmosphere  (for  H2O)  and  urban  aerosol  model  were  assumed.  In  this  case  the  atmosphere  clearly 
had  a  greater  water  column  than  model,  but  otherwise  agreement  is  excellent  throughout  the 
spectrum,  except  for  the  A-band  of  02  .  The  likely  cause  of  this  discrepancy  is  a  mild  underestimate 
of  the  width  of  the  instrument-function  used  to  convolve  the  spectrum,  or  the  overestimate  of  02 
absorption  coefficients  in  MODTRAN3.5. 


Conclusion 

We  compare  the  data  to  that  taken  by  a  collocated  MFRSR  and  model  calculations  using 
MODTRAN-3.5.  Agreement  is  generally  excellent  (better  than  2%,  an  optimistic  assessment  of 
our  calibrator’s  absolute  uncertainty).  The  modeled  diffuse  irradiance  is  substantially  too  large 
for  the  two  clear  days  of  the  lOP  in  comparison  to  the  observed  results  unless  MODTRAN’s 
“urban  aerosol”  properties  for  the  single-scattering  albedo  and  phase  function  are  used.  If  the 
“rural  aerosol  model”  is  used  then  the  total  surface  irradiance  predicted  by  the  model  is  about 
8%  too  large.  The  broad  spectral  behavior  of  this  potential  discrepancy,  and  it’s  presence  only  in 
the  diffuse  component,  clearly  identifies  it  as  associated  with  scattering  processes,  not  trace  gas 
absorption.  The  spectral  signatures  of  discrepancy  suggest  that  the  scattering  treatment  of  the 
diffuse  irradiance  and  uncertainties  about  aerosol  microphysics  may  be  the  reasons  for  “clear-sky 
anomaly” .  In  the  absence  of  more  data,  and  hopefully  in-situ  measurements  of  single-scattering 
albedo  and  asymmetry  parameter,  the  question  of  whether  there  is  a  “clear  sky  anomaly”  remains 
moot. 
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Architecture  of  the  MWIR  Data  Correction 
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Data  Processing  and  Analysis  Procedure 
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GOES-8  MWIR  Imagery  (with  Buoy  Locations) 
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SURFACE  TEMPERATURE  RETRIEVAL 
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MWIR  Atmospheric  Correction  Algorithm 
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ABSTRACT 

The  Performance  and  Analysis  Capability  for  Earth  Observation  Systems  (PACEOS™)  computer  simulation  software  was  devel¬ 
oped  as  an  engineering  tool  by  Lockheed  Martin  to  assist  in  the  evaluation  of  concepts  and  designs  for  space-based  remote  sensing 
systems.  This  development  took  place  over  the  past  ten  years.  PACEOS™  models  the  entire  remote  sensing  process  from  the  radi¬ 
ation  source  to  the  information  extracted  by  the  user.  The  simulation  is  divided  into  10  primary  sections:  Input  Scene,  Scene  Genera¬ 
tion,  Sensor  Performance,  On-board  Processing,  Data  Compression,  Downlink,  De-Compression,  CJround  Processing,  Metrics,  and 
Analysis,  PA(3EOS™  simulates  the  performance  of  earth  and  atmospheric  remote  sensing,  surveillance,  and  early  warning  systems 
and  is  primarily  used  for  concept  definition  and  evaluation  of  hardware/software  performance  to  evaluate  the  effects  of  sensor,  plat¬ 
form,  and  processing  parameters  on  the  quality  and  accuracy  of  the  information  provided  to  the  ultimate  users. 

This  paper  discusses  the  uses  of  the  MODTRAN  software  in  the  Scene  Generation  module  of  PACEOS™.  The  MODTRAN  software 
is  used  to  calculate  radiances  and  transmissions  that  are  applied  to  the  input  scene  as  needed  by  the  PACEOS™  software  for  electro- 
optical  systems.  PACEOS™  sets  up  and  runs  MODTRAN  and  extracts  the  parameters  needed  to  generate  the  scene.  The  PACEOS™ 
user  has  the  option  to  calculating  radiances  as  a  function  of  albedo,  emissivity,  or  temperature  from  an  input  image  or  from  an  image 
comprised  of  material  characteristics.  Options  to  vary  atmospheric  profiles,  sensor  zenith  angle,  and  sun/moon  position  are  avail¬ 
able,  as  well  as  the  ability  to  insert  cloud  layers  over  the  surface  scene.  RADTRAN  is  used  in  a  similar  way  for  microwave  scene 
generation. 

Keywords:  PACEOS™,  scene  generation,  cloud  layering,  material  mode,  non-homogeneous  atmospheres,  cloud,  bandpass  modifi¬ 
cation 


1.  OVERALL  PACEOS™  PROCESS 

PACEOS™  is  used  to  assess  the  performance  of  current  and  proposed  electro-optical  (EO)  sensor  systems  by  modeling  the  effects 
of:  1)  scene  phenomenology  and  propagation  to  the  sensor,  2)  sensors,  detectors,  and  platforms,  3)  on-board  processing,  4)  down¬ 
link,  5)  ground  processing  and  analysis,  and  6)  end  user  product  generation.  The  overall  PACEOS™  functions  are  illustrated  in 
Figure  1.1.  PACEOS™  is  primarily  used  to  generate  data  from  which  hardware  and  software  trades  can  be  made  based  on  their  effect 
on  user  information  requirements.  Interfaces  to  various  data  bases  and  external  codes  are  provided  to  assist  in  the  assessment  of 
the  performance  of  current  and  proposed  UV,  visible,  IR,  and  microwave  systems  for  surface  and  atmospheric  remote  sensing,  early 
warning,  and  surveillance. 

The  PACEOS™  software  can  be  run  interactively  or  in  various  batch  modes.  PACEOS™  is  designed  to  be  user-friendly  in  the  inter¬ 
active  mode  with  a  large  variety  of  pull-down  and  pop-up  menus  on  a  Sun  or  DEC  Alpha  workstation.  The  user  is  first  presented 
with  the  overall  menu  shown  in  Figure  1.2,  from  which  any  portion  or  an  entire  set  of  the  software  can  be  run.  With  each  primary 
section  chosen,  the  user  can  select  the  desired  level  of  simulation  detail,  which  determines  the  options  available  in  the  remaining 
menus.  If  the  user  is  less  familiar  with  a  given  portion  of  the  simulation,  a  level  with  less  detail  can  be  chosen  and  defaults  used 
for  some  of  the  less  critical  choices.  If  the  user  is  familiar  with  a  given  portion,  a  more  detailed  level  can  be  chosen  with  more  control 
over  the  simulation  parameters.  A  discussion  of  the  overall  sunulation  was  presented  m  a  1994  SPIE  paper,  “Geometric  Mapping 
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Figure  1.1  Overall  PACEOS  Functional  Areas 

in  an  Advanced  Remote  Sensing  Simulator”  ^  Other  details  on  PACEOS™  scene  generation  were  discussed  in  a  1995  SPIE  paper. 
“Input  Scene  Generation  for  System  Performance  Evaluation”^. 


2.  INPUT  SCENE  SELECTION 

The  PACEOS™  process  begins  with  the  selection  of  an  input  scene  from  data  bases  COTtainmg  different  scene  categories.  A  large 
variety  of  satellite  and  airborne  single  or  multi-band  scenes  are  available  to  choose  from,  as  well  as  scenes  made  up  of  materials 
from  either  a  materials  classification  code.  Environmental  Data  Record  (EDR)  algorithms  or  manually  generated  scenes.  Various 
routines  are  available  to  apply  geometric  corrections  to  the  scene.  Other  routines  allow  the  user  to  modify  the  intensity  of  pixels 
within  the  scene  to  represent  different  albedos  or  emissivities  and  to  add  targets  and  test  patterns  to  the  scenes.  Other  input  scenes 
have  be  generated  for  space  or  limb  viewing  and  complete  simulations  have  been  made  for  airborne  as  well  as  space-based  systems. 

MODTRAN  software  is  used  to  calculate  the  emitted  and  scattered  radiance  for  scenes  prq)agated  through  an  atmosphere.  PA¬ 
CEOS™  sets  up  most  of  the  needed  MODTRAN  parameters  based  cm  the  conditions  to  be  simulated  and  gives  the  user  the  opportuni¬ 
ty  to  view  and  modify  the  inputs  if  desired.  No  changes  internal  to  MODTRAN  have  been  made;  however,  a  large  amount  of  software 
written  in  C,  UNIX  and  PVwave  has  been  generated  to  setup  the  input  files  and  to  extract  the  calculated  temperatures  and  wave¬ 
lengths.  MODFIL  is  used  to  apply  spectral  filters  to  the  MODTRiW-calculated  radiances  and  transmissions.  The  remainder  of 
this  paper  discusses  the  setup  and  running  of  MODTRAN  and  the  use  of  its  outputs  in  generating  input  scenes.  This  methodology 
is  described  by  showing  two  examples,  one  for  a  material-based  earth  surface  input  scene  and  one  for  an  earth  limb-viewing  scene 
from  space. 


326 


Figure  12  Main  PACEOS  Control  Window 


3.  MATERIAL  BASED  INPUT  SCENE  GENERATION 

In  this  example,  the  radiance  distribution  of  a  5 12x5 12  pixel  input  scene  generated  using  tibree  surface  materials  with  varying  surface 
temperature  will  be  propagated  through  a  non-homogeneous  atmosphere  with  low  and  high  altitude  cloud  layers  using  the  spectral 
response  of  AVHRR  band  4  centered  at  approximately  10.8  [im.  The  sensor  zenith  angle  and  the  solar  azimuth  and  elevation  angles 
will  also  be  varied  across  the  scene.  The  initial  input  scene  used  for  this  run  is  shown  in  Figure  3.1. 

3.1  BAND  FILE  SELECTION  AND  CONTENT 

Most  of  the  sensor  and  platform  related  parameters  used  by  PACEOS™  are  contained  in  ‘‘band”  files.  These  files  are  generated  prior 
to  or  during  the  running  of  PACEOS™,  a  separate  file  is  generated  for  each  band  of  each  sensor.  For  the  simulation  results  shown 
in  this  paper,  band  4  of  the  AVHRR  sensor  was  diosen.  Parameters  in  this  file  include  values  that  define  the  beginning  and  ftnHing 
wavenumbers  and  the  resolution  of  the  MODTRAN  run,  the  viewing  geometry,  and  the  system’s  spectral  response  as  a  function 
of  wavelength  for  MODFBL.  These  values  are  used  along  with  default  values  to  generate  the  interactive  input  window  shown  in 
Figure  3.2. 

Any  of  the  values  presented  in  this  window  can  be  changed  by  the  user  prior  to  or  during  the  simulation.  In  this  case,  MODFIL  will 
be  run;  therefore,  the  bandpass  specified  for  the  MODTRAN  run  is  wide  enough  to  accommodate  all  the  wavelengths  defined  in 
the  MODFIL  input  fiOie.  All  of  the  MODTRAN  options  for  atmospheric,  haze,  and  dust  models  can  be  selected  from  this  window. 
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0  1  2 


Figure  3.1.  Initial  Material  Input  Scene 

ff  a  user-defined  atmospheric  model  is  chosen,  PACEOS™  will  present  other  windows  through  which  the  user  can  build  the  desired 
profile.  Oioosing  to  vary  temperature  across  the  scene  results  in  die  amplitude  of  the  scene  being  scaled  between  two  temperature 
values  (which  arc  defined  later  in  the  run)  and  the  setup  of  the  MODTRAN  input  file  to  run  cases  at  the  two  temperature  values. 
The  option  to  force  MODTRAN  to  run  was  installed  to  allow  the  user  to  skip  the  running  of  MODTRAN  if  a  identical  run  was  done 
previously  and  the  output  files  saved.  The  parameters  illustrated  in  Figure  3.2  result  in  MODTRAN  being  run  at  1  cm“i  resolution 
from  833  to  1030  cm"  i. 


3.2  MATERIAL  MODE  OPERATION 

When  an  input  scene  consisting  of  several  materials  is  chosen,  PACEOS™  sets  up  to  run  MODTRAN  sequentially  for  each  matftrial 
The  individual  material  output  scenes  are  merged  into  a  single  complete  image.  In  this  mode  of  operatirai,  it  is  necessary  to  have 
each  material  in  the  scene  identified  with  the  appropriate  index  value  for  that  material  in  the  MODTRAN  reference  background 
(refbkg)  file.  It  is  also  necessary  for  the  refbkg  file  to  have  reflectance  values  for  all  the  materials  in  the  scene  across  the  bandpass 
being  simulated. 

The  selected  input  scene,  shown  in  Figure  3.1,  is  comprised  of  three  materials:  ocean,  desert,  and  snow.  The  refbkg  file  was  extended 
to  include  values  for  these  materials  at  high  spectral  resolution  across  the  band  being  simulated.  10.8  pm  for  this  run,  and  the  materi¬ 
als  were  assigned  numerical  values  that  correspond  to  the  material  index  values  in  the  modified  refbkg  file. 

3.3  NON-HOMOGENEOUS  ATMOSPHERE 

The  non-hompgeneous  atmosphere  mode  is  chosen  in  PACEOS™  to  allow  variations  in  the  atmospheric  profile,  viewing  angles, 
and  sun/moon  position  across  the  scene.  Nine  points  are  selected  for  profile  specification;  one  in  each  comer  of  the  scene,  one  in 
the  center  of  the  scene,  and  one  midway  down  and  across  the  edges  of  the  scene  as  illustrated  by  the  numbers  overlaying  the  scene 
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Figure  32  Fully  Editable  MODTRAN  Parameters  Window 

in  Figure  3.1.  The  surface  temperature,  tbound,  was  varied  for  each  material  as  shown  in  Table  3.1.  The  values  that  were  varied 
for  each  of  the  grid  pdnts  at  which  the  conditions  are  specified  ate  shown  in  Table  3.2. 

Table  3.1 .  Specified  Surface  Temperature  (IQ  for  Materials 


Material 

Point  0 

Point  1 

Point  2 

PoirrtS 

Point  4 

Points 

Points 

Point? 

Points 

Ocean 

280 

280 

278 

274 

273 

272 

271 

271 

270 

Snow 

270 

270 

270 

271 

270 

270 

270 

269 

268 

Desert 

279 

279 

279 

273 

273 

273 

270 

269 

270 
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Table  3.2.  Parameters  Varied  at  Atmospheric  Profiie  Grid  Points 


Grid 

Point 

Atmospheric  Model 

Haze  Model 

Dust  Model 

Sun/Moon 

El,  Az  (degrees) 

Sensor  Zenith 
Angle  (degrees) 

0 

6  -  US  Standard 

4  -  Maritime 

8  -  Extreme 

1/2  Moon  30, 180 

180 

1 

6  -  US  Standard 

4  -  Maritime 

3  -  High/Fresh 

Sun 

70,180 

160 

2 

6  -  US  Standard 

2  -  Rural,  5  km  vis 

3  -  High/Fresh 

Sun 

90,180 

140 

3  -  Midlat  Winter 

2  -  Rural,  5  km  vis 

5  -  Moderate 

1/2  Moon  10, 180 

180 

m 

3  -  Midlat  Winter 

1  -  Rural,  23  km  vis 

5  -  Moderate 

Sun 

60,180 

160 

3  -  Midlat  Winter 

1  -  Rural,  23  km  vis 

5  -  Moderate 

Sun 

80,180 

140 

5  -  Subarctic  Winter 

None 

None 

Sun 

0, 180 

180 

■ 

5  -  Subarctic  Winter 

None 

None 

Sun 

20, 180 

160 

[s 

5  -  Subarctic  Winter 

None 

None 

Sun 

40,180 

140 

3.4  SCENE  RADIANCE  CALCULATES 

No  cloud  layers  were  added  during  this  run.  The  resultant  radiance  distributiou  is  shown  in  Figure  3.3.  The  effects  of  the  surface 
temperatures  and  atmospheric  models  used  can  be  seen  verticaUy  across  the  scene  with  higher  radiances  toward  the  top.  The  effect 
of  the  haze  can  be  seen  by  comparing  tire  lower  left  comer  to  the  center  right  region  with  the  snow  and  desert  matp.riak  In  the  lower 
left  comer  with  no  haze,  it  is  easier  to  see  the  differences  in  radiances  than  in  the  center  right  where  the  rural  haze  with  23  km  visibili¬ 
ty  was  added. 

Generation  of  this  image  required  many  individual  MODIUAN  runs.  There  is  no  limit,  other  than  the  time  required  to  setup  and 
run  the  simulation,  and  the  number  of  atmospheric  grid  points  that  can  be  specified. 

4.0  CLOUD  LAYERING  OPTION 

An  input  scene  with  a  low  and  high  altitude  cloud  layer  was  simulated  for  band  1  of  AVHRR  centered  at  0.87  [tm.  This  simulation 
was  run  with  the  same  non-homogenous  atmospheric  conditions  as  the  10.8  [un  scene.  A  stratus  cloud  layer  and  a  0.15  km  cirrus 
cloud  layer  at  a  base  height  of  12  km  were  added  to  the  scene.  The  spatial  distribution  of  each  cloud  layer  was  specified  by  previously 
generated  maps  which  were  scaled  by  PACEOS™  to  match  the  input  scene  size. 

4.1  VARIATION  OF  CLOUD  RADIANCES  AND  TRANSMISSIONS 

Various  algorithms  were  developed  for  PACEOS™  to  generate  realistic  cloud  radiances  and  transmission  distributions  and  to  overlay 
them  on  the  surface  scene.  The  process  involves  generating  a  map  using  fractal  routines  that  represents  the  type  of  spatial  distributicm 
desired.  Several  of  these  maps  have  been  generated  and  are  accessed  from  PACEOS™  for  both  low  and  high  altitude  cloud  layers. 
MODTRAN  is  setup  to  generate  radiances  and  transmissicais  with  and  without  the  clouds  present  and  interpolates  between  the  val¬ 
ues  based  on  the  weighting  obtained  from  the  cloud  distribution  maps  specified  for  each  run. 

Calculations  are  made  for  various  altitude  layers  in  the  atmosphere.  The  first  layer  calculated  is  from  the  surface  scene  to  the  top 
of  the  first  cloud  layer.  In  this  case,  the  surface  was  specified  to  an  altitude  of  0  km  for  all  atmospheric  grid  points  and  the  top  of 
the  lower  altitude  cloud  layer,  the  stratus  clouds,  was  specified  to  an  altitude  of  1  km.  Since  this  run  was  done  in  the  material  mode, 
the  radiance  distribution  for  each  material  was  calculated  separately,  and  the  sections  were  merged  together  to  form  the  radiance 
distribution  at  the  top  of  the  cloud  layer.  The  results  of  these  calculations  are  illustrated  in  Figure  4.1.  Many  additional  MODTRAN 
and  MODFIL  runs  were  made  to  generate  the  radiance  distribution  at  the  top  of  this  altitude  layer. 

The  cloud  distribution  chosen  for  this  layer  has  large  secticms  with  no  clouds  present.  In  these  legicms  the  surface  can  be  seen.  The 
rffect  of  the  moon  illumination  can  be  seen  in  the  very  low  minimum  radiance  and  in  the  upper  left  side  of  the  image  where  the 
clouds  over  the  ocean  are  very  dark. 
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F^ure  33.  Calculated  10.8  |un  Radiance  Distribution  at  Sensor  Aperture 
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The  second  layer  calculated  is  from  the  tqp  of  the  stratus  cloud  layer  to  the  top  of  the  cirrus  layer.  In  this  case,  separate  runs  for 
each  material  are  not  needed  since  the  radiance  at  the  top  of  the  lowest  layer  is  calculated  with  all  the  effects  of  the  materials  ac¬ 
counted  for.  The  radiance  distribution  at  the  top  of  this  layer  is  not  shown  separately.  Instead  the  radiance  distribution  at  the  sensor 
aperture  is  presented  in  Figure  4.2.  This  distribution  includes  the  remainder  of  the  atmospheric  path  radiance  and  transmission  ef¬ 
fects. 


Figure  4,2,  Radiance  Distribution  for  0,87  Scene  at  Sensor  Aperture 


Due  to  its  thinness  the  cirrus  cloud  layer  is  partially  transparent,  and  the  underlying  stratus  and  surface  features  are  visible 
the  scene. 


across 


5.  EARTH  UMB  SCENE  GENERATION 

Generation  of  limb  scenes  by  PACEOS'"^  require  some  differences  in  the  setup  and  running  of  MODTRAN.  The  scene  iliustrated 
h^e  is  calculated  using  the  non-homogeneous  mode  to  select  grid  points  at  which  the  MODTRAN  atmospheric  conditions  and 
illumination  and  viewing  geometries  are  specified.  A  5x5  set  of  grid  points  is  specified  and  setup  to  run  at  five  tangent  heights  with 
five  horizontal  variations  across  the  scene,  with  the  line-of-sight  passing  through  the  entire  atmosphere  from  space  to  space.  Ra¬ 
diances  calculated  for  each  grid  point  are  interpolated  to  form  a  5 12x5 12  pixel  image  representing  the  radiance  at  the  sensor  aper¬ 
ture. 

The  five  tangent  heights  used  were  1, 5, 10, 20,  and  40  km.  Horizontal  variations  are  simulated  by  choosing  the  Subarctic  \Wnter 
atmospheric  model  for  the  right  side  of  the  image  and  the  Tropical  atmospheric  model  for  the  left  side  of  the  image.  The  three 
profiles  specified  between  these  two  are  the  Subarctic  Summer,  Mid  Latitude  Winter,  and  Mid  Latitude  Summer.  This  simulation 
is  run  at  11  pm.  Changing  the  illumination  location  is  not  required  since  at  this  wavelength  the  majority  of  the  radiance  is  from  the 
surface  and  atmospheric  emittance  and  die  effects  of  the  solar  scattering  are  minimal.  No  clouds  are  added  to  this  scene.  The  calcu¬ 
lated  radiance  at  the  sensor  aperture  is  shown  in  Figure  5. 1 . 
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Figure  5.1.  Simulated  11  jJm  Limb  Scene 

6.  SUMMARY  AND  PLANNED  FUTURE  DEVELOPMENTS 

This  paper  illustrates  some  aspects  of  one  of  flie  modules  in  die  PACEOS™  software.  MODTRAN  is  used  to  calculate  radiances 
and  transmission  for  electro-optical  system  perfonnance  assessment  As  illustrated  by  the  scene  generation  process  shown  in  this 
paper,  a  wide  variety  of  scenes  can  be  generated  using  the  Scene  Generation  module  trfPACEOS™.  An  example  of  u^ing  PACEOS™ 
to  generate  simulated  data  for  a  1.6  pm  band  and  illustrations  of  additional  information  that  could  be  extracted  if  a  1.6  pm  band 
were  added  to  AVHRR  is  being  prepared  for  publication.  Hie  simulation  runs  rather  quickly;  and  the  generation  of  all  the 
presented  here  took  less  than  a  day  to  setup  and  run. 

The  accuracy  of  the  scenes  generated  by  PACEOS™  have  been  evaluated  by  simulation  of  operational  systems  such  as  LANDSAT 
and  TIROS  and  by  extracting  scene  parameters  with  EDR  algorithms  such  as  Sea  Surface  Temperature  and  comparing  them  to  the 
known  inputs.  Scenes  simulated  using  MODTRAN  in  conjunction  with  PACEOS™  have  been  found  to  be  of  sufficient  accuracy 
to  evaluate  the  effects  of  modifications  to  various  remote  sensing  system  and  sub-system  elements.  The  use  of  the  PACEOS™  simu¬ 
lation  software  within  Lockheed  Martin  and  in  support  of  other  oiganizations  is  increasing  as  the  modelling  capability  of  the  soft¬ 
ware  is  further  developed. 

PACEOS™  was  used  in  conjunctirai  with  MOD11RAN3  for  these  calculations.  Currently.  MOD’LRANS.S  version  1.2  is  Ivmg  inter¬ 
faced  to  PACEOS™  and  will  replace  MODIRANS.  In  addition,  new  interfaces  and  input/output  options  are  planned  for  PACEOS™ 
so  that  MODTRAN  calculations  can  be  made  and  scenes  generated  for  other  applications  with  more  fidelity  and  accuracy.  Other 


planned  upgrades  to  the  Scene  Generation  pardon  of  PACEOS™  include:  modeling  of  cloud  shadowing;  more  efficient  means  to 
simulate  hyper-spectral  systems;  more  effidait  means  to  account  for  surface  elevation  effects;  and  the  use  of  bi-directional  reflec¬ 
tances  for  the  materials. 
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AF  Optical  Transmission  Meeting 
Phillips  Laboratory  -  Hanscom  AFB  MA 
10-12  June  1997 


Abstract 


Introduction.  A  visible  transmissometer  was  designed  and  constructed  at  Penn  State  by 
the  authors.  During  April  1997,  this  instrument  was  used  to  take  transmittance  data  at  the 
Scripps  Institution  of  Oceanography,  as  part  of  the  Navy’s  Electro-Optical  Propagation 
Assessment  in  the  Coastal  Environment  (EOPACE)  program.  The  goal  of  EOPACE  is  to 
improve  the  understanding  of  surf-generated  aerosols  and  of  electro-optical  propagation  in 
coastal  regions. 

Hardware.  A  feedback-controlled  constant  current  source  drives  a  halogen  bulb,  which 
serves  as  the  white  light  source  (approximated  by  a  3000K  black  body).  The  light  is  then 
mechanically  chopped  and  filtered  to  eliminate  near-IR,  and  is  collimated  by  the  source 
telescope.  A  receiving  telescope  collects  some  of  the  light  at  the  far  end  of  the  optical 
path,  focusing  it  through  a  second  IR  blocking  filter  onto  a  blue-enhanced  silicon  detector. 
The  received  signal  is  then  amplified  by  the  receiver  and  extracted  from  DC  and  noise  by 
the  lock-in  amplifier.  The  lock-in  uses  a  reference  signal  from  the  chopper  to  perform  this 
extraction,  so  a  coax  cable  runs  from  the  source  to  the  receiver.  The  lock-in  then  digitizes 
the  signal  and  transmits  it  to  the  computer  via  RS-232  serial  interface.  The  computer 
processes  and  stores  the  data  (2  samples  per  second)  using  Lab  View  software. 

System  layout.  The  diagram  shows  the  locations  of  the  source  and  receiver  units  during 
the  intensive  operating  period  (lOP).  The  264m  optical  path  included  the  entire  surf  zone. 
Note  that  the  elevations  shown  are  approximate  maximum  distances  above  the  ocean 
surface  (based  on  mean  lower  low  tide),  and  the  distances  were  less  during  most  of  the 
lOP. 

Data  analysis.  The  next  viewgraph  is  an  example  of  the  raw  data  obtained  from  the 
instrument.  This  graph  shows  the  data  taken  on  April  7  between  1 :36am  and  midnight. 
The  vertical  axis  is  the  signal  level  in  volts,  and  the  horizontal  axis  is  arbitrary  data  point 
number  (0  to  163,000).  Note  that  there  are  many  “dropouts”  in  the  data;  these  were 
caused  by  interference  picked  up  by  the  350m  coax  cable.  By  taking  the  maximum  often 
consecutive  points  (5  seconds),  then  averaging  twelve  consecutive  maxima,  most  of  the 
dropouts  were  eliminated. 

A  wide  variety  of  meteorological  data  during  the  lOP  was  provided  by  Scripps  and  NRad, 
including  air,  sea  surface,  and  bottom  temperatures,  relative  humidity,  wind  speed  and 
direction,  wave  height  and  period,  tides,  and  barometric  pressure.  Data  firom  a  vertical 
array  of  temperature  sensors  was  also  available  from  another  Perm  State  experiment. 
These  parameters  were  compared  with  the  transmittance  data  in  an  attempt  to  find  a 
correlation.  The  best  correlation  found  was  with  wind  speed  and  direction.  On  April  7 
and  8,  a  moderately  strong  west  wind  (10-15  knots)  was  present  during  much  of  the 
daylight  hours,  and  correlates  with  the  times  of  low  transmittance.  On  April  5,  data 
collection  began  at  about  2pm,  and  transmittance  remained  roughly  constant  through  the 
remainder  of  the  day.  On  this  day,  the  wind  speed  was  only  about  7  to  8  knots.  It  is 


337 


theorized  that  the  west  wind  is  either  bringing  in  aerosols  from  the  open  ocean,  or 
affecting  the  surf  plumes  in  such  a  way  as  to  reduce  transmittance. 

Preliminaiy  modeling  (MODTRAN).  The  predicted  transmission  versus  wavelength 
curves  for  the  264km  path  were  plotted  using  four  standard  MODTRAN  aerosol  models: 
Rural  23km  Visibility,  Rural  5km  Visibility,  Maritime  23km  Visibility,  and  Navy  Maritime 
23km  Visibility.  Next,  the  combined  spectral  responses  of  the  detector  and  the  two  IR- 
block  filters  were  entered  into  the  filter  function  associated  with  MODTRAN.  The  ‘Tilter 
Function  Input”  viewgraph  shows  the  resulting  wavelength  weighting  curve.  The  black 
body  input  was  set  to  3000K  for  the  halogen  bulb.  Finally,  the  filter  function  was  applied 
to  the  four  models,  yielding  predicted  transmittances  of  82.9%,  95.9%,  95.7%,  and 
98.4%,  respectively.  Returning  to  the  data,  when  the  voltage-to-transmittance  conversion 
is  applied,  it  is  seen  that  the  measured  transmittance  for  the  57  hours  of  data  falls 
predominantly  at  about  65%  and  about  7%. 

Future  work.  Additional  transmittance  data  and  improved  understanding  of  surf  plumes 
will  be  required  in  order  to  better  understand  and  predict  surf  zone  transmittance. 
Comparison  of  data  from  this  visible  transmissometer  with  the  NRad  infrared 
transmissometer  (also  operated  during  this  EOF  ACE  lOP)  may  also  provide  some 
insights. 
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Visible  Transmissometer  Mike  Zugger,  Penn  State  University 


Transmissometer  Layout  for 
EOPACE  at  Scripps 
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Scripps  Pier  UCSD 


Transmissometer  data  -  April  7  (.5sec/all  points) 
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Meteorological  Data  Available 
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Temporal  and  Spatial  Characteristics  of  Surf  zone  Plumes^ 


Christophe  F.  Bas^  and  C.  Russell  Philbrick^ 
The  Pennsylvania  State  University 
Electrical  Engineering  Department 
University  Park,  PA  16802 


Abstract 

One  of  the  goals  of  the  investigation,  Electro-Optical  Propagation  Assessment  in  the  Coastal 
Environment  (EOPACE),  is  to  enhance  our  understanding  of  the  impact  that  surf  generated  aerosols  have 
on  the  propagation  of  optical  radiation  through  coastal  regions.  A  bi-static  LIDAR  system  uses  a  fan 
beam  to  illuminate  the  atmospheric  region  above  the  surf  which  permits  measurements  of  both  the  vertical 
and  horizontal  extent  of  aerosol  plume  structures.  The  scattered  radiation  is  imaged  by  two  sensitive  CCD 
cameras  and  a  digital  imaging  camera,  the  latter  being  available  only  for  the  most  recent  test.  Three 
measurement  phases  of  the  experiment  program  have  provided  data  on  the  surf  generated  coastal  plumes 
at  two  locations  and  under  a  range  of  different  meteorological  conditions.  General  spatial  and  temporal 
characteristics  of  these  "coastal  plumes",  are  described  in  this  report. 

1.  Introduction 

Judging  from  the  interest  in  scientific  investigations  such  as  MAPTIP  [1],  VOCAR  [2],  and  EOPACE 
[3],  the  coastal  environments  are  of  considerable  importance.  The  attraction  of  the  EOPACE  program  is 
its  focus  on  the  propagation  of  electromagnetic  waves  at  optical  frequencies  in  the  coastal  environment. 

Toward  this  study,  ARL/PSU  (Applied  Research  Laboratory  of  The  Pennsylvania  State  University) 
has  provided  measurement  support  using  three  instruments:  (1)  a  LIDAR,  (2)  a  transmissometer  [4],  and 
(3)  a  temperature  mast  (included  in  the  April  1997  experiment).  The  ARL/PSU  lidar  has  obtained  data  at 
the  two  Californian  experimental  sites:  (1)  the  Scripps  Institution  of  Oceanography  (SIO,  La  Jolla)  in 
1996  and  1997  and  (2)  Moss  Landing  (Monterey)  in  1996.  Upon  the  first  visit  to  the  SIO  (22  January  to  5 
February  1996),  tall  vertical  structures  were  revealed  by  the  bi-static  LIDAR  (Figure  1).  These  plumes 
provided  a  visual  depiction  of  the  scatters  which  decrease  transmission.  Hence,  both  tlie  composition  and 
characteristics  of  these  plumes  are  of  interest. 


Figure  1:  Typical  plumes  observed  over  the  surf  zone  at  La  Jolla,  CA,  when  illuminated 
by  a  vertical  Argon-ion  laser  sheet.  The  picture  was  taken  on  2  February  1996,  at  7:15 
PM  with  a  low-noise  CCD  camera.  The  exposure  time  was  10  seconds. 


*  Paper  presented  at  the  Conference  on  Atmospheric  Transmission,  AFGL,  Hanscom,  MA,  June,  1997. 
^  cfbl04@psu.edu 
^  crp3@psu.edu 
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2.  Experimental  set-up 

Figure  2  shows  the  set-up  of  the  third  experiment  which,  like  the  first  one,  was  conducted  at  the  SIO 
(31  March  31  to  11  April  11  1997).  The  unconventional  bi-static  LIDAR  used  a  cylindrical  lens  to 
transform  the  incident  Argon-ion  (514.5  nm)  CW  laser  beam  into  a  vertical  laser  sheet.  The  receivers 
were  low-noise  high-resolution  (16  bits)  CCD  cameras  and  one  digital  video  camera.  The  former  allows 
long  integration  times  to  retrieve  distribution  of  the  scatterers  [5]  and  the  latter  to  reveal  the  dynamics  of 
the  plumes.  The  temperature  mast  was  located  on  the  pier  in  the  middle  of  the  surf  zone  and  provided 
more  details  about  the  processes  just  above  the  interface.  Simultaneously,  the  SIO  kept  a  record  of  the 
meteorological  conditions  from  a  weather  station  at  the  end  of  the  pier.  The  distance  over  which  the 
waves  continuously  break  has  been  observed  to  be  typically  between  200  and  100  meters  from  the  beach 
end  of  the  300  meter  long  pier. 


Figure  2:  Experimental  set-up  at  La  Jolla,  CA.  L  refers  to  the  location  of  the  source  for 
the  laser  fan  (emitter),  A  and  B  refer  to  those  of  the  CCD  cameras,  C  and  D  refer  to 
those  of  the  digital  video  camera,  T  refers  to  the  one  of  the  temperature  mast,  and  R  is  a 
fixed  reference  point  of  known  location  which  appears  in  every  image. 


361 


3.  The  Coastal  Atmospheric  Boundary  Layer 

Considerable  effort  is  currently  underway  to  model  the  dynamics  of  coastal  environments  [6],  Those 
focusing  on  the  first  few  tens  of  meters  above  the  surf  zone  may  validate  their  models  using  the  ARL/PSU 
images.  Figure  3  shows  a  sample  of  a  cropped  image  from  one  of  the  CCD  cameras  and  its  corresponding 
projection  along  the  pier.  By  keeping  track  of  the  location  of  the  plumes  and  their  heights  under  similar 
atmospheric  conditions,  one  can  plot  the  distribution  of  those  heights  versus  the  distance  from  the  origin 
of  the  pier.  Figure  4  displays  such  a  plot  obtained  during  the  night  of  2  April  1997,  night  during  which  a 


Figure  3:  The  left-hand  side  image  shows  a  truncated  image  captured  by  the  CCD 
camera  positioned  at  A  (see  the  map  of  Figure  2)  on  2  April  1997,  at  11:43  PM. 
Exposure  time  was  10  seconds.  The  structure  in  the  background  is  that  of  the  weather 
station  at  the  end  of  the  pier  and  the  bright  dot  on  its  top  is  R,  the  reference  point.  The 
projection  of  this  image  onto  a  vertical  plane  which  is  parallel  to  the  pier  is  shown  on 
the  right-hand  side  plot.  The  intensity  contour  chosen  to  define  a  plume  was  set  to  33 15 
digital  counts. 


Figure  4:  Plume  height  versus  distance  from  the  origin  of  the  pier  during  the  night  of  2 
April  1997.  The  tilted  dashed  lines  trace  the  aperture  of  the  laser  fan.  The  solid  line 
corresponds  to  a  least-square  fit  of  the  data  points  (circles)  with  the  following  function: 
Height  =  Const^Distance  -  Origin  • 
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slight  offshore  wind  was  present.  As  the  air  mass  advected  from  land  glides  over  the  warmer  surf-zone,  a 
Thermal  Internal  Boundary  Layer  (TIBL)  develops  [7].  Several  equations  exist  to  model  a  TIBL,  but  most 
of  them  can  be  parameterized  as: 

Height  =  Const^Distance  -  Origin  (1) 

where  the  value  of  Const  depends  on  the  environmental  conditions.  The  offshore  wind  acts  as  a  sink  for 
the  net  flux  of  particles  emanated  from  the  breaking  waves.  These  conditions  typically  existed  from  late 
evening  to  early  morning  during  the  experimental  phase  (Figure  5a).  When  the  temperature  difference 
between  water  and  land  was  small,  the  wind  decreased  in  intensity  and  the  density  of  the  scatterers  is 
approximately  symmetrically  distributed  over  the  surf  zone.  Figure  1  shows  such  a  case  where  the  plumes 
are  standing  vertically.  This  particular  situation  occurred  during  transitional  periods  corresponding  to  a 
shift  in  direction  of  the  mean  wind  (from  offshore  to  onshore  and  vice-versa)  and  are  depicted  in  Figures 
5b  and  5d. 

When  land  temperature  is  greater  than  water  temperature,  onshore  wind  results.  With  it  comes  an 
TIBL,  but  located  at  the  land  and  water  intersection.  Documentation  of  a  possible  internal  boundary  layer 
generated  by  the  change  of  surface  roughness  between  the  open-ocean  and  the  surf  zone  has  yet  to  be 
satisfactorily  supported.  One  reason  is  due  to  the  geometry  of  the  bi-static  lidar  whose  vision  of  the 
plumes  in  their  fall  height  is  impossible  close  to  the  transmitter  (see  dashed  lines  in  Fig.  4).  We  suspect 
that  an  TIBL  exists  over  the  entire  ocean  surface,  but  that  it  is  only  discernible  when  the  aerosol  from  the 
surf  provides  a  tracer  which  can  be  observed  from  the  laser  illumination.  The  distinctive  feature  which  is 
believed  to  form  the  TEBL  is  the  temperature  gradient  in  the  first  several  meters  above  the  surface,  which 
is  governed  by  the  air/sea  temperature  difference.  If  this  hypothesis  is  true,  then  the  same  kind  of  plume 
characteristics  will  be  present  during  high  sea  states  or  when  the  aerosols  are  produced  in  the  wake  of  a 


Surf-zone 

Figure  5:  Cyclic  temporal  changes  over  the  surf  zone,  (a)  Over  night  (surface 
temperature  of  water  is  greater  than  the  one  over  land),  (b)  early  morning  after  sun  rise 
(surface  temperature  of  water  and  land  are  approximately  equal),  (c)  later  into  the 
morning  (surface  temperature  of  water  is  smaller  than  the  one  over  land),  and  (d)  early 
into  the  night  (surface  temperature  of  water  and  land  are  approximately  equal). 


363 


Table  1:  Summary  of  the  mean  meteorological  conditions  and  resulting  plume  heights 
over  the  surf  zone  during  the  three  lOPs,  _ _ 


Location 

SIO(D 

Moss  Landing 

SIO  (ID 

Ts..-T.ir  [C] 

3.3 

1.7 

2 

3.1 

1 

1.2 

PAU  il  1  ILiJ  ilJliJ  1  Hi  M  iMl 

25 

10 

11 

ship  in  the  presence  of  air/sea  temperature  difference.  Favorable  conditions  for  the  development  of  the 
coastal  TBL  were  between  late  morning  and  early  evening  (Fig.  5c).  This  time  frame  renders  the 
observations  of  the  plumes  even  more  challenging  because  of  the  presence  of  the  extremely  large  optical 
background  signal  during  daytime. 

Table  1  reviews  the  mean  meteorological  conditions  in  which  all  three  experimental  phases  were 
conducted  [8].  Departures  from  those  mean  values  (all  else  being  equal)  reveal  interesting  features.  First, 
when  the  wind  increases,  the  plume  height  decreases.  Second,  when  the  temperature  difference  between 
the  sea  and  air  decreases,  so  does  the  plume  height.  This  last  effect  agrees  with  Equation  1  for  which 
Raynor  et  al.  and  Venkatram  [7]  have  parameterized  the  proportion  constant  as  follows: 

Const  «  -  r^^l^/wind  (2) 

We  have  observed  that  the  tallest  plumes  occurred  with  calm  winds  and  large  temperature  gradients. 

We  have  argued  that  the  plumes  scale  with  the  depth  of  an  internal  boundary  layer.  This  gives  us  a 
better  idea  of  the  atmospheric  conditions  which  will  likely  entrain  the  scatterers,  provided  that  the  wave 
activity  produces  them  in  sufficient  quantity  to  render  the  plumes  visible.  A  reliable  source  of  ejected  and 
spayed  droplets,  i.e.,  scatterers,  are  the  breaking  waves.  Because  not  all  breaking  waves  produce  the  same 
amount  of  aerosols,  two  of  their  categories  are  of  interest:  plungers  and  spillers.  As  waves  turn  into 
plunging  breakers,  trapped  pockets  of  air  entrained  below  the  water  surface  eject  many  aerosols  (i.e., 
water  droplets)  as  the  bursting  air  bubbles  reach  the  surface.  Spilling  breakers  are  more  gentle,  surface 
bursting  bubbles  remain  mainly  at  the  surface.  Under  similar  atmospheric  conditions,  plunging  breakers 
produced  more  scattering"^  than  spilling  ones. 

Hence,  three  variables  seem  to  be  appropriate  to  predict  the  height  of  the  plumes:  wind  speed, 
temperature  difference,  and  breaker  type. 

4,  Temperature  Mast 

One  way  to  assess  whether  atmospheric  stability  influences  the  height  reached  by  the  plumes  is  to 
measure  the  thermal  vertical  structure  of  the  region  above  the  surf  zone.  This  aim  motivated  the 
installation  of  the  mast,  equipped  with  7  thermistors,  in  the  middle  of  the  surf  zone  (Fig.  2).  Figure  6 
sketches  the  mast  and  provides  the  height  for  each  sensor  relative  to  the  pier  level.  Once  the  data 
gathered,  another  analysis  has  been  undertaken  because  of  the  necessity  to  refine  and  enhance  some  of  the 
data  sets  required  to  estimate  the  Richardson  number  over  the  surf  zone.  Again,  it  is  necessary  to  note 
that  the  height  of  the  plumes  is  affected  by  the  intensity  level  chosen  to  define  a  “plume”.  This  number 
depends  on  the  intensity  of  the  scattering  by  the  plumes  and  the  exposure  time  to  capture  them. 

Relative  comparison  between  the  environmental  characteristics  and  the  mean  height  of  the  plumes 
consists  of  the  computation  of  (1)  the  mean  temperature  and  the  standard  deviation  of  the  mean 
temperature  over  all  the  sensors  at  each  point  in  time  (this  is  equivalent  to  assuming  that  each  sensor  sees 
one  realization  of  the  same  random  variable)  and  (2)  the  time-average  for  both  the  mean  and  the  standard 


We  imported  the  images  into  MATLAB  where  the  function  contour  delimited  the  plumes.  The 
subjective  part  of  the  analysis  was  to  define  the  intensity  levels  corresponding  to  what  we  defined  as  being 
“plumes”.  Once  a  number  had  been  selected  (for  a  given  exposure),  it  became  a  reference  for  the  other 
cases. 
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Mean  Flime  height  versus  mean  tenrpersture  standard  delation 


Figure  6:  The  temperature  mast  was  located  in  the 
the  middle  of  the  surf  zone  where  it  captured  seve¬ 
ral  profiles.  Heights  are  referenced  to  the  pier. 
Water  level  was  approximately  -10  meters. 
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Figure  7:  Mean  height  of  the  plumes  for  2-10  April, 
1997,  as  a  fimction  of  temperature  standard  devia¬ 
tion.  With  a  10  sec  exposure,  squares  refer  to  a  33 15 
contour-level  and  triangles  to  one  of  4200.  With  a  30 
sec  exposure,  diamonds  refer  to  a  4200  contour-level 
and  circles  refer  to  one  of  3500. 


deviation  of  the  mean  temperature  (which  is  equivalent  to  invoking  ergodicity  to  expand  the  total  number 
of  realizations).  Table  2  summarizes  the  results  of  this  processing.  Figure  7  plots  the  plume  height  versus 
the  standard  deviation  of  the  mean  temperature  for  this  third  experiment. 

Before  formulating  hasty  conclusions  on  the  problem  at  hand,  it  is  essential  to  realize  that  Table  2 
allows  comparison  within  the  same  group  only.  More  work  is  required  to  compare  results  between 
different  groups. 

First,  both  the  diamonds  and  the  triangles  reveal  that  the  taller  the  plumes  the  higher  the  standard 
deviation  of  temperature.  Second,  the  mean  temperature  results  from  a  complicated  process  which  takes 


Table  2:  Summary  of  the  results  fi*om  the  processing  of  the  plume  heights  and  the  data 
taken  on  the  mast  of  temperature  sensors. _ _ 


Day/Month/ 

Year 

Time  (PST) 

Symbol 
(Fig- 7) 

Mean  plume 
height  [m] 

Mean  tempe¬ 
rature  fC] 

Mean  temperature 
std.  deviation  [C] 

23:35-23:38 

13.3 

11.0 

0.29 

23:54-00:00 

circle 

13.6 

10.9 

0.33 

19:50-19:55 

diamond 

7.6 

16.7 

0.45 

21:43-21:50 

diamond 

8.3 

16.1 

0.51 

KES33EEI 

diamond 

12.1 

16.9 

0.65 

BES3SEEI 

7.6 

17.0 

0.81 

i^bssesi 

20:21-20:24 

7.4 

16.5 

0.45 

^  The  plumes  were  analyzed  over  this  time,  but  the  closest  set  of  temperature  data  was  taken  on  8  April 
from  23:54  to  midnight.  The  meteorological  and  oceanographic  conditions  were  stationary  from  one  time 
segment  to  the  next. 

^  The  plumes  were  analyzed  at  this  time,  but  the  closest  set  of  temperature  data  was  taken  from  01:19  to 
01:25.  The  meteorological  and  oceanographic  conditions  were  stationary  from  one  time  segment  to  the 
next. 
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into  account  both  the  air  and  the  droplets  inside  the  plumes.  The  temperature  sensors  located  upwind  and 
outside  the  surf  zone  give  an  actual  air  temperature  (relatively  free  of  surf-zone  generated  droplets) 
considerably  lower  than  the  mean  temperature  measured  by  the  sensors  on  the  mast  (by  2  degrees 
Celsius).  The  only  other  potential  thermal  source  is  the  water  itself  whose  "bottom"  temperature  (defined 
by  the  SIO  as  being  3.5  meters  below  the  mean  lower  low  water  level)  remains  around  17  degrees  Celsius. 
Hence,  the  measured  temperature  results  from  a  mixture  of  air  and  water  droplets,  the  latter  of  which 
interacts  with  the  thermal  sensors  through  the  processes  of  deposition  and  evaporation.  This  result 
confirms  that  the  content  of  the  plumes  is  sea  water  in  nature.  Contamination  of  temperature  sensors  by 
water  has  been  reported  by  Friehe  [9]  who  has  noticed  the  spiky  signature  of  the  measurements. 

An  increase  in  the  temperature  fluctuations  requires  that  there  is  “enough”  time  for  some  of  the  water 
to  evaporate  and  that  there  is  a  significant  amount  of  water  coming  suddenly  into  contact  with  the 
thermistors.  The  physical  processes  which  promote  this  behavior  need  to  be  further  analyzed. 

Out  of  curiosity,  the  analysis  of  the  plume  heights  for  a  different  intensity  level  has  been  investigated. 
Once  this  level  set  to  4200  (the  same  as  for  the  triangles  and  diamonds),  neither  the  circle  nor  the  square 
of  Figure  7  resolve  plumes.  Notice  that  their  standard  deviations  are  the  smallest  of  all  the  data  points. 

Therefore,  we  conclude  that  there  is  a  strong  correlation  between  the  mean  standard  deviation  of 
temperature  and  the  mean  plume  height.  However,  the  mean  temperature  is  not  a  constant  indicator  of 
the  plume  height. 

5.  Conclusion 

Our  study  shows  that  an  internal  boundary'  layer  provides  a  means  of  estimating  the  height  of  the  surf 
zone  plumes.  The  physics  behind  such  a  boundary  seems  to  include  wind  speed,  buoyancy  forcing,  and, 
possibly,  changes  in  surface  roughness.  The  study  of  boundary  layers  is  rendered  difficult  by  the  nesting’ 
of  boundary  layers  over  the  surf  zone. 

The  meteorological  conditions  seem  to  prepare  the  atmosphere  to  accept  any  aerosols  to  form  plumes, 
but  the  major  contributor  of  aerosols  is  the  surf  zone  itself  via  its  breaking  waves  which  eject  significant 
amounts  of  water  droplets.  The  brightest  and  tallest  plumes  are  observed  in  the  presence  of  larger 
temperature  difference  between  air  and  water,  minimum  wind,  and  plunging  breakers.  Further  studies  of 
the  effect  of  turbulence  is  necessary  to  determine  whether  it  promotes  tall  plumes  or  shreds  them.  Our 
next  step  is  to  simulate  those  effects  using  Large  Eddy  Simulations  [10]  and  compare  the  outputs  of  the 
model  to  our  experimental  data. 
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’  If  one  accounts  for  all  the  phenomena  generating  their  own  boundary  layer,  one  faces  an  infinite  number 
of  intricate  boundary  layers,  some  developing,  others  blending  in,  but  several  boundary  layers  coexist, 
especially  over  the  surf  zone. 
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Evaluation  of  atmospheric 
models  and  characterization 
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Special  radiance  calculator  convolves  material  spectra  with  radiative  environment 
and  sensor  spectrai  response  function  to  simuiate  at-sensor  radiance  or 
reflectance 
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HYPEX  Atmospheric  Correction 
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Coincident  coiiection  of  cloud,  ground,  and  , 
atmospheric  "truth" 
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24.34N,  81.41  W,  6:47  GMT,  29  Feb  1996,  Tropical  Winter  Model  Atmosphere 
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Mon  Jun24  11:27:06  166G 


Adjustment  of  Temperature  Profile 
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Transmission  Differences 

Adjustment  of  Temperature  Profile  +  8  degrees 
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WavAldngih  (micron*)  Wavelength  (microna)  Wavelength  (microna) 
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Photon  Research  Associates,  Inc. 
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HYPEX  Derived  Reflectances 
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Simulated  responses  for  sensor  design 


MICROWAVE  TRANSMISSION 
IN  LIQUID  WATER  CLOUDS: 
UNCERTAINTIES  AT  LOW  TEMPERATURES 


ALAN  E.  UPTON,  MICHAEL  K.  GRIFFIN 

OPTICAL  EFFECTS  DIVISION 
GEOPHYSICS  DIRECTORATE 
AIR  FORCE  RESEARCH  LABORATORY 


AND  ALAN  G.  LING 
RADEX,  INC. 

BEDFORD,  MASSACHUSETTS 


OUTLINE 

*  MICROWAVE  ATTENUATION  MODELS 

*  TEMPERATURE  AND  FREQUENCY  DEPENDENCES 

•  CLOUD  vs.  TOTAL  ATTENUATION 

•  EFFECT  OF  ATTENUATION  UNCERTAINTY  ON 

-  SATELLITE  BRIGHTNESS  TEMPERATURES 

-  LIQUID  WATER  CONTENT  RETRIEVALS 
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MICROWAVE  ATTENUATION  MODELS 
EOR  CLOUD  LIQUID  WATER 


USE  RAYLEIGH  ABSORPTION  APPROXIMATION 

CLOUD  DROPLETS  ARE  SMALL,  RELATIVE  TO 
WAVELENGTH  OF  RADIATION 

ATTENUATION  DEPENDS  ON  LIQUID  WATER  CONTENT, 
NOT  ON  DROPLET  SIZE  DISTRIBUTION 

DIELECTRIC  CONSTANT  OF  WATER  IS  REQUIRED 
RAY  (1972)  MODEL 

-  WIDELY  USED  STANDARD  (INCL  RADTRAN,  EASE) 

LIEBE,  HUFFORD,  AND  MANABE  (1991)  MODEL 

-  TWO  RELAXATION  TERMS 
"DOUBLE  DEBYE"  METHOD 

-  TWO  OPTIONS  FOR  TEMPERATURE  DEPENDENCE 
OF  FIRST  RELAXATION  FREQUENCY 

NO  DATA  BELOW  0°C! 
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TWO  MODELS  OF  PRIMARY  RELAXATION  FREQUENCY 
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Temperature  (°C) 


Frequency  (GHz)  Frequency  (GHz) 


RAY  -  LIEBE  BRIGHTNESS  TEMPERATURE  DIFFERENCE 
U.S.  STANDARD  ATMOSPHERE  WITH  CLOUD  BASE  =  4.4  km,  TOP  =  5.4  km  (-ZO^C),  MAX  LWC  =  0.3  g/m 
ZENITH  ANGLE  =  55°,  OCEAN  BACKGROUND  WITH  WIND  SPEED  =  5  m/s 
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100  150  200  250 

Frequency  (GHz) 


CLOUD  LIQUID  WATER  RETRIEVALS  FROM  SIMULATED  SSM/I  DATA 

ALISHOUSE  ET  AL.  (1990)  REGRESSION  FORMULA 
U.S.  STANDARD  ATMOSPHERE  WITH  CLOUD 
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True  Ray  Liebe  Ray 

Cloud  Top  =  — 20°C  Cloud  Top 


CONCLUSIONS 


THERE  ARE  LARGE  UNCERTAINTIES  IN  MICROWAVE 
ATTENUATION  BY  LIQUID  WATER  AT  LOW  TEMPERATURES 

THESE  UNCERTAINTIES  IMPLY  SIGNIEICANT  ERRORS 
IN  SATELLITE-RETRIEVED  CLOUD  PROPERTIES 

LABORATORY  AND  MODELING  RESEARCH  IS  NEEDED 
TO  REDUCE  THE  UNCERTAINTIES 


396 


ON  ERA 


ECRiN96 

(Etude  dll  Champ  Radiatif  Infra  rouge  des  Niiages) 


P.  Simoneau,  C.  Malherbe.  A.  Boischot  G.  Durand.  B. Rosier.  A.  Delannoy 

ONERA.  France 

P.  H.  Flamand,  L.  Sauvage.  H.  Dirand 
LMD.  France 


J.  Felon 
SA.  France 
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ECRIN  96 


Objectives 


1/  Code  validation  : 

-  NUALUM 

-  MONET 

2/  Radiance  fluctuations  on  water  ciouds 
3/  Microphysic  data  base 
4/  Cirrus  radiance  measurements 
5/  Restitution  of  atmospheric  radiatives  fluxes 
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3 


Measurement  strategy 


ARAT/Fokker-27 

Caravelle  116 

Merlin  IV 

Ceiling  altitude 

Cruising  airspeed 

Endurance 

5500  m 

80-100  ms-1 

3h-4  hours 

6000  m 

200-300  Kni 

5  hours 

7  000  m 

100-130  ms-1 

4-5  hours 

Instrumentation 

-  backscatter  lidar 

LEANDRE-1 

-  pyranometers, 

pyrgeometers, 

Barnes  PRT-5 

-  in  situ  probes 

-  IR  spectrometer 

-  IR  imager 

1 

-  pyranometers, 

pyrgeometers, 

Barnes  PRT-5 

-  in  situ  probes 

-  microphysics  probes 

ECRIN  96 


Measurement  strategy 


-  Satellite  information  :  AVHRR,  SPOT 

-  Radiosounding 


Radiance  measurements  : 

^  spectrometer,  IR  imageur 


In  Situ  measurements  : 

■»  LWC/t/U  probes.  FSSP, 

Cirrus  occurence,  aerosols  : 


Lidar 
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ONER  A 


Caravelle  116 


Cloud  and  atmospheric  parameters 

Instrument 

Carrier 

-  CVF  spectrometer 

'  spectral  range  1.5  -  5.3  pm 

-  spectral  resolution  AXyX  =  IP/c 

“  liquid  nitrogen  cooled 

-  fov  =  3  mrd 

- 10  spectra/s 

Spectrometer  SICAP 

Caravelle 

-  PtSi  camera 

“  waveband  3 -5pm 

-  fov  =  1.5  ^  14  ®  (depending  lens  focal  length) 

-  25  images/s 

Imager  CERAP 

Caravelle 
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3 


MeV'xc'w  ^  ti-r 


ONERA 

MERLIN  4 


Cloud  and  atmospheric  parameters 

Instrument 

Carrier 

-  down-  ( )  and  upward  shortwave  ( )  fluxes 

Pyranometer 

Merlin  IV 

-  down-  ( )  and  upward  longwave  ( F^^^ )  fluxes 

Pyrgeometer 

-  thermodynamic  :  temperature  (T),  pressure  (p),  and 

humidity  (q)  profiles 

-  microphysics  :  liquid  water  content  (LWC),  particles  size, 

concentration 

In-siiu  probes 

FSSP 

Merlin  IV 
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-9,  obcm 


ARAT 


Cloud  and  atmospheric  parameters 

Instrument 

Carrier 

-  cloud  base  and  height,  geometrical  thickness  AH, 

Backscatter  Lidar 

ARAT 

-  optical  depth  at  0.53  pm  (5vis), 

-  vertical  structure  of  extinction  (a)  and  backscaiier  (p) 

coefficients  at  0.53  pm, 

-  scattering  phase  function  p(7t)  or  backscatter  to  extinction 

ratio  (k),  p{0)  =  ^7tk 

-  ice/water  discrimination,  vertical  structure  of  linear 

depolarization  ratio  (A) 

LE  ANDRE- 1 

-  down-  ( )  and  upward  shortwave  ( F.|.^  )  fluxes 

Pyranometer 

-  down-  )  and  upward  longwave  ( )  fluxes 

Pyrgeometer 

-  down-  and  upward  emissivity  (e),  infrared  optical  depth 

(Sir) 

PRT-5 
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ON  ERA 


ECRIN96 


Date  in 
October  1996 

Mission  # 

Take-off 

Duration 

Objectives 

9 

1 

3.5  hours 

low  clouds 

10 

2 

3.6  hours 

low  clouds 

11 

3 

2.3  hours 

cirrus 

14 

4 

2.4  hours 

edges 

15 

5 

2.6  hours 

SPOT 
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E2_V101.xls 
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Luminances  Int6gr4es  (2) 
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Niiage_2_Vis6e-05.xls  Page  1 


Heure :  13:00:16.075 

Status :  Dteignation  Spectres  SICAP 

Num6ro(s):  2250  -  (R6f :  2774  h  2786) 


Longveur  d'onde  (fim) 


QMERA 
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A20 


ECRIN’96/voI38/pointA1/12h51-13h00  (48  km) 


Ld 
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Variation  de  I’altitude  du  sommet  (m) 


LccJiXit  t 
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Epaisseur  optique 


Nombre  d’apparitions 


ECRIN’96/vol38/pointA3/1 4h08-1 4h1 8 


Epaisseur  optique 


-c'cXcui  r 
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ECRIN  96 


Scientific  developments 


code  validations  (ID,  3D) 

3D  cartography 
comparison  Udar/ln  Situ 

comparison  in  situ  measurements/ AVHRR  informations 
(altimetry,  microphysic) 

fluxes  restitution 


publication  (B/WIS) 


ONERA 


ECRIN  96 


Conclusion 


3  flights  with  measurements  on  Sou  : 

-  radiance  (spectra,  images) 

-  in  situ  microphysic 

-  lidar  information  (cirrus,  boundary  layer) 

-  AVHRR  images 

-  radiosounding 


1  flight  with  measurements  on  CIRRUS  : 

-  radiance  (spectra,  images) 

-  lidar  information  (k^J 

-  AVHRR  images 

-  radiosounding 


1  flight  with  SPOT  appointment 
-  presence  of  ice 


Measurements  of  land  and  sea  radiance 

3 


A  lot  of  data  and  a  lot  of  WORK  !! 
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MIT  Lincoln  Laboratory 


Methodology  for  Cloud  Determination 

Based  on  Humidity  and  Temperature  Prof  lies 
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j^:  Latent  heat  of  sublimation  2834  J/gm,  R^:  0.4615  J/  gm-K 


Opaque  Cloud/ Transmissive  Cloud/  No  Cloud 


MIT  Lincoln  Laboratory 


A  Sample  No-Cloud  (Clear)  Case 


Relative  Humfdiij'  (K)  Temperature  (K) 

i— MIT  Lincoln  Laboratory 


A  Sample  Low  Cloud  Case 
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Relative  Hurrifdlt^  {K)  Temperature  (K) 

MIT  Lincoln  Laboratory 


Sample  Middle  Cloud  Case 


ZD  40  ao  SO  ^  1 0O  200  2Z0  240  260  Z80 

Rflla^lve  HuiriFdltj^  {%)  Temperature  (K) 

MIT  Lincoln  Laboratory 


A  Sample  High  Cirrus  Case 


40  «0  SO  1 0O  200  220  240  260  280  300 

Rfilatrlve  Humrdit^  (K)  Temp-einature  (K) 

MIT  Lincoln  Laboratory 


Diagnostic  Test  for  Vertical  Extent  of 

Cirrus  Cloud 


MIT  Lincoln  Lai 
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SURFACE  MEASURED  fC) 


SURFACE  TEMPERATURE  RETRIEVAL 
(NIGHTTIME,  OCEAN,  TRANSMISSIVE  CLOUD) 


(Oo)  iss  aaAiyaa  aimaivs 


BUOY  MEASURED  (°C) 


MWIR  Atmospheric  Correction  Algorithm 
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MODTRAN4  and  EASE:  Potential  for  Oimate-Change  Modeling 


G.P.  Anderson,  J.H.  Chetwynd 
PL/Geophysics  Directorate,  Hanscom  AFB,  MA  01731 
781-377-2335,  Fax 377-8900,  ‘ganderson@plh.af.mil’ 

A.  Berk,  L.S.  Bernstein,  P.K.  Acharya 
Spectral  Sciences,  Inc.,  Burlington,  MA  01803 
781-273-4770,  Fax 270-1 161,  ‘lex@spectral.com’ 

H.E.  SneU 

Atmospheric  and  Environmental  Research,  Cambridge,  MA 
617-549-2281,  FAX 617-661-6479,  ‘snell@aer.com’ 

1.  INTRODUCTION: 

“Hyperspectral”  is  a  relatively  new  descriptor  for  measurement  techniques  historically  used  by  the 
atmospheric  remote  sensing  communities.  As  opposed  to  “multi-spectral”  sensing,  which  includes  any  instrument 
with  a  finite  number  of  specific  channels,  filters,  or  bands,  hyperspectral  implies  employing  raough  (usually 
contiguous)  channels  to  provide  redundant  information  on  both  the  characteristics  of  die  sui&ce  and  die  atmosphere 
when  die  instrument  is  configured  for  nadir  viewing.  This  definition  (Alex  Goetz,  Univ.  of  Colorado,  private 
communication)  can  be  contrasted  with  others,  such  as  ‘any  instrument  with  better  than  4  cm’*  resolution.”  This  latter 
definition  can  be  inadequate  in  the  infiared  spectral  range  because  many  overlapping  molecular  systems  will  not  be 
sufficiendy  discriminated  at  such  ‘narrow’  resolution.  However,  the  4  cm’*  is  ofto  much  more  than  adequate  in  the 
visible,  when  typical  resolutions  of  10  ran  (e.g.  AVIRIS,  Green  et  al.,  1996)  at  600  ran  corresponds  to  a  fi'equency 
resolution  of  over  200  cm’* . 

Widi  die  development  of  MODTRAN4  (Bernstein,  et  al.,  1996),  a  flexible  tool  for  applications  of 
realistic  hyperspectral  sensing  is  now  available.  While  this  version  of  MODTRAN  is  grounded  in  the  prior  series  of 
AF  radiative  transfer  band  model  (BM)  algorithms  (LOWTRAN,  Kneizys,  et  al.,  1980,  1983,  1988,  throu^ 
MODTRAN3,  Berk  et  al.,  1989,  1995),  it  is  distinct  in  its  ability  to  employ  Beer’s  law  (T,^  exp(-kv,,-  n,)  to  describe 
local  layer,  species-specific  transmittance  for  input  to  the  radiance  calculations.  While  this  capability  is  not  always 
necessary,  it  allows  ^propriate  handling  of  multiple  scattering  (m.s.)  using  casting  non-BM  ^goritfams  (DISORT, 
Stamnes  et  al.,  1988,  and  Isaacs  et  al.,  1987).  MODTRAN4,  while  maintaining  the  basic  2  cm’*  spectral 
resolution,  can  now  complement  the  m.s.  routines  by  introducing  a  Correlated-^  (CK)  capability  which  is 
expressly  compatible  with  Beer’s  law  formulations.  MODTRAN4  also  provides  greatly  improved 
predictive  capabilities  under  cloudy  and/or  heavy  aerosol  loading  conditions  in  both  the  visible  and  IR  by 
allowing  the  explicit  definition  of  water  and  ice  cloud  vertical  profiles  and  spectral  data,  either  by  scaling 
and  combining  default  model  clouds  or  by  redefining  entirely  new  model  clouds  with  micro-layering 
options. 

It  can  be  shown  that  this  combination  of  improvements  will  permit  rapid  identification  of 
atmospheric  contaminants/signatures  in  window  regions  as  well  as  accurate  spectral  calculations  of 
heating/cooling  rates  in  the  presence  of  clouds  for  both  thermal  and  solar  spectral  regimes.  Finally,  in  the 
regions  of  molecular  opacity,  where  weighting  functions  peak  in  the  atmosphere,  it  is  expected  tiiat 
MODTRAN4  can  play  a  role  in  very  quick  ‘primitive’  retrievals,  avoiding  the  large  number  of  line-by-line 
(LBL)  calculations  necessary  for  initiating  derivative  (perturbation)  matrices  (Anderson  et  al.,  1993). 
While  the  error  estimates  and  residuals  associated  with  a  2  cm’*  algorithm  will  be  larger  than  those 
associated  with  LBL  retrievals,  the  speed  advantage  for  image  processing  might  warrant  this  iiutial 
approach. 


FASE  (SneU  et  al.,  1995  a,b)  can  now  be  used  for  thermal  infiured  hyperspectral  sensing,  most  notably 
for  atmospheric  compensation  for  interferometer  data  taken  firom  aircraft  or  satellite  platforms.  Based  on  the 
validation  of  the  HIS  (High  resolution  Interferometer  Sounder),  Wang  et  al,  1996,  FASCOD3  proved  capable  of 
fitting  the  measurements  within  a  few  percent,  given  accurate  sonde  data.  But,  this  paper  equaUy  showed  that  the 
high  resolution  measurements  can  actuaUy  reveal  the  state  of  the  atmosphere,  independent  of  sonde  (Revercomb  et 
al.,  1993;  Smifli  et  al.,  1995).  This  process  is  termed  ‘atmospheric  compensation’  when  applied  to  hyperspectral 
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imaging.  HIS  is,  of  course,  a  ‘soda  straw’  rather  than  an  imager,  but  the  inferences  can  be  extrapolated  to  the  next 
generation  of  instruments.  Certainly  surfece  properties  can  be  ascertained  (^.  the  sea  surface  albedo  in  the  Wang 
paper),  and  sea  surfece  temperature  can  be  readily  retrieved. 


2.  ADDITION  OF  A  CORRELATED-it  CAPABILITY  TO  MODTRAN 

Addition  of  a  CK  capability  to  MODTRAN  (Bernstein  et  al.,  1995)  provides  an  accurate  and 
fast  means  for  evaluation  of  the  effects  of  clouds  and  heavy  aerosol  loading  on  retrievals  (both  surface 
properties  and  species  concentration  profiles)  and  on  atmospheric  radiative  heating/cooling  calculations. 
These  radiative  transfer  computations  require  coupling  the  effects  of  gaseous  molecular  absorption  due 
primarily  to  water  vapor,  carbon  dioxide,  and  ozone,  with  particulate  multiple  scattering  due  to  volcanic 
aerosols,  ice  crystals,  and  water  droplets.  In  order  to  adapt  a  band  model  approach  for  use  in  scattering 
calculations  it  is  necessary  to  express  the  band  model  transmission  function  in  terms  of  a  weighted  sum  of 
Beer's  law  exponential  terms.  Thus,  a  method  for  determining  the  weighing  factors  and  monochromatic 
absorption  coefficients  for  the  MODTRAN  band  model  is  required.  An  abbreviated  discussion  of  the  CK 
approach  as  tailored  for  integration  into  MODTRAN  is  given  below;  for  a  more  complete  discussion  of  the 
CK  method  the  reader  is  referred  to  Lacis  and  Oinas  (1991). 

For  simplicity,  consider  the  problem  of  determining  the  average  transmittance,  as  defined  by 
Beer's  law,  for  a  homogeneous  path  over  a  finite  spectral  interval.  The  generalization  to  inhomogeneous 
paths  is  straightforward.  The  path  transmittance  can  be  exactly  determined  through  evaluation  of: 

j  02 

T(u)  =  - f  dof  exp(-k(co)u)  , 

0)2-0), 

where  ©  is  frequency,  k(o)  is  the  monochromatic  absorption  coefficient,  and  u  is  absorber  column  density. 
The  basis  of  the  CK  approach  is  that  evaluation  of  T(u)  by  integration  over  frequency  can  be  replaced  by 
an  equivalent  integration  over  the  distribution  of  absorption  coefficient  values  f(k)  in  the  spectral  interval 

CO 

T(u)  =  \dkf(k)  Qxp(-ku)  . 

0 

The  distribution  function  f(k)  is  not  smooth  or  monotonic;  it  generally  consists  of  a  series  of  sharp 
spikes  which  reflects  the  sharp  line  structure  of  k(©).  It  then  becomes  more  computationally  convenient  to 
work  with  the  smooth  and  monotonic  cumulative  probability  distribution  function 

g(k)  -  ]d'kf(k')  . 

0 

where  k(g)  is  given  by  the  inverse  of  g(k),  k(g)  =  g'‘(k). 

The  MODTRAN  band  model  for  a  single  species  is  based  on  four  parameters:  (1)  the 
integrated  line  strength  S  in  a  spectral  interval  A©  (A©  =  1  cm"'  in  MODTRAN),  (2)  Ae  effective  number 
of  equivalent  lines  n  (non-integer  values  of  n  are  acceptable)  in  the  interval,  (3)  the  average  pressure 
broadening  Lorentz  line  width  Yl,  and  (4)  the  Doppler  line  width  yp.  These  parameters  are  determined 
directly  from  the  1996  HITRAN  parameter  line  compilation  (Rothman  et  al.,  1996). 

3.  INITIAL  VALIDATION  OF  MODTRAN4  MULTIPLE  SCATTER  SOL/^ 

Development  of  MODTRAN4  was  necessary  because  MODTRAN3  predictions  of  multiple 
scatter  solar  radiances  in  spectral  regions  where  non-continuum  molecular  absorption  is  important  are 
inaccurate.  This  scenario  cannot  be  validated  by  comparisons  to  FASCODE  (Clough  et  al.,  1988)  because 
of  the  lack  of  solar  scattered  contributions.  Thus,  initial  validation  of  MODTRAN4  multiple  scatter  solar 
calculations  were  made  directly  to  measurements. 
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Figure  1:  A  comparison  between 
SICAP  measurements 
and  model  predictions  for  a  solar 
illuminated  cloud  top.  See  text 


For  example.  Figure  1  illustrates  a  comparison  of  calculated  radiance  predictions  to  airborne  measurements 
(Malherbe  et  al.,  1995)  performed  by  ONERA  and  CELAR  using  SICAP,  a  circular  variable  filter 
cryogenic  spectrometer  (1500  -  5500  nm,  2%  spectral  resolution).  The  aircraft  altitude  was  3.0  km,  the 
cloud  top  was  2.5  km,  the  sensor  line-of-sight  (LOS)  zenith  angle  was  104*,  and  the  solar  zenith  and 
relative  azimuth  angles  were  48"  and  137",  respectively.  Three  sets  of  calculations  are  shown:  (1)  results 
from  NAULUM  (Malherbe  et  al.,  1995),  a  new  radiative  transport  model  developed  at  ONERA,  (2) 
MODTRAN3  calculations  without  the  CK  approach,  and  (3)  MODTRAN4  calculations  performed  with  Ae 
CK  approach.  The  MODTRAN4  cloud  model  upgrade  enabled  the  cloud  profile  and  spectral  data  to  be 
explicitly  entered  for  both  the  MODTRAN3  and  MODTRAN4  calculations.  Both  MODTRAN  calculations 
were  performed  using  a  simpler  two-stream  multiple  scattering  model  (Isaacs  et  al.,  1987);  the  discrete 
ordinate  model  in  MODTRAN,  DISORT  (Stamnes  et  al.,  1988),  was  run  with  8-streanis  over  a  limited 
spectral  sub-region  and  product  similar  results.  The  fit  is  excellent,  but  if  the  sun  is  in  the  near  line  of 
sight  (eg.  1 1"),  the  fit  would  not  properly  capture  the  forward  scattering  lobe. 


4.  COOLING  RATES  IN  THE  PRESENCE  OF  CLOUDS 

An  important  part  of  model  development  is  continual  comparison  with  measurements,  and 
theoretical  predictions.  In  recent  years  the  ARM  community  participants,  particularly  Clough  and 
colleagues,  have  undertaken  the  extension  of  the  work  originally  begun  by  F.  Lu&er  and  R.  Ellingson,  as 
represented  in  the  dedicated  issue  of  J.  Geophys.  Res.,  96,  May  1991.  The  cooling  rate  discussions 
therein  have  been  augmented  by  more  recent  studies  by  Clough  et  al.,  1992,  Bernstein  et  al.,  1996,  and 
Kimball,  1992,  for  instance.  In  Figure  2a  the  cooling  due  to  HjO  alone,  in  the  0-3000  cm-1  spectral  range 
is  represented,  including  a  0.5  km  thick  cirrus  cloud  centered  at  6.5  km.  The  atmosphere  is  sub-arctic 
winter;  integrated  results  can  be  seen  in  2b.  This  result  is  preliminary,  but  MODTRAN4  now  allows  for 
the  inclusion  of  such  clouds  with  multiple  scattering  and  micro-layering.  Initial  results  show  a  clear 
dependence  on  the  mode  of  implementing  the  cloud  within  the  code,  with  as  much  as  a  0.7K/day  difference 
(MOD3  vs.  MOD4)  from  that  shown  if  2b,  as  well  as  0.5K/day  difference  between  separate  MOD4 
calculations.  Certainly  much  additional  theoretical  work  needs  to  be  done,  but  without  good 
measurements,  the  impact  of  cloud  cooling  will  remain  uncertain. 
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Figure  2:  Spectral  and  integrated  cloud  cooling  for  a  SAW  atmosphere,  with  an  embedded  thin  cirrus. 
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